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1. Introduction 



Since its discovery in 1998 Jl]||, dark energy has become one of the central problems 
in theoretical physics and cosmology [151J^,^,|6|J^,^,|9|, |l^|Tl1 , [T^ , |l^ , |l4| . Thousands of papers 
have been written on this subject, while it is still as dunting as ever to understand the 
nature of dark energy. 

We start with a brief history of the dark energy problem in time order [|3|,|T5|1 . 



1917: Einstein added a cosmological constant term in his field equations, for the following 
reasons: Firstly, for isolated mass not to impose a structure on space at infinity in closed 



universe; secondly to obtain a static universe [16 



1920s: Pauli realized that for a radiation field the vacuum energy is too large to gravitate. 
As we shall review, the vacuum energy density of a radiation field is proportional to the 
cutoff to the fourth power. Pauli (unpublished) showed that using the classical electron 
radius as an ultraviolet cutoff, the curvature of the universe will be so large that the 
universe "could not even reach to the moon" . 

1931: Einstein removed the cosmological constant |17j] because of the discovery of the 
cosmic expansion. Although in 1923 on a postcard to Weyl he already wrote: "If there is 
no quasi-static world, then away with the cosmological term!" , it seems that he does not 
believe in "no quasi-static world" until [17]. 

1960s: To explain why there are so many quasars centering around red-shift around z = 
1.95, some people suggested to use the Lamaitre model (with A > 0, k = 1), and around 
a = 1/(1 + z) = 1/2.95, the universe is like the Einstein's static universe ||18|| . 
1967: Zeldovich reintroduced the cosmological constant problem by taking the vacuum 
fluctuations into account. This introduced the old cosmological constant problem, Zel- 



dovich used the word "fine-tuning" [19,20]. 



1987: Weinberg "predicted" a non- vanishing and small cosmological constant [^T|, and 
two years later published by now the famous review article 0. 

1998: Based on the analysis of 16 distant and 34 nearby supernovae, Riess el al. first 
discovered the acceleration of expanding universe |IJ. Soon after, utilizing 18 nearby 
supernovae from the Calan-Tololo sample and 42 high-redshift supernovae, Perlmutter et 
al. confirmed the discovery of cosmic acceleration @. 

2000s: String theorists reintroduced the anthropic principle when discovered the string 
landscape []22],|3 . 



2011: Because of the great discovery of cosmic acceleration, Adam Riess, Brian Schmidt, 
and Saul Perlmutter win the Nobel prize in physics 2011. 

We will review a set of theoretical ideas and phenomenological models of dark energy, 
we will also review some numerical works done so far. Thus the body of this review 
consists of two parts, the first part reviews some of major theoretical ideas and models of 
dark energy, and the second part is devoted to reviewing some numerical works. 

Summary of review on theoretical work: We will review briefly the history of the 
problem of the cosmological constant. After this, we follow Weinberg's classical review to 
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divide the old approaches reviewed by Weinberg into five categories. Then, we add three 
more categories in order to include most of the more recent ideas and models. 

Summary of review on numerical work: First, we will review the mainstream cos- 
mological observations probing dark energy, including the type la supernovae (SNIa), the 
cosmic microwave background (CMB), the baryon acoustic oscillations (BAO), the weak 
lensing (WL), the galaxy clusters (CL), the gamma-ray burst (GRB), the X-ray observa- 
tions, the Hubble parameter measurements, the cosmic age test, the growth factor (GF), 
and so on. We will introduce the basic principles of these observations, describe how these 
observations are included into the x 2 statistics, and introduce the related research works 
of the observational groups. 

Next, we will provide a brief overview of the present and future dark energy projects. 



According to the Dark Energy Task Force (DETF) report ||24|| . the dark energy projects 
can be classified into four stages: completed projects are Stage I; on-going projects, ei- 
ther taking data or soon to be taking data, are Stage II; near-future, intermediate-scale 
projects are Stage III; larger-scale, longer-term future projects are Stage IV. Some most 
representative dark energy projects of Stage II, Stage III, and Stage IV will be introduced 
in this work. 

Lastly, we will review the current numerical studies on dark energy. Faced with the 
plethora of DE scenarios, the cosmologists currently have two options: (1) they can test 
every single theoretical model against observations and make a comparison of different 
dark energy models; (2) they can reconstruct the evolution of dark energy based on the 
observational data and explore the nature of dark energy in a model independent manner. 
The related research works of these two routes will be introduced in section 15 and section 
16, respectively. 

The conventions used in this review are as follows: we use the metric convention 
(—,+,+,+), and use natural units c = U = 1. The Newton's constant and the reduced 
Planck mass are kept explicit and they are related by M 2 = 1/(8ttG). 

2. The theoretical challenge 

Einstein was the first to introduce the famous cosmological constant term in his equa- 
tions. In the Einstein equations 

R tl „-^g^R = 8irGT flv , (2.1) 



one adds a term to the stress tensor on the R.H.S. of the above equations 

Tfj.u T[iv - ^q A 9^, (2-2) 

where A is a constant, the cosmological constant. The Einstein equations can be rewritten 

as 

R„u ~ \g^R + A#m- = SttGT^. (2.3) 

Before trying to understand the nature of A, we note that it has a couple of ready 
interpretations. First, if we take the additional term in ( |2.2| ) as coming from some ideal 
fluid, whose stress tensor is given by 

Tpv = (pa+Pa)uhU»+Pa9hu, (2-4) 

where p\ is the energy density and p\ the pressure, then the cosmological constant can be 
interpreted as a fluid with 

1 . . 

Pa = -Pa, Pa = -^q A - (2-5) 

This is certainly an unusual fluid, since if A > 0, the pressure is negative, and the strong 
energy condition is violated since (T^ u — ^Tg lJLV )u^u 1 ' = —A < 0. Of course the strong 
energy condition is not something sacred. The null energy condition is marginally satisfied. 

Second, since the cosmological constant term is proportional to g^ v thus is Lorentz 
invariant, it can be interpreted as the vacuum energy. And indeed Lorentz invariance 
forces upon us the condition p = —p. In principle, before we understand the origin of the 
vacuum energy, the energy density can be positive, negative and zero. 

In the Friedmann-Robertson- Walker cosmology, the two Friedmann equations read 

3MpH = p m + pa !-, 



6M 2 - = 2p A -p r 
y a 



a2 (2.6) 



where = 1/(SttG) and M p is usually called the reduced Planck mass, H = a/a is the 
Hubble "constant", p\ is the vacuum energy density p\ = M^A, and p m is the matter 
energy density (with p m = 0) . Einstein introduced a positive cosmological constant moti- 
vated by a static universe. To have a static universe H = and without the cosmological 
constant, we deduce from the first equation above p m = 3Mp/a 2 for a positive spatial 
curvature k = 1. But this does not work, since the second Friedmann equation we have 
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d/a < for p m . Namely the universe will collapse due the attractive force of matter. Thus, 
to have both H = and d/a = we need to have 2p A = Pm and pa = M 2 /a 2 . This is 
why the cosmological constant has to be positive and fine-tuned to balance matter. The 
spatial curvature also has to be positive. 

However, due to the cosmic expansion, Einstein himself later abandoned the cosmo- 
logical constant. Nevertheless, as the cosmological constant is the simplest extension of 
the original Einstein equations, the theoretical possibility is left open. 

Up to now, we have only discussed the classical story. Without quantum mechan- 
ics, a very small cosomological constant poses no problem, one simply regards it as a 
parameter in theory. However, in quantum mechanics, we know that vacuum fluctuations 
make contribution to the energy of a vacuum, thus the vacuum energy density receives 
two contributions, one may be called the bare vacuum energy, the classical one without 
quantum contribution, the other comes from the zero-point fluctuations of all quantum 
fields (in quantum field theory), and the zero-point fluctuations whatever dynamic degrees 
of freedom there are. 

For a quantum field with a given mode of frequency u, the zero-point energy is ±^cj, 
with the plus sign for a bosonic field and the minus sign for a ferminic field. The total 
zero-point energy is then | ^2i(±)<^i- ^ n the continuum limit, we have for a free field 

i 

this integral is divergent, thus we need to introduce a physical cutoff. A cutoff is reasonable, 
for example, if the zero-point energy density is infinite, then our universe is infinitely curved 
and the space size is infinitely small. Thus by examining the Friedmann equations, we know 
that this cutoff must be the Planck energy at largest. 

Let the cutoff be A > m, the energy density of a bosonic field is then 

f x k 2 dk ri2 2 ,i f x k 3 dk A 4 , nn . 
Take, for example, A = M p , then 

1 



P = 



2 x 10 71 GeV 4 . 



2 10 7T 4 G' 2 

But in reality, p A < p c « (3 x 10" 12 GeV) 4 « 10" 46 GeV 4 . If there is no supersymmetry, 
| Pa | should be greater than p^ from a single bosonic field, namely pb < p c , but what we 
have is 

^ « 10 117 . (2.9) 

Pc 

6 



This is the famous cosmological constant problem. Because if the cosmological constant is 
much larger than p c , the universe will never look like what we observe today. 

With exact supersymmetry, the bosonic contribution to cosmological constant p\, is 
canceled by its fermionic counterpart. However, we know that our world looks not super- 
symmetric. Supersymmetry, if exists, has to be broken above or around lOOGeV scale. 
Even if we take the Asm =100GeV cutoff, below which physics is believed to be well- 
described by the particle physics standard model, still 

pb ~ Mm = i 6 G eV 4 , ^ w 10 52 . (2.10) 

167T 2 p c 

One could use fine-tuning to solve the problem in some sense, by introducing a bare 
cosmological constant and letting it cancel with the quantum contribution using renormal- 
ization. However, one has to make two independent numbers cancel by the accuracy of 
one part in 10 117 or at least 10 52 . This is extremely unlikely to happen. Thus it remains a 
problem why the cosmological constant is not large. This is known as the old cosmological 
constant problem ||19j| . 

The cosmological constant problem remained the above statement until crucial exper- 
iments came in. Riess et al. and Perlmutter et al. in 1998 discovered the accelerating 
expansion of our universe Jl]|| • The simplest explanation of this phenomenon is the return 
of a positive cosmological constant. Let p c = 3MpH 2 be the critical energy density, data 
available today tell us that pa = 0.73p c and p m = 0.27 p c . Thus not only there is a small 
positive cosmological constant, but also its energy density is the same order of matter en- 
ergy density. These discoveries lead to a new version of the cosmological constant problem, 
or the problem of dark energy. 

The modern version of the cosmological problem divides itself into two parts: 

(a) Why pa ~ 0, namely why it is so small? 

(b) Why pa ~ Pm, this is called the coincidence problem. 

Our writing of the theoretical part of this review is motivated by Weinberg's review 
article [Q, we will first review the theoretical efforts before the discovery of acceleration 

0. 0. following Weinberg's classification: 

1. Supersymmetry and superstring. 

2. The anthropic principle. 

3. The tuning mechanisms. 

4. Modifying general relativity. 
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5. Ideas of quantum gravity. 

Twenty years have passed since Weinberg's classification, all the new ideas more or less 
still belong to the above five categories. For example, category 1 and category 2 combine 
to become the string landscape + the anthropic principle, and there are new ideas in the 
category of symmetries. The tuning mechanisms now include ideas associated with brane- 
world models. The category of modifying gravity now includes f(R) models, f(T) models, 
MOND and TeVes models and the Dvali-Gabadadze-Porrati (DGP) model. There are also 
some progress in the category of quantum gravity. 

Further, we can add three more categories: 6. The holographic principle; 7. Back- 
reaction of gravity; 8. Phenomenological models. Thus, we will also write about: 

1. Symmetry. 

2. The anthropic principle. 

3. The tuning mechanisms. 

4. Modifying gravity. 

5. Quantum gravity 

6. The holographic principle 

7. Back-reaction of gravity. 

8. Phenomenological models. 



3. Weinberg's classification 

In this section, we will briefly review Weinberg's classification of ideas about the 
cosmological constant problem as given in Q. 

3.1. Super symmetry and superstring 

In any supersymmetric theory in 4 dimensions, the supersymmetry algebra contains 
at least J\f = 1 generators Q a and their conjugate Q^ a [^5| such that 



{Qa,Q^} = K) Q/3 P^, (3.1) 

where Oi are the Pauli matrices and ao = 1. If supersymmetry is unbroken, then 



Q a \Q) = QUQ)=0, 
8 



(3.2) 



the supersymmetry algebra leads to 



H\Q) = P°\n) = 0, (3.3) 



namely, the unbroken vacuum has exactly zero energy |26[ . 

In a supersymmetric quantum field theory, there are a number of chiral super-fields 
for which the potential V is determined by the super-potential W . W is a function of 
complex scalar fields cf) 1 in the chiral multiplets, 

V = J2\ 9 iW\ 2 , (3.4) 

where djW = dW/d<p l . For a vacuum with unbroken supersymmetry V = thus d^W = 0. 

Of course any supersymmetry must be broken in our world, thus in general 
YliaiQoLiQ'l,} = 2if > 0. If there is translational symmetry in space-time, we must 
have H ~ Vp\ thus pa > 0. So a positive cosmological constant is a consequence of 
broken supersymmetry in quantum field theory, as long as this is the whole contribu- 
tion to the effective cosmological constant. This conclusion agrees with observations since 
1998, but what is the exact value of p\! For a typical quantum field theory, one ex- 
pects pa ~ Mg USY 3> p c . Thus we would say that supersymmetry does not solve the 
cosmological constant problem. 

It was commonly thought that it is better to explain pa = first, then take the next 
step to explain why pa ~ p c - 

The potential in a super-gravity theory is determined by both the super-potential and 
the Kahler potential Kty^fi) 



V = e K 

where we already set 8ttG — 1, and 



G lJ DiWDj(l>-3\W\ 



(3.5) 



DiW^diW + diKW, G i3 = d4<dpW, G^G- jk = 8l (3.6) 

In a class of the so-called no-scale super-gravity models [^S[ , one can fine-tune param- 
eters to break supersymmetry, meanwhile keep the vacuum energy vanishing. There are 
three classes of complex scalar fields, C a , S n and T. Let 

K=-3ln(T + f-h(C,C)) + K(S,S), W = Wx(C) + W 2 (S), (3.7) 

9 



then 



V = e K (T + T -h)- 3 N ah d a Wd b W + G mn D m WD n W , 



(3.8) 



where iV = d a dih. Since both N and G are positive definite, so V > 0. To have V = 
there ought to be 



if W 0, then D a W ^ and by definition supersymmetry is broken. 

In such a model, T is not fixed, this is why a model like this is called no-scale super- 
gravity model. The problems for this model include: 

• The coefficient of the first term must be -3, this is fine-tuning. 

• The form of W is fine-tuned. 

• Quantum corrections usually spoil those fined-tuned coefficients. 

The no-scale super-gravity models continue to attract attention today. 

3.2. Anthropic principle 

The terminology of the anthropic principle is due to Brandon Carter, who articulated 
the anthropic principle in reaction to the Copernican Principle, which states that humans 
do not occupy a privileged position in the Universe. Carter said: "Although our situation 
is not necessarily central, it is inevitably privileged to some extent." |[29|| . 

As Weinberg formulated in his review H, there are three different kinds of anthropic 
principle: Very strong version, very weak version and weak version. 

The very strong version states that everything in our universe has something to do 
with humankind, this is of course absurd. 

The very weak version takes the very existence of our humankind as a piece of exper- 
imental data. For instance, in order not to kill a person with the products of radio-decay, 
the life-time of a proton must be at least 10 16 years. 

Now the weak version. This is the version that prevails in certain circle of people. In 
this version, it is postulated that there are many regions in the universe. In these regions 
physical laws are in different forms. It just happens that in the region we are dwelling 
all physical laws, physical constants and cosmological parameters are such that clusters of 
galaxies, galaxies and our solar system can form, and humankind can appear. It appears 



d a W = D m W = 0. 



(3.9) 



The point where V = is a minimum, so there is no instability problem. For 



D a W = d a W + d a KW = d a KW ~ d a hW, 



(3.10) 
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that all the conditions are fine-tuned but they are not, because all different sorts of regions 
exist in the universe, such a universe is called multiverse. We found ourselves in our part 
of the multiverse simply because this is all we can observe. 

In such a weak version, we are not supposed to expect every physical law and every 
physical constant be designed for the existence of humankind. For instance, we do not 
know whether the fact that the life-time of proton is more than 10 33 years has anything to 
do with us, or whether it is connected with other conditions which are necessary for the 
existence of us. 

Dicke |3(],|31| may be one of the first persons making use of this weak version of the 
anthropic principle, when he was considering Dirac's large number problem. Let T be the 
age of the universe. The death date of the sun is greater than T, and T is greater than 
the formation time of the second and the third generations of stars. Thus T is about 10 10 
years. This is a nice explanation of a large number. However, other time scales such as 
l/(m 3 /p) ~ 10 10 years are completely incomprehensible. 

Weinberg's anthropic consideration of the cosmological constant PI[] is considered to 
be a genuine prediction of the anthropic principle by some people. The simplified version of 
his argument is the following. Let z be the red-shift when galaxies form, the matter density 
is Pm(z) = (1 + z) 3 p Q m ~ 100p^. The formation of galaxies crudely requires p\ < 100p^. 

It is assumed that the primordial density fluctuation is about dp/p ~ 10~ 5 in Wein- 
berg's calculation. But in a full anthropic calculation, even this number is not to be 
presumed j32[] (for a complete discussion on "scanning parameters", see H33j ). If we let 
Sp/p be a free parameter, we usually have p\ < Xp Q m and X 3> 100!, thus the currently 
observed value of dark energy can not be considered a consequence of the anthropic argu- 
ment. 

Weinberg also pointed that the age of the universe is a problem. He argues that 
if Pa/Pw = 9 then T = I.IHq . Actually, other people, including de Vauconleurs |3jJ] , 
Peebles fl35| , Turner, Steigman and Krauss |3f| also considered this issue. 

The anthropic argument is also applied to a possible negative cosmological constant. 
Let pa < 0, then the Friedmann equation 3M 2 H 2 = p m + Pa implies that when p m = \pa\, 
the universe starts to contract, thus \p\\ < p® m . 

11 



3.3. Tuning mechanism 

The idea is to use a scalar to self-tune the stress tensor, which is a function of this 



scalar ||37|| . Assuming 

VV^ ~ T% ~ i?, (3.11) 

where i? is the scalar curvature. If the trace of the stress tensor vanishes when <fi rolls to a 
certain value <po, <p will stay at this value and the effective cosmo logical constant becomes 
vanishing. However, it can be proven that 

Tftv = e^g^Co (other fields), (3.12) 

so this tuning mechanism can not be realized (otherwise this will force Co = 0). 

Another possibility is <fio = — oo, but in this case the effective Newtonian constant 
vanishes. 

3-4- Modifying gravity 

This is a rather popular theme now, but at the time when Weinberg wrote his review, 
there was only one proposal mentioned. This is the uni-modular metric theory |3"8| . |3l| (see 
5[],^T| for recent progress). 



The idea is very simple. Removing the trace part from the Einstein equations 

Rnv - ^g^R = StvGT^, (3.13) 

we will have 

iV " ~ A 9^R = 8kG(T^ - ^g^T), (3.14) 
where T = Tfi. Assume be conserved, namely V^T^j, = 0, we deduce from ( |3.14| ) that 

d^R = SirGdfjT, (3.15) 

thus 

8nGT = -R + 4A, (3.16) 



where A is an integral constant, when we use (|3.14j ) instead of (|3.13|) as a starting point. 



Substituting this back to the traceless Einstein equations, we obtain 



Rfiu - ^g^R = SnGT^ - Ag^, (3.17) 
12 



so A is indeed the cosmological constant. That A is an integral constant is due to the fact 
that there are one fewer equations in the traceless Einstein equations. In fact Einstein 
himself considered this theory for a while. 

The traceless equations can be considered as a consequence of the requirement that 
the determinant of the metric is 1, so we have g^ v 5g^ v = or 5^/g = 0, we can 
introduce a Lagrangian multiplier into the action to enforce this. 

As an integral constant, now A is considered as a free parameter in the theory to be 
determined by initial conditions, or to set a framework for utilizing the anthropic principle. 

3.5. Quantum cosmology 

An accurate definition of quantum cosmology does not exist since there is no theory 
of quantum gravity to provide an appropriate framework. 

The so-called quantum cosmology as advocated by Hawking can at most be consid- 
ered as a qualitative picture of more accurate underlying theory. One starts with the 
Hamiltonian constraint 

(3) i? - 2A + N- 2 {E l3 E t:i - E 2 ) - 2N- 2 T 00 = 0, (3.18) 

where N is the lapse function as in goo = —N 2 1 and E l i = 5/5hij, is the scalar 
curvature of the three spatial geometry in the ADM splitting. The wave equation governing 



the wave-function of the universe is called the Wheeler-De Witt equation [f42|,pl3 

(- &i,lh jm& +(3) R ~ 2A " 2To °) * = °' (3 ' 19) 

where ^ is the wave function of the universe, it is a functional of the three geometry 
and other fields on the three dimensional spatial slice, however, it contains no time, and 
time has no place in the constraints. Actually, the Hamiltonian constraint is a consequence 
of the requirement of time-reparametrization invariance. As a consequence, we can not 
impose the usual normalization condition on \&: 

J |*(%,</>)| 2 [dhd(f)} = 1. (3.20) 

Weinberg proposed to use one of the dynamical variables buried in hij or <fi to replace 
the role of time. 
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There are infinitely many solutions to the Wheeler-De Witt equation, if we assume 
this equation be well-defined. Hartle and Hawking proposed in |44]] to select one out many 
by using the path integral 

*(h ij ,<f>) = J[dgd4>]Me- SE , (3.21) 

where we assume that the three dimensional space £ is the boundary of the four dimen- 
sional space M and g is a Euclidean metric on M, Se is the Hilbert-Einstein action on 
this Euclidean four space obtained by an prescription of Wick-rotation: 

Se = f V9( R + 2 A) + S B (<f>), (3.22) 

where Se{4>) is the Euclidean action of the matter fields 0. Since the path integral is 
determined by the value of and ^ on S, there are no additional boundary conditions. 
This is called by Hartle and Hawking the no-boundary proposal of the wave function of 
the universe. 

To have a probabilistic interpretation, consider a physical quantity A, a function of 
hij and ^ on S. We define the probability of A assuming value Aq be 

P(A ) = J [dgd(/>)S(A{h^) - A )e- SE . (3.23) 

However, A is a fixed parameter in action Se- In order to compute the probability distri- 
bution of A, we need to make it variable. Hawking introduced a four form field strength 
in [15] to make the effective A a dynamic variable (see also [HB1) 



F^ Xp = 4d {ll A uXp] , (3.24) 

with the action 



S(A) 



^ J V^F^ Xp F^. (3.25) 



2 x 

The equation of motion for F is 

d^ vXp = 0. (3.26) 

Since F is totally asymmetric, let F^ uXp = Fe^ vXp j y/—g, the equation of motion leads to 
d^F = 0, namely F =constant, and 

S(A) = lf V^F 2 , (3.27) 



2 _ 
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compared with the action 



S(g, A) = — I V^(R ~ 2A), (3.28) 



we may conclude that A — > A — 4nGF 2 . This is incorrect, since we need to consider the 
energy of this solution E = f d 3 x^F 2 , and this leads to 

A — >• A + 4nGF 2 . (3.29) 

In the Euclidean action 

Se = T^gJ ( R + 2A )y/g + 2^41 / V9F^XpF^, (3.30) 
let F^ Xp = Fe^P/y/g, we have 

Se = J( R + 2 ( A + 47tGF 2 ))v^, (3.31) 

this agrees with the consideration of energy. 



Now, according to (j3.23|) , we compute the probability distribution of A: 



P(Ao) = J [dgdA]d(A - Ao)e" SB , (3.32) 

where A = 4-kGF 2 (for simplicity we let the bare cosmological constant be zero). The 
most contribution to the above path-integral comes from a classical solution to equation 
3M 2 H 2 = \F 2 . We find 

Se = -Q^py ( 3 -33) 

thus 

P(A )~e^. (3.34) 

Hawking concluded that Aq — > +0 is the most probable value. But we now know that this 
is a wrong prediction! 

Coleman later pointed out that Hawking's consideration is not the complete story. He 
suggests that one needs to take wormholes and baby universes into account [[|7| . Each type 
of baby universe is characterized by its physics properties, and we use % to label them. Let 
a\ be the creation operator for the baby universe of type i. Local creation of wormholes 
has the effect of modifying the action: 

S -> S = S + ^(a, + at) J Oi(x), (3.35) 
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(3.36) 



this effective action works for the parent universe, and Oi is an local operator. 
Consider a state without any wormhole, so ai\B) = and 

\B) = f Hdaifia^a), 

i 

(di + a\)\a) = ati\a}, 

and 

i 

When acting on \a), the action becomes 

S^S + J2 a i J Oi, (3.38) 

i 

thus Qij becomes a coupling constant in the parent universe. The cosmological constant can 
be regarded as a coupling constant, its corresponding operator is the lowest dimensional 
operator 0\ = \f—g. In this framework, coupling constants including the cosmological 
constant become dynamical automatically. No additional mechanisms such as four form 
fluxes are needed. When observation is made on spacetime geometry, the observers will find 
the wave function of the universe to be in an eigenstate of \a). Coleman suggested that any 
apparently disconnect universes could be actually connected by wormholes. In this sense 
all the other apparently disconnected universes, which Hartle and Hawking ignore, should 
be summed over. This argument will not affect real constants in nature but does affect 
the effective "constants" which come from the baby universe creation operators ( |3.35|) . 
Using techniques well developed in quantum field theory (the summation of all vacuum to 
vacuum Feynman diagrams is the exponential of the summation of connected diagrams), 
one obtains 

P(a) = exp (^J [dgdA]e~ SE ^j , (3.39) 



P(A)=exp(exp(|^)), (3.40) 



and it leads to 

P( \) = fixn I fixnf 

GA 

again this predicts A = 0, a wrong prediction. 

We have briefly reviewed Weinberg's classification, now we turn to more recent ideas 
and models about dark energy. 
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4. Symmetry 

4-1- Supersymmetry in 2+1 dimensions 

Witten pointed out in 1995 that supersymmetry in 2+1 dimensions may help to solve 



the zero cosmological constant problem [ 48|j . If such a theory exhibits the same phe- 



nomenon as the type IIA string theory in 10 dimensions, an additional dimension may 
emerge and becomes the third spatial dimension. The 2+1 dimensional supersymmetry is 
smaller than the 3+1 dimensional supersymmetry. As a result the unbroken 2+1 dimen- 
sional supersymmetry (in 3+1 spacetime with an emergent spatial dimension) does not 
contradict with observations and meanwhile still forces the cosmological constant be zero. 

In 2+1 dimensions, one can argue that there is no boson-fermion degeneracy. The 
existence of a particle of mass m creates a deficit angle in 2 spatial dimensions, 9 ~ mGs, 
where G3 is the 3 dimensional Newtonian constant. The appearance of the deficit angle 
makes the definition of charge impossible, thus the energy degeneracy becomes impossible 
too. 

Becker, Becker and Strominger used an Abelian Higgs model to realize Witten's idea 
49(1 , but unfortunately there is no emergent third dimension in their model. Their model 



is a supergravity theory with field content (</>, A M , N, x, A) + (g^, ij)^). BBS showed that 
the spectrum of vortices in this theory does not exhibit boson-fermion degeneracy. 

4-2. 't Hooft-Nobbenhuis symmetry 

't Hooft and later 't Hooft and Nobbenhuis proposed to consider symmetry transfor- 



mation 50 



x n iyii^ p v -iVy- (4.1) 
As the simplest example, Consider a scalar field 

S = f d*x(-±{d<f>)*-V{<f>)), (4.2) 



2 

with the stress tensor 

T^ = d^0 + ^£(0). (4.3) 
Under transformation (O), or = iy^, = id^, 



c y = -c = -l(dy<p) 2 + v, 

TjL = = dim + g^uC 
17 



(4.4) 



We have in particular Tq = — T 00: however since d D 1 x = i D 1 d D 1 y, the Hamiltonian 
is not simply reversed in sign 



Let 



H* = J I&d^y = -(-i^H. (4.5) 

= / d D ~ l p {a(p)e^ + a\p)e-^) , 

^(^ = y d °~^P (-ia(p)e ipx + ia\p)e~ ipx ) , M 
P° = (P 2 + m 2 ) * ■ 



We have 



<f>(iy t it) = J d^qiay^e^ + a^e-^), 

n y = in(iy, it) = J dP^q {-ia y (q)e iqx + ia\{q)e- ipx ) 



(4.7) 



and 



<f> = (<t > -m 2 )*, a y (q) = (-i) D - 1 a(p), a\{q) = (-i) ' 1 ^ (p), (4.8) 

if a^(p) is the Hermitian conjugate of a(p), al is longer the Hermitian conjugate of a y . 

Note that if we demand T^ u — y —T^ u be a symmetry, then T o|fi) = 0. Upon intro- 
ducing gravity, let g^ v = g^(x = iy) thus ds 2 x = -ds 2 y , ->■ —R^v, or R^ = -R* u (iy). 
Start with the Einstein equations with a cosmological constant 

Rfiu - \g^R + A<7^ = SttGT^, (4.9) 

we obtain 

R% ~ \al„R v - A</», = SttGT^, (4.10) 

where we used T^ v — y —T^. Demanding \Q) be invariant, we deduce A = 0, This trans- 
formation maps a de Sitter space to an anti-de Sitter space. 

't Hooft and Nobbenhuis pointed out that a scalar field and an Abelian gauge field 
can realize this symmetry transformation but 

• The delta function 5 3 (y) need be treated carefully in the second quantization scheme. 

• m 2 — y — m 2 , leading to tachyon. 

• This symmetry can not be realized in a non-abelian gauge theory. 

• The boundary conditions need be treated carefully. For instance, the boundary conditions 
at x = oo makes H y < in quantum mechanics, and makes it equal to —iH x in quantum 
field theory. 

This concludes our discussion on the 't Hooft-Nobbenhuis symmetry. 



18 



4-3. K apian- Sundrum symmetry 

This symmetry is quite similar to the 't Hooft-Nobbenhuis symmetry. Kaplan and 
Sundrum propose that to each matter field ijj there is a ghost companion ijj, the 



Lagrangian is 

( M 2 

£ = - A ) + £matter(^, D „) ~ £ ma tter (?/>, D (4.11) 

where the form of ^matterCVO an d the form of £ ma tter (VO are identical. We see that the 
name ghost is appropriate since the kinetic term of ijj has a wrong sign. 
Ignoring gravity for a while, there is a symmetry between if), if): 

P (4-12) 

Under this symmetry, H — > —if. since if = H(if)) — H(ijj). Namely PH = —HP or 
{P, H} = 0. For a vacuum |0), P|0) = |0), then 

(0|{P,fT}|0) = 2(0|ff|0) = 0, (4.13) 

thus if |0) is an eigen-state of H, H\0) = 0. 

When gravity is introduced, this symmetry is broken. Let P : — > g^, the Hamil- 
tonian 

S 2 

H = -G lJ > lk —— 2A - 2T 00 + (3) R (4.14) 

ohijdhik 

does not anti-commute with P. If the wave function \1/ has P*B = W then 



(4.15) 



S 2 

UP^ = {-G lhlk + (3) R-2A- 2T 00 )P* 

dhijdhik 

= (-G lJ > lk - J] - +^ i? - 2A + 2T 00 )^, 
dhijdhik 

we deduce Too^Z = 0. 

As an example, consider a scalar <p, then the ghost companion is (f> 

£ = \(d^f - \m 2 4> 2 - A</> 4 + \{d^) 2 - \m 2 ^ 2 - A</> 4 . (4.16) 
Fixing g^, the path integral in a quantum theory is 

[d<f>d4>]e iS ^- iS W. (4.17) 
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The propagator of <p is forward with positive energy, thus ie prescription is used, while the 
propagator of is backward for positive energy, — ie prescription is used. We then infer 
S e ff(4>, 4>) = S e ff(4>) — S e ff{4>). In particular, if there is a term A J \J—g in S e ff(4>), this 
term is canceled by a term in S e ff{4>). 

Now consider the effect of quantum gravity. The interesting aspect of this symmetry 
is that the effect of gravity will introduce a small A. Since there is no ghost companion of 
g^u, the quantum fluctuation of the metric introduces a term A J y/—g. Let the cut-off be 
(U, we must have 

A ~ /I 4 . (4.18) 

Since A ~ (2 x 10~ 3 ef) 4 , |i<2x 10~ 3 ev , or jU -1 > 30 microns. We know that Newto- 
nian gravity is tested above this scale, so it is possible that quantum gravity may break 
Newtonian gravity below this scale. 

We see that the local qauntum effects of ip and ip cancel exactly, while non-local effects 
may not cancel, so it is possible to have a contribution to A: A — > A + /J, 6 /M 2 , this is much 
smaller than [i 4 . 

If tp and V> are coupled, the vacuum is unstable. For example, for g 2 (p 2 4> 2 , there is a 
process |0) — >(/) + (/) + (/) + 4>, the amplitude is divergent: 

= g 2 f n d4 p^ n ^ - m2 )^* 2 - ^ 2 ^ 4 (E^ + **)) = °°- ( 4 - 19 ) 

Let s = (pi +P2) 2 , the integral in the above amplitude is divergent, we need to introduce 
a cut-off for s, s max , and a cut-off for p\ < e, then 

Po^pfc ~ ^ 2 e 2 s max . (4.20) 

Even without direct coupling between a field and a ghost field, their coupling to the 
metric also triggers instability. For example, we estimate 

^ - - 2 X ^J^F^f^ 3 >< ^Gyr)-K (4.21) 

This term is negligible. 

It is interesting to notice that the force between a field and its ghost companion is 
repulsive. Next, we consider the effect of potential of scalar fields. Suppose there are two 
local minima in V(t(j), then there are two local maxima in — V(i(j), if ift runs to the global 
minimum, $ runs to the global maximum, the values cancel exactly. 
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4-4- Symmetry of reversing sign of the metric 



Recai Erdem considered this kind of symmetry ||52|| . Consider action in D spatial 
dimensions 

this is a generalization of the Einstein-Hilbert action. We demand this action be invariant 
under the reflection of the metric gAB — >• ~9ab- Since R — > —R, we have yj (—l) D gR — > 
-y/(-l) D + D+1 gR. Namely 

= (_ 1 )^/2 ) (423) 

or (— l)( r> + 1 )/ 2 = — Ij we deduce D + 1 = 2(2n + 1). When n = 0, the dimensionality of 
space-time is 2, but -D + 1 = 4 is not a solution. 

Since (— 1)(- D + 1 )/ 2 = —1, we know that \J—g changes sign, and the cosmological 
constant term is not invariant, thus forbidden by this symmetry. 

This symmetry is preserved in the action of a scalar field 

-J d D + x x^-g AB d A 4>d B <\>. (4.24) 

The action of a fermionic field is not invariant, unless we demand r ) A — » — r ) A . The worst 
thing is that the stress tensor Tab = Qa^b^ — \gAB^ has a invariant part so the vacuum 
expectation of T 00 is not vanishing. 

4-5. Scaling invariance in D > 4 



There may be many such approaches, a typical one was proposed by Wetterich |k| . He 
postulates a dilatation field £, when we rescale the metric gAB — > X~ 2 gAB, X ~~ >" X^ D ~ 3 ^ 2 ^ 
the action 

s(0 = J ^g{-\eR + f aw) ( 4 -25) 

is invariant. If there is a potential V(£), we require X~ D V(X ( - D ~ 2 ^ 2 ^) = V(£). IfV(£) = 
then = 2IP/(D - 2). When £> = 4, 6, v is an integer. If D > 6, F = 0. 

Let Q be the internal volume, O = J d D ~ 4 xVh 1 let x = the effective four 

dimensional action is 



r = -V^gx z R + V=g(dx)- (4-26) 

There is no term Ax 4 , and x = Xo is a solution. However, there is a problem of stability. 
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5. Anthropic principle 

If the anthropic principle is the reason for A ~ 0, one of the necessary conditions 
is that, A is either a continuous variable, or if it is discrete, the interval AA must be 
sufficiently small. 

Bousso and Polchinski are probably the first to point out that the second possibility 
is realized in string theory J22|(see also |54[] for a review). This is the beginning of studies 



on the string landscape. 

5.1. Bousso-Polchinski scenario 

There are many totally anti-symmetric 3-form fields in string theory, for instance, 
there is a membrane coupled 3-form A^p in M-theory. In the type IIA string theory, the 

(3) 

field C\iJ p is coupled to D2-branes. When IIA (B) theory is compactified on a Calabi-Yau 
manifold CY3, other C fields induces a number of 3-form fields. Let £ p _2 is a p — 2 circle 
in CY3, we have 

f C p+1 -> C^ p , (5.1) 

different p and different £ p _2 gives rise to a different 3-form field in 4 dimensions. 

For instance, in IIA string theory, on £2, on £4; in IIB string theory, 
on £3, but there is no £5 and £1, so and do not introduce 3-form fields. 

Let us focus on a single 3-form field Cp Vp , its four form strength is F4 = 4dC with an 
action 

S= J ( R ~ 2A o) - ^4! F 4 2 ) + S brane , (5.2) 

where Z is a normalization constant. The equation of motion 

a M (v^F^)=0, (5.3) 

has a solution 
this leads to 

PA = PA + f C 2 . (5.5) 

It may be not obvious that the constant C is quantized. Indeed if space-time is really 
four dimensional, C is a continuous parameter. If space-time is higher dimensional, C is 
quantized. For example, in the 11 dimensional M theory, we have 

1 



S = 2nM^ I ^(R - j^Fl). (5.6) 
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A M5-brane is coupled to A 6 , the dual of A3. The coupling is 

2ttM? 1 Ja 6 , (5.7) 

this coupling leads to 



2nM^ / F 7 = 2irn, (5.8) 
J M 7 

this is the standard Dirac quantization condition. Now, F 7 is dual to F 4 

F p...»7 = ±eK-^- v *F Vl ... VA . (5.9) 

Let ni...n 7 E M 7 , then F? = e^ 1 -^ 7 F , and 27rMf 1 V 7 F = 2-ku, thus F is quantized 
F = n/iM^Vr). The reduced action is 

S = 2-nMtM J V^iR - ^If^ 2 ), (5-10) 

thus 

— - = 2ttM 1 9 1 V7 = Z, F = — s = — . (5.11) 

For a M2-bane, the charge is e = 2isM\ 1 , namely the coupling is 27rM 3 1 J A3 = e J A3. 
We have F = ne/Z. Now, for solution F Ml - ^ 4 = Ce^—^/y/^g, we have 



I 77 p 

^o=^ M4 F^=C<=-, (5.12) 



we see that C is quantized. Our conclusion is that, if the action is 



and the charge e J A 3: then C = ne/Z. 
In formula 

Z n 2 e 2 
PA = PA + -^C 2 = p Ao + —, (5.14) 

the dimensions of the constants are [e] = M 3 , [Z\ = M 2 . If e 2 ~ M%, Z ~ MjL 7 , then 
AA > M~ 3 L- 7 . Let M^L' 7 ~ (2 x l(T 3 eV) 4 , then L~ 7 ~ 10 7 GeV, or L -1 ~ lGeV. 
This is possible, but there is a problem: suppose the bare cosmological constant Ao ~ — M 4 
thus n 2 ~ M^Z/e 2 so AA oc ne 2 /Z ~ M 2 e/xf~Z. Use e ~ M 3 , \fZ ~ Mp /2 L 7 / 2 we infer 
AA ~ v /MpL" 7 / 2 ~(2x 10" 12 GeV) 4 , L~ 7 ~ 10" 64 GeV 7 , or L" 1 ~ 10" 9 GeV = leV, this 
is too small, or L is too large. 
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To solve the above problem, Bousso and Polchinski propose to consider multiple 3-form 
fields C^ vX , a = 1, . . . , J. Let Z t = 1 and e — > q i: then 

1 J 

A = Ao + -$>^ 2 . (5.15) 

i=i 

Thus, the additional term is a distance squared in the n-dimensional Euclidean space with 
fundamental lattice spacing qi. The volume of a fundamental cell is 

J 

AV = Y[ gi . (5.16) 

i=l 

Now, consider a shell in between r and r + Ar, the volume of this shell is 

Ojr J_1 Ar. (5.17) 

If this volume is greater than AV, then there is at least one lattice point falling into this 
shell. 

Let A = -r 2 /2, and 

9 h'' 



a = ^(E^-- 2 )= aa - ( 5 - 18 ) 

Let {rii} fall in between r and r + Ar, there must be rAr = AA, Ar = r _1 AA,then 

0jr J-i Ar = o jr ^-2 A A = Oj(2|A |) (J - 2)/2 AA. (5.19) 

This volume must be no smaller than the volume of a fundamental cell, so 

J 

Y[ qi <O J (2|A |) (J - 2)/2 AA. (5.20) 
The physical A is equal to A A, the smallest allowed value is 

wee see that if q t < a/|A |, it is easy to have a very small A. fij, the volume of the unit 
sphere in J-dimensional Euclidean space, is 2n J / 2 /r( J/2). Let qi = 100 _1 a/|A |, we find 
that when J ~ 100, we have roughly 

A 



|Aol 



10" 120 . (5.22) 
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This is how the Bousso-Polchinski scenario solves the cosmological constant problem. 

Just how to generate the right quantum numbers n^? There is the Brown- Teitelboim 



mechanism [ |5"5| , |5"6|| to use. The BT mechanism is similar to electron-positron creation in a 
electric field in 1+1 dimensions. 

In 1+1 dimensions , let there be a pair of charges q and —q, d x E = —qS(x) +q5(x + L), 
the solution is 

(E x < 0; 
E(x) = <E-q 0<x<L; (5.23) 

[e x > l. 

So before the pair-creation, the electric field strength is E everywhere, and becomes E — q 
in between the two charges after creation. It decreases until E reaches the minimum 
-E m i n = E — [E / q]q < q. The rate of pair-creation is 

Pocex P(-^)- ( 5 - 24 ) 

Similarly, a spherical membrane is created in 4 dimensions in a background of F4 and 
F4 — > F4 — C = F4 — q, and the cosmological constant is shifted as 

h^h+ 1 -(F4±q) 2 -\Fl (5.25) 

But this process is very slow and in the end we need the help of the anthropic principle. 
The difference between membrane creation and pair-creation is that the membrane creation 
can either decrease F4 or increase it, however the probability of decreasing A is greater 
than that of increasing A. 

In the original BT model, there is an empty universe problem. This is because in the 
original BT, there is only one kind of membrane, thus q is required to be very small in 
order to have a very small A. But before the membrane creation A — AA is also very small, 
thus the decay rate of this previous universe is extremely small such that inflation in this 
phase makes the universe almost empty. 

This problem is eliminated in the BP model, since before the last transition, the 
previous A is large and the inflation period is short. The major problem remains in the 
BP model is the moduli stablization. 
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5.2. KKLT scenario 



In the Bousso-Polchinski scenario, charges qi as well as normalization constants Zi are 
all moduli dependent. To solve this problem, Giddings, Kachru and Polchinski proposed 
to introduce D7-branes as well as fluxes of and in compactification []57 . 



A D7-brane is specified by a complex function r over CY3 ||58|1 , and form a 
complex 3-form G (3) = F^ - tH^\ Let 

p(s) = F (5) _ c {2) A H (3) _|_ o F (3) = (1 + *) A A . . . A cte 3 ] , (5.26) 



then 

1 

Imr 



dA + ^—dr A iteA = 0, A = etf 6) G (3) - i«G (6) , (5.27) 



where *(6) is the dual operation in the 6 manifold CY3, A is the warp factor in the warped 
metric 

ds 2 = e 2A{y) r]^dx^dx" + e- 2A ^ds 2 CY , (5.28) 

and y are coordinates on CY3. This helps to fix many of the moduli except for all the Kahler 
moduli, for example, the complex scalar p associated with the scale of the Calabi-Yau CY3. 
This is the scalar appearing in a no-scale supergravity theory. 

The warp factor in the GKP model offers us a means to solve the hierarchy problem, 
as in the Randall- Sundrum scenario |59| . 

Consider a complex structure moduli, corresponding to a conifold (z is the size of S 3 
in S 2 x S 3 ), let 

1 r 1 



/ F (3) = 2ttM, — / # (3) = -2itK, (5.29) 

Ja=s* 2na'J B K ' 

then z is stablized to 

2ttK a 2nK 

Z ~ e M 3s , e^min = e 3M Ss . (5.30) 



KKLT first considered a single Kahler moduli |p0 |, and pointed out that there are two 
effects to help to fix the radial Kahler moduli. 

(a) If there is a four dimensional complex sub- manifold in CY3, then wrapping the Eu- 
clidean D3-branes on this sub-manifold forms D-instantons. The contribution to the su- 
perpotential has a form 

W ~ exp(27rip), (5.31) 
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(b) Fluxes induce D7-branes, and the field theory on D7-branes is iV = 1 super Yang-Mills 
theory, and we have 



9ym 

There is a superpotential 



4tt 

— = Imp. (5.32) 



* c 

Let 



W = AeM^)- (5-33) 



W = W + Ae 2niap , K=-2In{-i(p-p)), (5.34) 

then DW = leads to 



W = -Ae- a ^(l + ^aa c ) 1 (5.35) 



where a c = Imp, and we get 



V = -a 2 A 2 e- 2aac /(6a c ) < 0, (5.36) 

thus the cosmological constant is negative and the universe is an anti-de Sitter space AdS^, 
supersymmetry is unbroken. 

If fluxes are not balanced, we need to introduce anti-branes, namely D3 branes. The 
existence of anti-branes breaks supersymmetry, and all moduli will be fixed. Let ao = 
exp(2Ao) be the warp factor, then D3 contribute to the energy density by a factor 

where T3 is the D3-brane tension. Thus, V — > V + AV. Although AV is small, but 
V + AV > 0, the anti-de Sitter space is modified to become a de Sitter space. 

Next, KKLT argue that although the de Sitter space is meta-stable, the life time is 
much longer than the age of our universe, but smaller than the Poincare recurrence time 
t ~ exp(S'o) ~ exp(10 120 ). All the stable and meta-stable anti-de Sitter spaces and de 
Sitter spaces form the string landscape. 

There have been a lot of efforts invested to study statistics of the string landscape. 
There are in principle infinitely many meta-stable vacua of the Bousso-Polchinski type, 
but many of them are not reliable. In the KKLT models, are there upper limits for the 
flux numbers of and H^l Some estimates that there are at least Blf 



10 500 (5.38) 

meta-stable vacua. 

Susskind |62| invented terminology "string landscape", and argued that our universe 
may be multiverse consisting numerous regions in which physics varies. The implication 
of the multiverse is to be discussed in the following subsection. 
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5.3. Populating the landscape and anthropic interpretations 

The vast string landscape itself does not lead to an anthropic interpretation for cos- 
mological constant. Instead, the landscape must be populated 0. In other words, one needs 
a mechanism to produce different universes in a (at best connected patch of) multiverse. 

There are different approaches to populate the landscape. For example, Hawking 
and Hertog |34[ pointed out that the wave function of universe is one populating method. 
Different observers live in their different histories, and they are summed over in the no 
boundary path integral. In other words, observers in different universes live in different 
decohered branches of a single wave function. On the other hand, here we shall mainly 
discuss another better studied scenario: eternal inflation. 

When the quantum fluctuation of of the inflaton 5ip = ¥- is larger that the rolling 
of if in a Hubble time, namely 5ip = ¥- > = fH -1 , eternal inflation occurs |65lJS6| . 
Whether eternal inflation really happens is a matter of controversy [ |67| , |68| , |69| , [70| . For 
instance, the weak gravity conjecture [[FT] may prohibit it to occur |]69|j . 



On the other hand, de Sitter space itself may not be eternal [|f2| , |73[| , it has a finite 
life time, for instance its life time can not be longer than the Poincare recurrence time if 
we view this space-time has finite dimension of the Hilbert space. It is arg ued in 
that a universe with conditions all the same as our universe except the CMB temperature 
is higher is more likely, the probability of its occurrence is ~ e Si ~ s where S = entropy of 
the pure de Sitter space, Si = entropy of a particular universe. It is also pointed out that 
a de Sitter space is a resonant state in the multiverse. 

Anyway, eternal inflation is still a possibility which is semi-classically well defined. 
Before a more complete quantum theory of gravity clarifies all the subtleties, we have to 
take eternal inflation seriously. If eternal inflation indeed happens, our universe is in a 
small part of the eternal inflating universe. The situation is like our earth is a small part 
of our observable universe. This provides a playground for the anthropic principle. 

The validity of the anthropic principle is very controversial, there are a number of 
problems need to be addressed, including: 
• The measure of the multiverse. 

It is intuitive to imagine that a "typical" vacuum in the landscape is kind of vacuum 
that is "realized" in the multiverse most frequently. However, it is rather difficult to 



1 Except that, if one (much more aggressively) assumes every self-consistent mathematical 
structure is automatically "populated" by a higher level of the multiverse. See 
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realize this idea. A number of measures have been proposed in the literature, including 
the volume based measure [[74] (see Jf5| with scale- factor cutoff), the local measure |76|,[77 



the Liouville measure [78,79 ,K0f, the stochastic measure iRlfl, and so on. However, there is 



so far no principle to determine which measure to use. On the other hand, the measure 
problem is also closely related with other problems, as discussed below. 

• How to correctly apply Bayesian statistics? 

The measure problem addresses the a priori probability from the theoretical point of 
view. After that the anthropic probability must be assigned to have anthropic predictions. 
The anthropic probability is proportional to the product of the a priori probability and a 
conditional probability for the theory to be observed, according to the Bayesian statistics: 

, , . Pfselectionltheory x)P(theory x) 

P theory x selection = — y — 1 , - J — ' . . - J — J —. 5.39 

2_^ y P(selection|theory y)P(theory x) 

The "selection" here can be understood as an anthropic effect. Unfortunately, there is 
no principle stating how to properly define the probability from anthropic selection effect. 
For example, Page [|67l , |82| , jS3f assumes typical observers to make observations, Hartle and 
Srednicki J83J argues it is exactly us human observers instead of a typical observer that 
should be used here. On the other hand, some of us [^5j argue a compromise treatment 
seems more natural, that the anthropic probability should be assigned by the probability 
of existence, instead of the number of observers. Alternatively, Bousso, Harnik, Kribs and 
Perez j86| also proposed a more operable entropic measure to mimic the anthropic measure. 

• The problem of Boltzmann brains. 



In the concluding part of [|87[], Boltzmann described the idea of his assistant Schuetz: 
we may come from thermal fluctuations in an extremely large universe, whereas the whole 
universe is in a thermal equilibrium state. Following this idea, one eventually arrives at a 
paradox: the so called Boltzmann problem. 

Now we are not going to talk about whether our galaxies, stars and planets originate 
from thermal fluctuations or not. Instead, we refer ourselves as "human observers", who 
live in a low entropy environment and use entropy increasing to maintain their lives. 
However, if the universe is large enough and in a thermal state (indeed in an asymptotically 
de Sitter universe with Gibbons-Hawking temperature), there will also be other observers 
originated completely from thermal fluctuations, who are by pure chance be in a low 
entropy state themselves, and isolated from the thermal equilibrium environment. These 
observers are called Boltzmann brains, or freak observers. 
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In principle, we can not assert we are not Boltzmann brains but by the next moment 
we will be almost sure, by observing that ourselves are not returning to the thermal 
equilibrium. If we were Boltzmann brains, this probability is exponentially small. However, 
why are we not Boltzmann brains? 

As an example, if our universe starts from a big bang and end up asymptotically de 
Sitter, in a finite comoving volume the number of human observers is finite, because the 
available matter entropy difference will eventually be used up. i However, the number of 
Boltzmann brains is infinite, as long as de Sitter space does not decay. If we were typical 
observers, we should have inferred that we were Boltzmann brains, not human observers. 
This contradicts with the observation in the last paragraph, that we are not Boltzmann 
brains. Similar paradox happens in an eternally inflating universe. 



Even worse, Bousso, Freivogel, Leichenauer and Rosenhaus |38[ argued that time will 
end for eternal inflation with the assumption of typicality. This is too bad since no observer 
is even worse than too many observers. 
• Is there any definite prediction of the anthropic principle? 

It is in debate whether anthropic principle could make predictions. Some may say that 
Weinberg has already predicted the cosmological constant from anthropic principle and 
others may argue that what Weinberg actually did is that he cannot calculate cosmological 
constant from first principle, and simply assumes a probable value from existence of human. 

More generally, in principle one can construct a measure for the whole landscape, and 
predict where is most probable for human beings. If this prediction agrees with current 
experiments, we may take the view point that anthropic principle indeed make predictions. 
But again, the answer depends on the definition of prediction, because the prediction that 
anthropic principle does is logically different from predictions of traditional science. 

To understand anthropic principle, it is good to return to the much better understood 
question why we live on the earth, with an environment surprisingly suitable for human 
beings. However, even in this much easier case, we can not yet make predictions or answer 
whether we live on a most typical planet suitable for intelligence. The way for anthropic 
principle to make predictions for fundamental physics, is thus much longer. 



2 To be more precise, for a human observer, one at least must be able to remember something. 
To prepare the memory to remember anything, there has to be entropy increasing. When entropy 
is maximized, there is no available entropy difference to remember anything. 
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• There is evidence that there are just too many meta-stable vacua (Witten: M2-branes 
— > SU(n) x SU(n)/ (discrete symmetry) — > C^ +n _ 1 x C^ l+n _ 1 vacua. Take m = n ~ 10 100 , 
then the number of vacua of this type is (C^) 2 ~ 2 4n ~ 10 10 °°, a googolplex!) 

Although there are many problems with the anthropic principle, we can not yet rule 
it out. 



6. Tuning mechanisms 

6.1. Brant versus bulk mechanism 



We will mainly explain the work of Kachru, Schulz and Silverstein [q9|. They work 
in brane world embedded in a 5 dimensional space-time. They introduce a tuning scalar 
field, but there is always a singularity in the bulk. As Witten correctly pointed out, such 
a singularity can not be accepted based on general physics principle (Otherwise one can 
introduce just about anything in our world, such as a monopole in the usual Maxwell 
theory). 

Let us start with the action in 5 dimensions 

S=f d 5 xV^G(R - ^(V</>) 2 - he *) + J d A x^{-f{cj>)), (6.1) 

where /(</>) = Ve b< ^. If for any V, one can always find a flat 4 dimensional space-time 
solution, then the tuning mechanism is successful. This is because V is quantum corrected 
(with proof of the form of f(4>) be invariant). If /(</>) is not invariant, one need to show a 
flat solution still exists. 

The equations of motion derived from ( |6.1| ) are 



V^G(Rmn - \g MN R) - ^V^Vm^Vat^- ^GW(V0) 2 )+ 
1 



2 [Ae°^ v—GGmn ~ V^99^ M S" N 5(x 5 )} = 0, (6.2) 

V^G(-A<P - ahe a(t> ) - bV8{x 5 )e h(t> \f^j = 0. 
3 

Let the metric be 

ds 2 = e^^rj^dx^dx" + dx 2 5 , (6.3) 
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then 

-<P" + —A'</>' - aAe a ^ - bV5(x 5 )e b ^ = 0, 

6(A') 2 -^(^) 2 + ^A^ = 0, (6.4) 

3A" + ^ + ^e b *V6(x 5 ) = 0. 

To solve the above equations, let A' = ac/)' ', consider the following cases separately. 
• A = 0. 



When x 5 ^ 0, from 
we deduce a = ±1/3 and 
Thus 



"4 '3 

This solution has a singularity at x 5 = — 3c/4. 
Solution 1, let 

\ x 5 <0; 



thus 



6«V 2 = (6.5) 



b" ± ^<// 2 = 0. (6.6) 



3 4 

= ±-ln|-x 5 + c\ + d. (6.7) 



« = ^ 3 i n (6-8) 

1 -4 x 5 > 0, v ' 



f Jln||x 5 + ci| + di a; 5 < 0; , , 

^ = ^3^. +C2|+rf2 Xs>0 _ (6-9) 



The continuity at x$ requires 



4 1 1 I 
From equations of motion we obtain 

2 r 36 1 



3 3 
ln|ci| + di = — - In j c 2 1 + d 2 . (6.10) 



(6.11) 



The solution for the parameters always exists no matter what value of V is. For an 
arbitrary /(</>), we have 

|(0 / 2 (O)-0' 1 (O) = / / (0(O)), 

3 (6.12) 
3(« 2 0' 2 (O)-a 1 / 1 (O) = --/(0(O)). 
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Or 

-§(cr 1 +c 2 - 1 )=f(7ln|c 1 |+d 1 ), 

(6.13) 



c 2 1 -c/ = -^/(-ln|ci|+di), 

solution to these equations always exists. 

Since the effective Newtonian constant is 

G~ J cfe 5 e 2A( * 5) ~ J dx 5 \^x 5 + a\* . (6.14) 

This contant diverges if there is no cut-off on X5, this is so when c\ < for £5 < 0, thus 
we let ci > 0, then the lower limit of the integral of £5 is x$ = — 3ci/4. This lower limit is 
a singularity, however. Similarly, for £5 > 0, c 2 < 0, there is also a singularity. 
For an arbitrary function f(<f>). we require 

/(0(O))>O, -|<£<|. (6.15) 

• Fluctuations 
Let 



<//«/ = e 2A ?7 MI , + V,, (6.16) 



then 



3 

there is also a singularity, although — > 0, but 



4 1 

hpv oc |-x 5 + c|2, (6.17) 



1 -%c\ d> — > —00 weak coupling: ,„ 10 , 

^5 -> S 3 , t (6.18) 

4C2 <p — )> 00 strong coupling. 



Solution 2, in this case «i = «2 



and b — =f|. Or 



One finds 



<f>M = <tl l *\l X5 t Cl \tt 1 X5<0 ' (6.19) 

' ±|ln||x 5 + c 2 | +d 2 x 5 >0, v ' 



/'(0(O)) = =f|/(0(O)). (6.20) 



ci = -C2 = c, d!=d 2 = d, e T ^ = liL (6-21) 

v C 



A^0. 
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In this case a = — and 

9a ' 



2, r a =F VB „ „ A , 



The junction condition leads to 

i— 4 

V=-12aVB, b=—, (6.23) 

9a 

where V is a function of a and A, and is fined tuned. This is similar to the Randall- 
Sundrum scenario. 

When a = (A ^ 0), let h = <f>' , g = A', we have 

h! + 4hg = 0, 

6<7 2 - \h 2 + ^A = 0, ( 6 .24) 

, 4 n 

3^' + = 0. 

The case A = may be justified in string theory provided there is supersymmetry in 
the 5D bulk. It is expected that we leave /(</>) to be a general function due to quantum 
corrections. 

The major problem of this new tuning mechanism is the existence of singularities. 
The difference of this mechanism from that of Randall- Sundrum is that we need to fine 
tune relationship between the 4D cosmological constant and the 5D cosmological constant 
in the latter mechanism, while here the scalar field <fi massages the correcstions to V to 
the bulk. 

There are also a number of other brane world approaches to dark energy. For example, 
ipTJ considers codimension two branes. Other approaches include 



6.2. Black hole self- adjustment 

This mechanism is due to Csaki, Erlich and Grojean [|93 
Introducing a black hole in the 5D bulk with metric 



y 2 



/ 

where 



ds 2 = -h(r)dt 2 + — dT? k + /i _1 (r)dr 2 , (6.25) 



Mr) = (6.26) 
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and dE 2 , is the metric on the maximally symmetric space with spatial curvature k. Now 
suppose there is a three-brane, we need to require a Z 2 symmetry (so that there is one 
black hole to the left and another black hole to the right), for r < r , the metric is given 
by the functions h(r), r 2 and /i _1 (r), and for r > tq: 



r2 2 ~,2 

rix 7 _1 / \ , _i / ' n \ ' n ~ / ' n 



Mr) = M^), /^(rH^RR, f = (^) 2 , (6.27) 



9 

namely, the metric when r > tq is 

Cfe2 = _^(!l )dt 2 + ( i)2 rfE 2 + h -l( r lA dr \ ( 6 .28) 

If the brane location ro = -R(t) is a function of time, the equation of motion 

p + 3(p + p)^=0 (6.29) 
has a static solution R = R = p = 0. The junction conditions are 

6 v / /i(r ) = K§pro, 18/i'(r ) = -«f (2 + 3«;)/9 2 r , (6.30) 
where w = p/p. We have 

fi = -j-4(l + w ) p 2 r 2 . (6.31) 

If /it < 0, there is a naked singularity in the black hole solution. So we require p > thus 
u; < — 1 and the positive energy condition is violated. Also 

72 1 

p= —A. ( 6 - 32 ) 

this is fine-tuning. In order to avoid these drawbacks, we need to introduce new parameters. 
Assume that this is a 5D Abelian gauge field A M and A r is even under Z 2 and odd 
under Z 2 - The black hole solution with charge Q is 



, r z p Q z 



Q 4 < ^i 2 - 

The junction conditions are 

-2 ., r\2 



(6.33) 



oc( r ^ , Q \ 222 
36 (72 " J2 + 74) = K 5P r , 



'0 
,2 ,, on2 



36(^-4-^-) = -«I(2 + 2-)p 2 4 

6 / n I r 



(6.34) 



' 
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The solution to these equations is 

a* = 2(r 2 + ±-4w P 2 )4, 

6b (6.35) 
Qi = 2(r 2 + -4(l + 3w)p 2 )r 6 .. 

The condition Q 2 > leads to 

^° = v^l^r w< 4 (6 - 36) 

The condition for the existence of a horizon is w > or w < - 1, and p > p- where 

P- = A(^3(l + 6w-3«; 2 + v / (l + «') 3 (l + 9«;))i (6.37) 

When w < —1, p_ < p < po, we find 

/ r ° - 1 1 

dr^[^R - -F 2 - A + £ ma tte r 5(v / ^r(r-r ))]. 



Since 



we have 



Matter ~P~{ ^ 

R = -h"{r) - 6^ - 6^ + (2h' + 12-)8{r - r ), 
n r z r 



(6.38) 



(6.39) 



5 = nfr 2 H h{r H ) = 0, (6.40) 

this tells us that the 4D effective cosmo logical constant is zero A ef f = 0. 
• Other cases 

Take r = R(t) and 

( R e m k = ; 

R(t) = I Swh(Ht)/H k = -1; (6.41) 
( Cosh(Ht)/H k = 1, 

when A4 > and 

i2(t) = cos(Ht)/H, (6.42) 
when A4 < 0. The junction conditions are 

1 2 2 / ^ 7—2 Q \ 

R 2 = 36 K5P ~ { J? + ' ~ R± + R$ h (6.43) 
p + 3H (1 + w)p = 0. 



36 



For w = — 1, p = 0. 
* de Sitter 

We have w = — 1 and 



fi = o, g 2 = o 
K 2 p 2 _ r 2 = ^ > o 



1 2 2 2 2 ( 6 - 44 ) 

2 „2 7 — 2 7lt2 



* Anti-de Sitter 
When w = — 1, 
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# 2 = r 2 - -U 2 p 2 > 0, 

36 5P (6.45) 



p = Q^ = 0. 



When = — i 



36 5P (6.46) 
H 2 = T 2 . 



When = 0, 



p = 0, g 2 = 7 l K 2/Vo- (6-47) 
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For all the above cases a fine-tuning is required. 

Finally, we note that Lorentz symmetry is violated. Since the metric is 



ds 2 = -hdt 2 + r —dYl + h- 1 dr 2 , (6.48) 

t 

the speed of gravitational wave is not equal to c, the speed of light. On the brane, the 
speed of light is c = ^ = ^P-. When transverse to the brane, the speed of light is ^ = h. 
Detailed calculation shows that the speed of gravitational wave depends on the parameter 
E/\p, where E and p are conserved quantities. 



7. Modified Gravity 

Modified gravity is now a huge category, 
classes of models most related to dark energy, 
and briefly mention other directions. 
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Here we shall choose to introduce several 
Namely f(R), MOND, DGP type models 



7.1. f(R) models 

Well before dark energy and inflation have been proposed, there have been already 
attempts to replace the Ricci scalar R in the gravitational action by a general function 

f(R) H- 

It is well known that for the standard Einstein-Hilbert action, one can either treat 
metric itself as dynamical variables, or treat both metric and connection as variables when 
doing variation. The former is known as the metric formulation, and the latter is known 
as the Palatini formulation [p5|]96|l I. 

To generalize to f(R), there are thus two possibilities, the metric generalization or the 
Palatini generalization. Turns out that these two are not identical for f(R). Here we shall 
mainly discuss the metric formulation. The readers interested in the Palatini formulation 



are referred to the review 9£ 



It is also helpful to note that these f(R) models can be related to standard gravity with 
a non-minimally coupled scalar field, using conformal transformation [99 J100I1 . A conformal 
transformation is not a symmetry of general relativity nor its f(R) generalization. However 
there are ways to match observables such that calculation on one side can be used on the 
other side. For these aspects, the readers are referred to the review ||1U1| . 



Now we shall review the f(R) gravity models applied to dark energy @. To be specific, 
we here review two simple models: 
• The CDTT model 

The whole thing started with the paper of Carroll, Duvvuri, Trodden and Turner [ |103| . 
The idea is to use other Lagrangian terms to generate accelerated solutions. Although these 
terms are more complicated and more unnatural than the Einstein's cosmo logical term, we 
can not logically rule out these possibilities and need to test them. The simplest example 
is 

S = Yq^Q j d 4x V~9( R - ~jo) + ^matter, (7.1) 

The equation of motion is 

4 1 4 

(1 + ^)iV - -(1 - jp)gp V R + - V ifl V u) ]R~ 2 = SnGT^. (7.2) 



3 The so called Palatini formulation of general relativity is perhaps also discovered by Einstein 
See a historical review [97] for the story for readers who do not speak German or Italian. 

4 Here we just introduce the four-dimensional f(R) models. For five-dimensional f(R) gravity 



model, see |102 |. 
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Consider matter as an ideal fluid = (p + p)u^u u + pg^ u . Since 



R = 6 



+ 



6 (h + 2H 2 ^j , 



(7.3) 



we obtain the modified Friedmann equation 



3M 2 tf 2 - 



/i 4 M 2 



12(H + 2H 2 ) 3 



(2HH + 15H 2 H + 2H 2 + 6H 4 ) = p. 



(7.4) 



The other Friedmann equation is not independent and can be obtained from the above 
equation and the continuity of equation of p. 
We make a field redefinition 

,4 



D oap 



and 
where 

then 

S(g, V?) 
where 



9^ = eV^P -1 ^ = e^g^, 



Tp V = (p + p)u^u u + pg^, 



p = e - 2aifi p, p = e~ 2a ^ 



1 



16ttG 



d xy / —gR+ / d x\/—g 



P, 



1 



(7.5) 

(7.6) 
(7.7) 



— {dip) 2 -V {ip) 



+ ^matter (7- 



(7.9) 



V(ip) = n 2 M 2 e- 2atfi Ve ai fi - 1. 

Thus, the modified theory is equivalent to a tensor-scalar theory, although matter is not 
canonically coupled to the new metric. The new cosmological equations are 

3H 2 = M~ 2 {p + Pfp ), 
<p + SH(p + V - 1 - 3wV6p = 0. 



(7.10) 



The last term in the equation of motion of <p arises from the coupling between matter and 
<p. Thus, the theory is equivalent to the theory of quintessence with coupling to matter. 
The evolution of matter is 



39 



(7.11) 



There are three separated cases. 

1. Eternal de Sitter, V = 0, but this solution is unstable. 

2. Power-law acceleration. 

When e av > 1, V ~ n 2 M*e~J* M * Xv , then 

a(t) oc t%, a(t) oc t 2 . 



(7.12) 



The problem of this model is obvious, it is similar to the Brans-Dicke theory thus in 
general violate the equivalence principle. 



Chiba [p.04|| considered more generally the model 

S = \m 2 v J d^x^f(R) + S matter (g), 



which is equivalent to 



S 



-M 2 
2 p 



then 



S = \M'p I d'rJ-a 



R-^(d ( p) 2 --L( ip f-f) 



matter • 



As for the Brans-Dicke theory, we require 



A UJ + 1 

7-K2.8X10" 4 , 7 = ^7- 

u + 2 



(7.13) 



(7.14) 



The variation of the above action with respect to tp leads to (p = R if f'((p) ^ 0. 
Redefine 

9iav = f'{<p)9vLv, (7.15) 



(7.16) 



(7.17) 



For the CDTT model, u = 0, so the equivalence principle is violated. For the Starobinsky 
model f(R) = R + M~ 2 R 2 , we need M ~ 10 12 GeV, but the constraint (7- 1| < 2.8 x 10" 4 
is derived for m < 10 -2 GeV. 
• Modified source gravity 

To avoid the above problem, a new model was proposed ||105|| . Again, the idea of 
conformal transformation is used. Let ip = ^ In f'((p), and 



^ = e 2 %„ (7.18) 
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then the action becomes 

1 



d x\/ — g 



(7.19) 



where Xm denotes matter fields collectively, and 



(7.20) 



There is still a problem to be consistent with the solar system, to avoid this problem, some 
introduced Palatini formulaition, or remove the kinetic term of ip. The model without 
the kinetic term of ip is the modified source gravity theory, ip in this model becomes a 
Lagrangian multiplier 



S = I d A x\f—g 



+ S m (e-^~g, X m). 



(7.21) 



When there is no matter, tjj is a just a number, and V(tp) becomes a constant. 
Rescaling back to g^ u 



S 



d x\f^g 



1 



2 -- P 



+ S m (g,Xm)- 



(7.22) 



In the above action, the coupling between g^ v and matter is simpler, and there is a super- 
ficial kinetic term for ip. The Einstein equations 



M^G^ = T^( X ) + T^), 



(7.23) 



where 



T^(V) = -2d^d v ^ + 2V M V„V - <7 M „[e 4 V(V0 + (difj) 2 + 2VA# 
The equation of motion for ip 



(7.24) 



(7.25) 



In the following, we use U(ip) = e^V^). There appears kinetics for tjj, but, taking the 
trace of the Einstein equations we have 



1 



-R 



6 6M2 



(-T + 4t/) + (^) 2 + A^, 



(7.26) 
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this combined with the e.o.m for ifi we find 



U'ty)-4U(il>) = -T, (7.27) 

this is an algebraic equation. We consider several cases in the following i. 

1. The solar system 

In this case T^ v {xm) = thus 

4D($ Q ) - U'{$ ) = 0, (7.28) 

the contribution of -00 to vansishes, so 

M 2 G^ = -e 2 ^U(i/j ) gflu , (7.29) 

this is the Einstein equations with a cosmo logical constant, provided the cosmo logical 
constant is small enough, the solar system is fine. 

2. Interaction with matter 

Integrating out tp results in a complicated and un-renormalizable action, so it is better 
to view the model as a low energy effective model. 

3. Cosmology 

Let Tjf — diag(— p, p,p,p), the Friedmann equation 

3H 2 + - I = M~ 2 e~ 2 ^ [p + U{ip)] - 3V> 2 - GHifj. (7.30) 



Let x = lno, ip = Hi[), 



H 2 = (l + ^r 2 [^e- 2 ^(p + um - i, (7.31) 



■3M"2 Xl a 2! 



we see that the effective Newtonian constant is 87rG e ff = M p 2 e 2 ^(1 + ip ix ) 2 . We have 
the freedom to choose the form of U(ip), amounting to choosing i/)(p). Let 



The thermal dynamical properties of f(R) gravity is different from Einstein gravity, which 



is discussed in [106] 
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so when p — > oo, ijj — > and when p — > 0, if) — > — oo. Since U' — 4U = —T = p thus 



4(ap + p) p 



(7.33) 



C/(V)=ap e 4 ^- 7 ln(l-e^)], V>, 



4 v /J ' ^ 4(ap + p) 

The Friedmann equation is 

2 _ 4apo+4p 2 y/p «po + p-^ln^ _ ^ 
1 ap + 4p J 3M p 2 V"Po + P « 2 ' 



Let p^ = 3M 2 H 2 - p then 



1 dlnppE 

w de = — 1 — — — . (7.35) 

3 a In a 



Using the above equations in fitting the data, one finds 

a = 0.98, h = 0.72. (7.36) 

7.2. MOND and TeVes theories 

The MOND theory (the modified Newtonian dynamics) was proposed by Milgrom as 



a substitute for dark matter ||107|| . It proposes to modify Newtonian dynamics at very 
large scales in order to explain the rotational curves of galaxies. 

Since MOND is aimed to modify Newtonian dynamics, the most appropriate starting 
point is an equation for the Newtonian potential. The typical acceleration involved in 
the rotational curves is of order ao = 10 _8 cm/sec 2 . To make use of this fact, the basic 
assumption of MOND is the eqaution 

p(a/a )a = -V$, (7.37) 

where $ is the Newtonian potential, p(a/ao) is a function with the property p — Y 1 when 
a ^> ao- If p(x) = x for x <C 1, then we have 

a 2 GM . . 

— = -W = — 7.38 
a H 

since a = v 2 /r, the above leads to v = const. 
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Bekenstein and Milgrom ||108 | found that the MOND theory coupled to matter can 



be derived from a AQUadratic Lagrangian (AQUAL), 



where / is given by f'(x) = n(y/x). The equation of motion is 



with 



V-(//( )V$)=4ttGp, (7.40) 

a 



m = \ \ 3 < y J| (7.4i) 



Relativistic MOND 



The MOND action AQUAL can be easily generalized into the relativistic case [p. 08 
Again a scalar field \1/ is introduced in this model . Matter is coupled to e 2 ^g^ u , the 
Lagrangian of ^ and the action of a particle in the background are 

£ * = -tttW^V^cUv!/), 

8tvGL 2 {7A2) 



The relativistic MOND model is also sometimes called relativistic AQUAL, or simply 
AQUAL. In the low velocity limit, the action of a particle is 

v 2 

e*ds ~ (1 + $ - * - —)dt, (7.43) 

thus 

a = -V($ + ^), (7.44) 

combined with the action of we obtain MOND. 

There are two problems with relativistic MOND: 
1. When f = f, causality is violated. The discussion of causality and related modification 
of MOND, can be found in fT0| . 



2. Light is almost decoupled from dark matter deflects light more than relativistic 
MOND does. This unfortunately conflicts with the gravitational lensing experiments. 
Thus the relativistic MOND is ruled out. 
• PCG (The phase coupled gravity) 
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In addition to another field A is introduced with Lagrangian [110jT"l I 



£(* , A) = - l -[A^ A» + V ~ 2 A 2 ^^ ) + V(A% 



(7.45) 



the equation of motion is 



AA - r]- 2 A(d^f - AV'(A) = 0, 
V^^Vtf) = r] 2 e*M5(x). 



(7.46) 



If r) is very small, this model reduces to relativistic MOND, but without violation of 
causality. Still, it contradicts gravitational lensing experiments. 
• Deformed metric 

Since both relativistic MOND and PCG contradict gravitational lensing experiment, 



a deformed metric theory was introduced [112]. Let 



~ 9llv = e-™(A(X)g^ + B(X)d^d u V), 



X = -\(d^) 2 . 



(7.47) 



The second is introduced to deflect light more. To make sure of causality is npt violated, 
the sign of B is chosen such that the deflection of light is not enough. 
It is thus motivated to introduce £7 M to replace and 



9nu = e- 2 %. - 2U^U U sinh(2*), 



(7.48) 



where U„ is time-like. 



• TeVeS 

TeVeS is by far the most complicated theory among the MOND like models [ [113 
Besides the metric, there are new fields ip, a, t/^, with U 2 = —1 and 



9fiv = e ~ 2U,jJJ v sinh2(^, 



and actions 



S 



1 



16ttG 



(7.49) 



S(<p,a) = ~ J d A x[o 2 h^ip^ip„+ l -Gr 2 F(ka)y 



-g, 



(7.50) 



K 
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The dynamics of g^ v remains unchanged, however matter is minimally coupled to g^, for 
example 

-mj^^i^it. (7.51) 

The equations of motion 

= 8nG(f^ + (1 - e-^U a f a(p U u) + r M „) + Q^, 

T\muv = CT 2 (W <P,v + - ■ •)) ( 7 - 52 ) 

Q^ = K[g^U M U [l/>p] + ...}. 

In TeVeS, a is a Lagrangian multiplier. This model can avoid problems of violating causal- 
ity and contradicting gravitational lensing experiments, and reduces to MOND in a certain 
limit. Taking a proper function F, it can also generate the effect of dark energy. 

The major problem of TeVeS is that if g^ v is taken as the physical metric, then the 
tensor, vector and scalar perturbations always have propagation speed greater than speed 
of light. Also, the appearance of dark energy in this model is unnatural. 

7.3. DGP model 

This theory is actually quite simple. It postulates that there are two independent 
gravity theories, one on a 3+1 dimensional brane and another in the 4+1 dimensional bulk 
Til! The action is 



S=^2 I d 4 xdyV^GR + 7^2/ d 4 x^g~R + S m , (7.53) 

where the metric on the brane is g^ u = G Mi ,|brane- Dvali, Gabadadze and Porratri found 
out the 3+1 propagator. To do so, they first studied the 3+1 propagator of a scalar with 
action 

S = M S J d 4 xdyd A (pd A <p + M p 2 J d 4 xdy5(y)d fl <pd fl <p. (7.54) 
The Green's function satisfies 

(M s d A d A + M^(y)d^)G(x,y;0) = d 4 (x)5(y), (7.55) 

Take G(x, y; 0) as the retarded Green's function, the potential generated by the scalar is 

V(t) = J G R [t,x,y;0)dt. (7.56) 



The idea of brane bulk energy exchange as an origin of dark energy is also considered in [115] 
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Let 



G satisfies 



G R (x,y;0) = J -^e^G R (p,y), (7.57) 



(M V - d 2 y ) + MyS(y))G R (p, y) = 8(y). (7.58) 
In the Euclidean space (after Wick rotation), the retarded function is 



we have 



where 



V(r) = ^7-{sin — Ci( — ) + - cos — [tt - 2Si( — )]}, (7.60) 



Ci(z) = 7 + lnz+ / — (cost-1), Si(z)= [ — sin*, (7.61) 
where 7 = 0.577 . . ., the Euler constant and r = M 2 /(2M 3 ). When r <C r 

1/ M = -8^^ + (- 1 + ^ 1 V^l + °M' (7 ' 62 » 

it has the correct 1/r form in 3+1 dimensions. When r ^> vq 

v ^ = -WmiVt + 0{1 ^ < 7 ' 63 > 

it goes like 1/r 2 , the correct potential form in 4+1 dimensions. 

Thus, ro is a crucial scale, the world behaves as a 3+1 dimensional one when r is 
much smaller than this scale and as a 4+1 dimensional one when r is much larger than 
this scale. Let r ~ 10 28 cm. we have M ~ 10 12 (cm)" 1 ~ lOMeV. 

Back to gravity, there is one problem, namely there is one more degree of freedom. 

Let 

Gab = Vab + h AB , (7.64) 

and take the gauge 

d A h AB = \d B K (7.65) 
the equation for the Green's function is 

(M 3 d A d A + M 2 p 5{y)d^)h^ = Mp(y)d^h 55 + [T^ - ^T]%). (7.66) 
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The Fourier transform of h 



One reads the tensor structure from the above formula 

2 2 3 

Although the tensor structure is not what we want, when r <C ro, one still gets the correct 
Newtonian potential and when r ^> tq one gets the Newtonian potential in 4+1 dimensions. 
• Cosmology 

To study cosmology, let us consider the following form of time-dependent metric 

ds 2 = -n 2 (t, y)dt 2 + b 2 (t, y)dy 2 + a 2 (t, y)dT, 2 , (7.69) 
with the Einstein equations 

G A b = k 2 T ab , T A b = diag(-p B ,Ps, • • • + ^diag(-p,p,p,p,0), (7.70) 

o 

Solving the Einstein equations with the junction conditions 

4 (H 2 + \)-^P = <H 2 -^- PB -- A + 4)*. ( 7 - 71 ) 



2k 2K a 2 ' 6' v 2'° a 4 a 2 

where e = ±1. Take ps = C = 0, we find 



Consider two limits of the above equation. 

This is the case in the very early universe when the energy density is large, we have 



r2 , k P 



this is the usual Friedmann equation. 
2. p < ^ 
We have 



» 2 + ^ = i^(l + ^ ("4) 
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Take e = 1, then 



H 2 + 



PA 
3M p 2 



(7.75) 



this is the Friedmann equation with a cosmological constant p\ = 3M 2 /rQ, or A -1 = Tq. 
Thus, the universe is automatically accelerated in the late stage in the DGP model. 

What is the physical meaning of e? Consider n, a and b as functions of y near y = 



n(t,y) = 1 + e\y\a(a 2 + k) ». 

a(t,y) = a + e\y\(a 2 + k)K 
b(t,y) = l. 

Let O ro = l/(4rgif£), O fc = -k/(H$al), we have 



(7.76) 



this is different from the usual Friedmann equation 



+ + — 1- 



(7.77) 



(7.78) 



O ro can be regarded as the fraction of dark energy. 

Later, Nicolis, Rattazzi and Trincherini proposed a Galileon model of gravity as a 
modification of DGP, aiming at curing the ghost instabilities of the DGP self-accelerating 
solution ||116|| . A simplified version of ||116|| , as a scalar field model is reviewed in the 
section of phenomenological models. 



7.4- Other modified gravity theories 

There are perhaps too many modified gravity theories to review here. In this section 
we briefly mention some other modified gravity models. 
• Brans-Dicke and scalar tensor theories 

In the early 1960s, there have already been modified gravity theory known as Brans- 
Dicke gravity [ p.!7| . The gravitational part of the action is 



S 



d Xy 



-g 



(7.79) 



where a; is a constant. The theory is later generalized to f(<p, R), known as the scalar tensor 
theory ||118|| . In addition, this theory has also been extended to the case of 5-dimensional 

EH. 
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• Gauss-Bonnet gravity and Lovelock gravity 

Gauss-Bonnet gravity [ [L2U| is proposed in the mid 1980s, aiming to derive a low 
energy effective gravitational action from string theory In D dimensions, the action of 
Gauss-Bonnet gravity has an Einstein-Hilbert part, plus a correction 

J d D x^~g Q, g = W vp °R^ p(J - AR^R^ + R 2 . (7.80) 

This is the only ghost free combination at the R 2 order. However, in four dimensions the 
Gauss-Bonnet term is a total derivative. In order to have cosmological implications, in 
12 H| , Nojiri, Odintsov and Sasaki coupled the Gauss-Bonnet term g to a scalar field, and 



obtained a model of dark energy. 

Well before Gauss-Bonnet gravity becomes known in string theory, its more general 



form is already there in the early 1970s. This is known as the Lovelock gravity [|122f . The 
Lagrangian of k-th order Lovelock gravity is the summation 



C = (7.81) 



m=0 

where C m is called the Euler density: 



l~m - * c 1 d 1 -c m d m tt a 1 b 1 ' " \,i ra l'-° Z J 

The term m = corresponds to a cosmological constant; m = 1 corresponds to the 
Einstein-Hilbert action and m — 2 corresponds to the Gauss-Bonnet term. 

Unfortunately, in four dimensions the m > 3 terms are simply zero. Thus Lovelock 
gravity cannot be directly applied to dark energy. Nevertheless, 3rd order Lovelock dark 
energy models have been proposed from a dimensional reduction [ |123| | . 
• Horava-Lifshitz gravity 

One of the biggest problems in quantum gravity is that the Einstein gravity is non- 
renormalizable. Horava ||124|| suggested that (at least power counting) renormalizability 



can be obtained by giving up Lorentz invariance at high energies. The idea is to intro- 
duce higher order spatial derivatives while keep the time derivative at second order in 
the equation of motion. For cosmological applications, Saridakis proposed a model of 
Horava-Lifshitz dark energy ||125| . 
• f(T) gravity 
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Alternatively, torsion might also be useful in modified gravity. Bengochea, Ferraro 
and Linder ||126| , |127|| considered the cosmology of a so called f(T) gravity theory. The 
gravitational action is 

M 2 r „ 

S = — v - d A x eT, (7.83) 



2 

where e = dete M a , and e^ a is the tetrad field satisfying g [lv e^ a e v \ ) = r] a b- The scalar T is 
defined as 

T = S/ V T% U , TV = -e p a {d»e a v - d u e%), (7.84) 

= \{K»\ + 5%T°\ - 5%T^ a ), K»% = -\{T»\ - T v \ - T/"). (7.85) 

It is shown that the f(T) theory is also a candidate of dark energy [ |126| , |127|| @. 
• Conformal gravity 

Mannheim [ 129 ] suggests that the gravitational action could be modified using the 
conformal Weyl tensor as 



-a g 



j d 4 x^9C^ pX C^ pX . (7.86) 



It is shown that this model is on the one hand fits well the galactic rotation curves [130 



on galactic scales, on the other hand could behave as a component of dark energy on 



cosmo logical scales ||131[| . 
• Fat graviton 

Sundrum ||132|| proposed a scenario of fat graviton, in which the graviton has a size 
instead of localized. Unlike the other approaches where an action is proposed at the first 
place, the fundamental principle of a fat graviton scenario is yet unknown. Instead, the 
work [|132|1 aims to address how can a graviton be fat, without violating known physics. 
Especially, with the help of non-locality, the fat graviton is able to "know" whether a loop 
diagram is a vacuum loop or a loop attached to external legs. In this way, the equivalence 
principle is preserved. However, it is not clear what detailed rules a fat graviton obeys. 
Also, the size of graviton is required as an input parameter. To solve the cosmological 
constant problem, the graviton size I > 20 microns seems to be a large number in the 
particle physics point of view. 

There are also a number of scalar field models which can be viewed as modified gravity 
models. This is because the scalar part of gravity is coupled to matter. We shall review 
this class of models later, as part of the phenomenological models. 



7 In a latest work [128], a f(R,T) gravity model, where the gravitational Lagrangian is given 
by an arbitrary function of R and of T, is proposed. 
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8. Quantum cosmology 

As reviewed in Section |3|, the approach of wave function of the universe either results 
in a probability P(A) = exp(^), or its exponential. The predicted cosmological constant 
is not only too small, but also conflict with the thermal history of the universe. Several 
ways to solve the problem have been proposed. 

For example, Firouzjahi, Sarangi and Tye [|133| , |134|1 brings the decoherence effect to the 
quantum creation of universe, and found the preferred classical universe is born at string 
scale, which is good for inflationary cosmology. As another approach, Hartle, Hawking and 



Hertog [|135| , |136| argues that from the Bayesian point of view, the quantum creation prob- 



ability should multiply the spatial volume factor to get the observed probability. Again, 
inflation after the born of the universe is preferred. 

These approaches on the one hand bring the Hartle-Hawking wave function back to be 
consistent with observations, but on the other hand make the whole subject more relevant 
to inflation, and less relevant to dark energy. 

Nevertheless, dark energy is very probably an effect from quantum gravity. Indeed 
a lot of attempts on dark energy are related to quantum theories i. Here in this section 
by quantum cosmology we shall restrict our attention to approaches closely related to the 
wave function of the universe. Other quantum approaches are reviewed in other related 
sections. 

8.1. Cosmological constant seesaw 

The seesaw mechanism provides a connection between particles between large mass 
and small mass, which is extensively used in neutrino physics. Non-technically speaking, 
the (type I) seesaw mechanism relies on the fact that a matrix 

M=(° X ) (8.1) 
\x y J 

has eigenvalue 

y ± vV + 4x 2 
± 2 
When j/>i, the smaller eigenvalue is 

2 

X 



For example, it is argued that cosmological constant is an intrinsic scale of quantum gravity 



[137 1, and can be given by the energy-density of quantum black holes [ 138 1 . 
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which is suppressed by the larger matrix element y. Note that A_|_A_ = — x 2 . Thus when 
adjusting y, one eigenvalue goes up, another goes down, like a seesaw. In neutrino physics, 
the analog of matrix M is the neutrino mass matrix. 

The seesaw mechanism is aimed to explain why the neutrino mass is non-zero and 
small. This saturation appears very similar to the cosmological constant problem. Based 



on this observation, Motl and Carroll proposed on their blogs ||139| , |140|1 the idea that the 
smallness of the cosmological constant may be explained by the seesaw mechanism as well. 
The observation is that the present cosmological constant, the Planck scale and the lowest 
possible SUSY scale has a hierarchy 

M A ~ M| USY /M p , (8.3) 

where M\ = (pa) 1//4 is the energy scale of the dark energy. This equation is similar to 
(0)- 

This idea is realized using the wave function of the universe by McGuigan |141jj . It is 



noticed that the cosmological constant term acts exactly like a mass term in the Wheeler 
DeWitt equation. Consider a model with a set of massless free scalar fields ip 1 minimally 
coupled to gravity, the homogeneous and isotropic background action can be written as 



S= - I dta 3 N 



(8.4) 



2 

where iV is the lapse function. In (|8.4j) several coefficients are different from the literature. 



This can be understood as a different convention used in ||141||. Write down the conjugate 



momentum and quantize the equation of motion, one gets the Wheeler-DeWitt equation 

+ iJf£+*?y^ -*)• = <>. < 85 > 

where V = a 3 is the volume of the universe, which can be understood as a time variable. 
The scalar fields (p , on the other hand, acts as spatial coordinates. Note that the term 
proportional to A acts effectively as a mass term in this equation. 

Up to now, everything comes as conventional. Now consider two universes with dif- 
ferent cosmological constant and also with "coupling" : 



(-wim + ¥^ + M ?v^ ~ Al ) $1 + ^»* a = 0> 



$.6) 
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irwim + + M iv^ " Al " A2 ) $2 + v/x ^ 2$1 = °' (8 - 7) 



where the effective mass matrix now is 



M 2 



v^T 2 



Ai vA2 

This form is similar to ( |8.1| ). Now one can apply the seesaw mechanism, such that when 

Ai ~ (lOTeV) 4 , A 2 ~ Mp, (8.9) 
the smaller eigenvalue becomes 

(lOTeV 8 ) 

~ n^ - ' (8 ' 10) 

the correct order of magnitude for the current cosmological constant. However, it remains 
curious how the coupling constant is derived from first principle. In [ |142| , Linde proposed 
that some interactions between universes could follow from averaging effects. It remains 
interesting to see whether these interactions are of the type which we review in this subsec- 
tion. Moreover, it remains not clear in which form will interacting universes classically be. 
It is clear that particles should stay in their mass eigenstates for energy and momentum 
measurement but the analog is not clear for quantum universes. 

8.2. Wave function through the landscape 

In quantum mechanics, there is a well-established mechanism of resonant tunneling, 
which is not only tested by experiments but also widely applied in the industry of electron- 



ics. Tye [ |143| applied the resonant tunneling mechanism to the string landscape, which 
provides a possible solution for the cosmological constant problem. 

As a pure quantum effect, the probability for tunneling is typically exponentially 
suppressed. An intuitive understanding is that a bump in the potential, even higher 
than the energy of the quantum state, can not completely block the wave function. An 
exponentially small part of the wave function is leaked into the classically forbidden regime, 
which connects to another classically allowed regime. 
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Fig. 1: A resonant tunneling in quantum mechanics. The analog is applied to 
cosmology for the cosmological constant problem. 



However, things get changed when there are multiple bumps in the potential. As 
illustrated in Figure Fig. 1, When the incoming wave coming from region I has an energy 
equal to a bound state in region III, the probability for tunneling through the whole 
potential to region V is no longer exponentially suppressed by both of the barriers. Instead, 
the tunneling rate will equal to the larger exponential suppression factor (of II or IV) 
divided by the smaller exponential suppression factor. Especially, when these two factors 
are of the same order, the tunneling rate becomes of order one. 

The reason for the odd thing to happen is as follows: In region V, the wave function 
only has outgoing component. In WKB approximation, the mode is ~ e -f pdq . This 
outgoing component can be linked to two modes in region IV, namely an exponentially 
decaying mode (here decay means becomes smaller from III to IV) ~ e I pdq and an 
exponentially growing mode ~ e-l ' pdq . In conventional tunneling events without resonance, 
only the decaying mode is considered because the other one contributes an exponentially 
small part to the wave function in region III. Similarly, there are two modes in region II, 
only the decaying mode need to be considered for tunneling without resonance. However, 
when the incoming energy satisfies a bound state condition in region III 



the decaying mode in region IV can only link to the growing mode in region II, and the 
growing mode in IV can only link to the decaying mode in II, by requiring the wave 
function and its derivative are continuous at points a, 6, c and d. Thus there must be 
one barrier in which the tunneling amplitude is suppressed and another barrier in which 
the amplitude is enhanced exponentially. This is why resonant tunneling can happen in 
quantum mechanics. 





(8.11) 
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Further calculation shows that, for resonant tunneling to happen, the incoming energy 
should be exponentially close to ( |8.11|) . Otherwise the enhancement disappears. For a 



random incoming energy, the expectation value of tunneling coefficient T(I — > V) can be 
expressed as 

rq^v)~ r(RIII)T(IIl4V) . (8.12) 

K ' T(I -»■ Ill) + T(III -> V) 1 1 

It is also straightforward to generalize the above two step resonant tunneling to n steps. 
Tye ||143|| proposed that the universe may have experienced such resonant tunnelings 



in the string landscape. As reviewed in previous sections, the string landscape is extremely 
complicated and in some sense random. Using ( |8.12|) , one can show that the tunneling 
probability is enhanced. Moreover, the string landscape is many-dimensional. The total 
tunneling rate Y can be calculated by 

T ~ n d r , (8.13) 

where Tq is the tunneling rate for a single tunneling event, d is the dimension of the 
landscape and n is the number of steps for resonant tunneling ||144||. This is in sharp 



contrast with the case discarding resonance, where the tunneling rate is V ~ (d x To) n . 

In this picture, the universe starts with a site in the landscape with large cosmological 
constant and tunnels though a series of sites to the state with the current cosmological 
constant, with reasonably large probability. It is assumed that the universe can not tunnel 
to any of the vacua with negative cosmological constant. The argument is that if the 
whole universe tunnels to a vacuum with negative cosmological constant, the entropy 
would decrease; while if a portion of the universe tunnels to such a vacuum, a singularity 
will develop. 

The resonant tunneling picture gives an elegant understanding for the cosmological 
constant problem in the context of cosmic landscape, as long as one considers the mini- 
superspace model, where a closed universe as a whole could tunnel to other sites with 
resonance. In this case the quantum gravity problem reduces to a quantum mechanical 
problem. However, it is quite difficult to figure out what resonant tunneling looks like in 
quantum field theory. For example, in []145|| , it is proved that if the tunneling process, 



in terms of bubbles, satisfies a static boundary condition, there will be no such resonant 



tunneling in scalar quantum field theory. Later, it is shown ||146|| that for contracting 



spherical bubbles, resonant tunneling is possible. Moreover, in [147|, it is shown that the 
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requirements in | 145|| are too restrictive. A most probable escape path does not necessarily 
satisfy the no-go theorem in [ |145| , and a phase diagram is given to show that some of the 
potentials have resonant tunneling in scalar quantum field theory. For the case of gravity, 
and for the estimation of tunneling probabilities in a landscape, further investigations are 
deserved. 



In [p. 44[ i other rapid tunneling events are also considered, for example the Hawking- 
Moss tunneling without slow roll approximation. One reaches the conclusion that rapid 
tunneling can be common in the complicated string landscape thus the tunneling process 
become analogous to an electron in a random potential. Using renormalization group 
approach, Tye shows that after a number of successful tunneling events, the universe 
will settle down with localized wave function due to Anderson localization 1114811 . Similar 



ideas on Anderson localization of universe in a landscape is also discussed in | 14£ ] and 
15C|| . Besides Anderson localization, decoherence of the quantum wave function may also 



eventually localize the universe to one site in the landscape f|l 5 1 1 . 



9. Holographic principle 

The essence of reality have been one of the central problems in philosophy and physics 
since thousands of years ago. Among all the discussions and debates on this subject, one of 
the most profound allegory is the Plato's Cave. In the recent decades, the inverse problem 
of Plato's Cave, under the name holographic principle ||152| , |153|| , has become one of the 
building blocks for a modern understanding of theoretical physics. 

Holographic principle asserts that the world can be understood as a hologram. In 
other words a theory of gravity is dual to a boundary field theory without dynamical 
gravity in one less dimensions. 

The story of holographic principle tracks back to the investigation of black hole 
physics. Starting from ||154|J155|| , it is realized that classically, a stationary black hole 
in four dimensions is characterized by its mass, angular momentum and charge. A para- 
dox arises here that matter that collapses into a classical black hole appears to lose almost 
all its entropy. This is a contradiction to the second law of thermodynamics. Moreover, 
Hawking proved that [ |156| the event horizon of a black hole never decreases with time. The 
paradox becomes sharper in this sense because it becomes hopeless for the "lost" entropy 
to come out. 
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The paradox was conceptually resolved by Bekenstein ||157| , |158| , |159|| , by the conjecture 
that a black hole has a entropy proportional to its horizon area. The conjecture was soon 



proved by Hawking [16C], with the form 

Sbh = A (9.1) 

Later, 't Hooft and Susskind | 152| , 153 | realized that the black hole entropy can be 
understood as a dimensional reduction, or holographic principle. The entropy is viewed 
as degrees of freedom measured in Planck units, which lives on the surface of the strongly 
gravitating system. 

The right hand side of equation ( |9.1| ) is not only the entropy of a black hole, but also 
an entropy bound for any form of matter localized in a spherical region 

^matter < ~^Q' (*^ - ^) 

For example, consider a spherical symmetric shell of photons which falls towards the center. 



Before the photons forms a black hole, the left hand side of (|9.2|) is always smaller than 
the right hand side. When the photons hit the Schwarzschild radius, the bound may be 
saturated when the photons has maximal possible entropy, in which case the wave length 
of the photon equals to the Schwarzschild radius. 

However, as pointed out by Cohen, Kaplan and Nelson [ 161|| , the thermal entropy 
of an effective field theory can never saturate the entropy bound ( |9.2| ). This is because 
thermal energy E and entropy S of a field theory scales with the length scale of the system 
(infrared cutoff) and the temperature as 

M ~ L 3 T 4 , S ~ L 3 T 3 . (9.3) 

At the point just before the formation of a black hole, the field theory has Schwarzschild 
mass, however the entropy takes the form 

S ~ (M P L) 3 / 2 ~ 5 (9.4) 

which is always smaller than the black hole entropy. Thus, for thermal field theory matter, 
Schwarzschild mass behaves as a tighter bound than the entropy bound. In other words, 
let the ultraviolet cutoff of the system be Auv ? the maximum energy density in the effective 
field theory p ~ Ay V must satisfy 

L 3 A^j V ~ E < LM*. (9.5) 

The equation ( |9.4|) could also be written as follows: the maximal allowed energy density p 
satisfies 

p = 3c 2 M 2 L -2 , (9.6) 



where c is a number introduced in 162 
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9.1. Holographic dark energy 

The reasoning of last subsection can be applied to the vacuum, which leads to a 
holographic model of dark energy. The question is how to choose the infrared cutoff. As 
pointed out by Hsu ||163|| , the simplest choice L = 1/H does not work because it has a 



wrong equation of state. Li ||162|| pointed out that if one take L = 1/Rh, where Rh is the 



future event horizon defined as 



;°° dt f°° da , „. 



the energy density (|9.6|) becomes 



p de = 3c 2 M%R 2 h , (9.8) 

which does behave as dark energy, with the Friedmann equation 

3M 2 tf 2 = p de + p m + . . . , (9.9) 

where p m + . . . denotes matter and other components in the universe. 
To see this, note the index of equation of state w de can be defined as 

p' de + 3(1 + w de )p de = 0, (9.10) 

where prime denotes derivative with respect to In a. Take derivative of ( |9.8| ) and use ( |9.10j ) 
to substitute p' de , one can get 



i 2 i 2ynz , Q11 , 

Wde = -3-3HR- h =~S-^c— (9 - 11} 

This equation is about the nature of holographic dark energy itself, which is independent 
of what form of matter is present in the universe. Here Q d e is the relative energy density 
of holographic dark energy, defined as 

Pde C 2 

Qde = 3m/# 2 = Rjlp- (9 ' 12) 

Q de turns out to be a convenient variable for solving equations of motion for holographic 
dark energy. 

Before going to the equation of motion, now one can already find out qualitative 
behavior of holographic dark energy from ( |9.12|) . When the holographic dark energy is 
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sub-dominant (Orf e <C 1), Wd e — — 1/3 thus 0^ e ~ a~ 2 . When the holographic dark energy 
is dominant (Q^e ~ 1), Wd e — — 1/3 — 2c/3 thus the universe experiences accelerating 
expansion as long as c > 0. 

Now consider the universe dominated by holographic dark energy and pressureless 
matter. Take derivative of ( |9.12| ) with respect to In a, and make use of the Friedmann 
equation (|9.9|), one get an equation for fi^e as 



fide 

The equation can be solved as 



°' de -(i-n de )(i + ^^\. (9.13) 



w , s , n , cHi + ynz) cin(i-yng 8in( c + 2y/o^) 

ln(o/a ) = lnO de H ^— - ^— ^ , (9.14) 



where ao is a constant. The phenomeno logical implication of this solution will be discussed 
in later sections. 



It is also noticed in [162] that during inflation holographic dark energy is diluted. 
To have the correct fraction of dark energy at present time, one requires about 60 e- 
folds of inflation. In other words, holographic dark energy provides an explanation of the 
coincidence problem, as long as inflation only last for about 60 e-folds, before which dark 
energy also dominates. The detailed implication for inflation from holographic dark energy 



is considered in [164 



As an energy component with negative pressure, one might question the stability of 
holographic dark energy against fluctuations [ |165|| . However, the dynamics of holographic 
dark energy is actually not governed by that of a perfect fluid. Instead it is governed by 
the dynamics of the future event horizon. The fluctuation of the future event horizon can 
be written as 

5p de = -Zpde-^ 1 - (9.15) 

This fluctuation can be analyzed using cosmic perturbation theory. Spherical symmetric 
perturbations are calculated in [|166|| , and it is shown that the perturbation approaches a 



constant at late times thus the background is stable against the fluctuations. 

Phenomenologically it is interesting to investigate interactions between holographic 
dark energy and matter components |167[ j . The interaction can be added to the continuity 
equation as 

p' de + 3(1 + W de )pde = 36pi, (9.16) 
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where pi can be set to dark energy density p^ e , matter density p m , critical energy density 
p c , or their combinations. 

Finally, one should note that holographic dark energy solves both the old and new 
cosmo logical problem in a consistent way. This is unlike a number of other models, in which 
one have to assume the solution of the old cosmological constant problem, i.e. assume the 
vacuum energy to be zero, and propose a small dark energy on top of that. 

9.2. Complementary motivations 



Besides motivation from ||161|| , there are also a number of other theoretical motivations 
leading to the form of holographic dark energy, among which some are motivated by 
holography and others from other principles of physics. We shall briefly review some of 
the motivations in this subsection. 

• Casimir energy in de Sitter space 

The Casimir energy of electromagnetic field in static de Sitter space is calculated in 
|168| , |169|| . The Casimir energy can be written as 

-^Casimir = « H> ( 9 - 17 ) 

where the absolute value of u is the energy with respect of the time of the static patch. 
-Ebasimir can be calculated using heat kernel method with ( function regularization. The 
result is 

£ — - 1 ( ln " 2 - ■< - wm) (| - 1 ln (f )) + ^ 

where L is the de Sitter radius, 7 is the Euler constant and r'(— 1/2) ~ —3.48. Here a 
cutoff at stretched horizon is imposed, which has a distance l v away from the classical 
horizon. Note that the dominate term scales as Efcasimir ~ L/l p . Thus the energy density 
scales as pcasimir ~ M^L~ 2 , which is the form of holographic dark energy. 

• Quantum uncertainty of transverse position 

At Planck scale, gravity becomes strongly coupled and the classical spacetime picture 
breaks down. It is suggested by Hogan [ |170| , |171|| that the Planck scale quantum gravita- 



tional effect could be modeled by a random noise. A particle moving a distance l v will 
have a kick of the same order l p in the transverse direction. As the kick in the transverse 
direction is random, the summation of n kicks results in a random walk of distance \fnl p . 
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Thus when a particle moves distance L, there is a uncertainty in the transverse direction 
of order 

AX = ^fZT p . (9.19) 



This relation looks quite like the energy bound (|9.6| ). Indeed one can show that these two 
bounds are related. When one take L to be the scale of the universe, AX as an infrared 
cutoff gives the energy density for holographic dark energy. 
• Entanglement entropy from quantum information theory 

The vacuum entanglement energy is considered in the cosmological context by [[172 



The entanglement entropy of the quantum field theory vacuum with a horizon can be 
generically written as 

S En t = (9.20) 

where (3 is an order one constant and / is the ultraviolet cutoff from quantum gravity. The 
entanglement energy is conjectured to satisfy 

dE Ent = T Ent dSEnti (9-21) 

where T^nt = l/(27i\R/j) is the Gibbons-Hawking temperature. Integrate equation ( |9.21| ), 
one gets 

£ En t = (9.22) 
where iVdof is the number of light fields present in the vacuum. Thus the energy density is 

Pde = 2c 2 MlR~\ c= v^l. (9 . 23) 
H p h ' 2tvIM p k 1 

Here c is in principle calculable in the quantum information theory. 
• Dark energy from entropic force 

Verlinde conjectured ||173|| that gravity may be an entropic force, instead of a funda- 



mental force of nature. [ 174 | investigated the implication of the conjecture for dark energy. 
It is suggested that the entropy change of the future event horizon should be considered 
together with the entropy change of the test holographic screen. Consider a test particle 
with physical radial coordinate R, which is the distance between the particle and the "cen- 
ter" of the universe where the observer is located. The energy of the future event horizon, 
using Verlinde's proposal, can be estimated as 

Eh ~ N h T h ~ Rh/G, (9.24) 
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where Nh ~ Rft/G * s number of degrees of freedom on the horizon, and ~ 1/Rh is 
the Gibbons-Hawking temperature. Following Verlinde's argument (instead of Newtonian 
mechanics), the energy of the horizon induces a force to a test particle of order ~ 
GEhm/ R 2 , which can be integrated to obtain a potential 

V h 1- = -c 2 m/2, (9.25) 

where after the integration one can take the limit R — > Rh, and c is a constant reflecting 
the order one arbitrarily. Using standard argument leading to Newtonian cosmology, this 
potential term for a test particle will show up in the Friedmann equation as a component 
of dark energy pa e = 3c 2 M 2 .R^ 2 . Again it is the form of holographic dark energy. 
• Holographic gas as dark energy 

The nature of a general strongly correlated gravitational system is not well understood. 
In ||175|| it is suggested that the quasi-particle excitations of such a system may be described 



by holographic gas, with modified degeneracy 

u = u Q k a V b M* b - a , (9.26) 

where uq is a dimensionnless constant. Inspired by holography, when taking T oc U -1 / 3 
and S oc V 2 ^, one needs b = (a + 2)/3 and the energy density p can be written as 

where S, T, V are the entropy, temperature and volume of the system. Applying to 
cosmology, S = 8n 2 R 2 M 2 and T = l/(2irR), one obtains 

p = 3^m 2 R~ 2 . (9.28) 
a+4 p K ' 

This has the same form as holographic dark energy with 

(9.29) 



2 a + 3 
"~ a + 4' 



There are some other alternative motivations for holographic dark energy. For exam- 
ple, the relation between holographic dark energy and vacuum decay is discussed in [ p.76 |. 



In addition, there are also many extended versions of holographic dark energy, see \\L77 
and references therein for more details. 
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9.3. Agegraphic dark energy 

In this subsection we review another dark energy model motivated form holographic 



physics, named agegraphic dark energy by Cai [|178|| and later Wei and Cai ||179|| . As 
discussed in ||178| , |179|| , there is a subtlety in the original version of agegrephic dark energy 
model [[1 78 1, where cosmic time is used as the age cutoff. Thus here we mainly review the 
so called new agegraphic dark energy model [ p.79[ . 

The agegraphic dark energy model is based on the Karolyhazy ||180|| uncertainty prin- 
ciple 



St = (3t 2 J 3 t 1/s , (9.30) 

where (3 is an order one constant and t p is Planck time. It is noticed in [|178| and [p. 8 If that 
this relation has close relation with holographic principle and black hole entropy bound. 
Based on the Karolyhazy relation, Maziashvili derived an energy density of the vacuum 
energy 

3r? 2 m 2 

(9.31) 



Pde 



3n 2 m 2 



t 2 



where n is a numerical factor as introduced in [ 178 |. Wei and Cai proposed that, when the 
time in the above formula takes the form of conformal time 



n 



dt 

a 



da 
'a Y H' 



(9.32) 



the energy density ( |9.31|) can be well behaved as a dark energy component. 
The equation of motion for Qde takes the form 



dQ r 



■(1 - o 



de) 



1 



2 V^de 



n 



a 



da a 

Note that the scale factor a appears explicitly. 

The index of equation of state for agegraphic dark energy takes the form 

2 V^de 



(9.33) 



Wde 



•1 + 



(9.34) 



3n a 

thus there can be accelerating solutions. Especially at late times when a — > oo, Wd e — 1- 
Unlike the case of holographic dark energy, here Qde, n and a should satisfy an addi- 
tional consistency relation. For example, in the matter dominated universe the conformal 
time r\ oc a 1 / 2 , thus pa e oc a -1 . Comparing with ( |9.34j ), one obtains Ode = n 2 a 2 /4. 
Similarly, in radiation dominated universe the corresponding relation is Qde — n 2 a 2 . 

With the consistency relation, agegraphic dark energy is a model with a single pa- 
rameter. This is different from holographic dark energy where there are two parameters c 
and fi m . 
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9.4- Ricci dark energy 

Another natural choice for cosmological infrared cutoff is the intrinsic curvature of the 
universe. Based on this, Gao, Chen and Shen ||182|| proposed a model of Ricci dark energy 
i. In the Ricci dark energy model, the energy density is 

Pde = -(H + 2H> + ^=-—R, (9.35) 

where R is the Ricci scalar. 

With the proposed form of energy density, the energy density can be solved from the 
Friedmann equation as 

Pde = -^O m oe" 31na + / e- (4 -* )lna , (9.36) 
z — a 

where fo is a integration constant. The pressure can be solved from energy conservation, 



Pde = ~{3^~3j f ° e (9 - 37) 

The model parameter a and fo are to be determined by data fittings. For some extended 
versions of Ricci dark energy, see [ |184| , |185| , |186|| . 



10. Back-reaction 

A universe without nonlinearity is simple and simply dull. Physicists like to start 
their calculation from the linear case because it is mathematically easy and often solvable. 
However, nobody likes to live in such a universe where two waves always propagate through 
each other without any impact. 

In modern cosmology, one of the most important assumption is the cosmological prin- 
ciple, which states that the universe is homogeneous and isotropic on (observablely) large 
scales. In the standard setup, the cosmological principle is encoded into the FRW met- 
ric and all the subsequent conclusions, especially the Friedmann equation, are under this 
assumption 

9 In |183f| , Nojiri and Odintsov proposed a generalized framework of holographic dark energy, 
which contains Ricci dark energy as a special case. 

10 The Friedmann equation, coupled to the matter equation of motion, becomes already non- 
linear at the background level. While the non-linearity we are discussing here is on the perturba- 
tion level. By definition, linear perturbation will not back-react the background. Back-reaction 
from perturbations is possible to show up only when non-linear fluctuations are considered. 
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The cosmological principle is indeed supported strongly by cosmological experiments 
on CMB and large scale structure (LSS). However, at late times and on small scales (sub- 
Hubble scales), the universe is not homogeneous at all. There are all kinds of structures 
in the universe. On the other hand, on scales larger than the observable universe (super- 
Hubble scales), it is not known whether the universe remains homogeneous and isotropic 
or not. Due to the nonlinearity of gravity, these sub-Hubble or super-Hubble scales may 
back-react on to the scale of the observable universe. 



As early as in 1931, Einstein [17] already mentioned that the matter distribution is 
in reality inhomogeneous and the approximate treatment (cosmological principle) may be 
illusionary, when he was trying to explain why the Hubble's value of the Hubble parameter 



is about ten times too large [ r5|] . Half a century later, Ellis re-examined the effect of 



dumpiness on the average, under the name fitting problem |187| , p.88l - This starts the 



modern story of back-reaction as an effective component of dark energy. 

10.1. Sub-Hubble inhomogeneities 

In this subsection, we consider back-reaction from sub-Hubble scale fluctuations. Typ- 
ically, the aim of the sub-Hubble scale back- reaction theory is that, assuming there is no 
cosmological constant (i.e. the old cosmological constant problem is solved), and provide 
cosmological acceleration from back-reaction from sub-Hubble scale inhomogeneities. Most 
details of this section can be found in IU89II . We assume in this subsection that the universe 



is dominated by pressureless dust. 

To consider small scale fluctuations, it is useful to derive local versions of the Fried- 
mann equation and the continuity equation. To do this, it is helpful to decompose the 
derivative of the velocity field into components, where each component has clear meaning 
in the sense of fluid mechanics: 

Ufj.-u = -flV' +°V + ( 10 - 1 ) 

where is the four velocity of the fluid, and h^ v = — u ll u v is the spatial projection 
of the metric. 

On the right hand side of ( |10.1| ), the scalar part 9 is called volume-expansion scalar, 
defined as 9 = # , which measures the local expansion of the fluid. In the familiar 
homogeneous and isotropic FRW universe, 9 is reduced to 3H. The symmetric part is 
defined by a^ v = u^.^ — 9h^/3 and describes the shear of the fluid. The anti-symmetric 
part is defined by = u^.^ and describes the vorticity. 
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Take one more covariant derivative on equation ( |10.1|) , using the commutative rela- 
tion for covariant derivative to relate the derivatives to curvature, and using the Einstein 
equations to relate the curvature to stress tensor, one can obtain 

6+ \e 2 = -4%Gp-2a 2 + 2u 2 , (10.2) 

This equation is known as the Raychaudhuri equation [ |19U| | , which is widely used in general 
relativity. Another useful equation from the combination of the Einstein equations and 
the decomposition equation ( |10.1| ) is 

i# 2 = 8tvGp- l-R 3 + a 2 - w 2 , (10.3) 

which is the local version of the Friedmann equation, where R3 is the spatial curvature on 
the slice orthogonal to a 2 = a^ v a^ u /2, and u 2 = u^ v u^ u /2. Finally the continuous 
equation for matter component takes the form 

p + 6p = 0. (10.4) 

One can also derive equations for the time evolution of a^u and ui^ from the Einstein 
equations, but we will not need them here. 

Now we are about to average these local analog of the Friedmann equations. The 
spatial average operation is defined as 



. „. . . f d 3 X\/q3(t,x)f(t,x) 

(/>(«) = J y 3k j_^ =M, (10.5) 

where gs is the determinant of the three dimensional metric. 

Using this definition, the local equations (|10.2| ), (|10.3| ) and ( |10.4| ) can be written as 

3- = -4nG(p) + Q, (10.6) 
a 

3H 2 = 87rG(p)-~(R 3 )-^Q, (10.7) 

d t (p)+SH(p) = 0. (10.8) 



The equations (|10.6|) , (|10.7|) and ( |10.8|) are called the Buchert equations [|191|| , where the 



averaged scale factor a is defined as 



1/3 



«W = i „ r^r^ ' H=~. (10.9) 
/ d 3 x^g 3 (t ,x) I a 
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Note that a and H used in this subsection denote averaged variables, not to confuse with 
those in other sections. Q is defined as 

Q=l((9 2 )-{9) 2 )-2{a 2 ). (10.10) 

In the Buchert equations (10.6) and ( p.U.7| ), the back-reaction variable Q is the novel 
term compared from the familiar FRW equations. Q can be thought of emergent in the 
sense of coarse graining. When Q > 0, Q will behave as a effective component in the 
universe, which drives the late time acceleration. 

So far we have not said anything about how to choose spatial slices to do the aver- 
age. In fact choosing spatial slices properly is extremely important for the back-reaction 
calculation to make correct predictions, both theoretically and phenomenologically. 

Theoretically, in general relativity there is no preferred choice for the spatial slices. In 
linear cosmic perturbation theory, there is an elegant gauge invariant way to do calculation. 
However, here the back-reaction one considers is beyond linear order and a gauge invariant 
formalism is not available on sub-Hubble scales. 

Phenomenologically, cosmological experiments are typically carried out by measuring 
redshift and distance. When the universe is perturbed one has to make sure whether the 
calculated scale factor is indeed the one which is used to calculate these quantities. 

For example, as pointed out by Ishibashi and Wald |p.92[ , if one averages two discon- 
nected decelerating universes, one could get the conclusion that the "whole" universe is 
accelerating. This absurd conclusion shows that great care is needed to select the averaging 
hypersurface. For this purpose, one is led to consider the light propagation in a perturbed 
space. 



In ||193|| , Rasanen showed that the redshift in a dusty universe can be calculated as 

1 + z = exp j jT° drj Qfl + a^e"^ j , (10.11) 

where the integral is along the null geodesic, and e M is along the spatial direction of 
the geodesic. One can argue that when there is no preferred directions in the sky, the 
integration of cr All/ e M e l/ should be suppressed by averaging effect, as long as one choose the 
spatial slice with statistical homogeneity and isotropy. 

Similarly, it is shown in ||194j| that for angular diameter distance and the luminosity 
distance, similar average effects occur for statistical homogeneous and isotropic slicing. 
These results show evidence that the statistical homogeneous and isotropic slicing is the 
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one that should be used when the calculation of back-reaction need to be compared with 
observations. 

Having determined the choice of slicing, one can focus on the calculation of the back- 
reaction variable Q. Unfortunately the calculation is very difficult. The Newtonian cal- 
culation might be oversimplified, as in Newtonian gravity the two terms in Q cancels up 
to a surface term ||195|| . For calculations in general relativity, only toy models are doable 



currently. Some calculations report there is strong cancellation between the two terms in 



Q, although they do not exactly cancel [196]. Others show that it is not the case [197 



Thus whether or not small scale back-reaction could become important deserves further 
investigation. 



Besides the approach of averaging, there is also a Lemaitre-Tolman-Bondi ||198| , |199| , |200 | 
approach of sub-Hubble back-reaction. The idea is that there are voids in the large scale 
structure. We may live in a void in which the observables such as the luminosity distance 



need to be calculated with more care. For examples, in ||201| , |202| , |203| , |204|1 , it is shown that 
the LTB type back-reaction could mimick dark energy. 

10.2. Super-Hubble inhomogeneities 

Compared with the sub-Hubble theories, the theories of super-Hubble back-reaction 
has a wider variety of goals. Some of the works aim to provide a screening mechanism 
for the cosmological constant, and completely solve the cosmo logical constant problem. 
While some other works, like sub-Hubble theories, aims to give acceleration assuming the 
cosmological constant is zero. The difference is that, screening mechanism gives a negative 
energy density contribution from backreaction, while the latter gives a positive energy 
density from backreaction. 



For the screening mechanism, Mukhanov, Abramo and Brandenberger ||205| , |206|| have 
set up a gauge invariant formalism for perturbations up to second order. The quantities 
they use are gauge invariant after spatial integration (averaging). For scalar field mat- 
ter tp = (po + 5tp, the second order perturbated energy momentum tensor is calculated. 
Interestingly, the second order energy momentum tensor, in the slow roll approximation, 
has an equation of state p = —p > 0, in other words, it is a cancellation or screening of 
cosmological constant. There are also similar results from gravitational loop calculations 
by earlier studies by Tsamis and Woodard [ |207| , |208| . 



However, Unruh ||209|| (see also ||210|| ) pointed out that gauge invariance of spatial 



averaged variables does not guarantee that the calculated effect is accessible or observable 
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by a local observer. This objection is supported by Geshnizjani and Brandenberger ||211 



from an explicit calculation of a scalar field coupled to gravity. It is shown that without 
isocurvature fluctuations, the expansion 9, which is locally accessible, do not receive any 
backreactions from super-Hubble perturbations. The expansion 9 including back-reactions, 
as a function of ip = ipo + <V, nas the same function dependence as the background 
expansion 9q, as a function of tpo: 



9 = y/3GV(4>). (10.12) 

Thus for super-Hubble perturbations, a local observer will observe if and can not find a 
difference in dynamics compared with the background. 

However, with isocurvature perturbations, the situation changes. In [ 212|| , Geshnizjani 



and Brandenberger considered back reaction in two scalar field model and found that in 
this case the back-reaction does not vanish. The reason is that there are different choices 
for local clocks. For different choice of clocks (proper time clock and energy density for a 
scalar field, as considered in ||212|| ), the results of back- reaction are different (actually the 



correction changes sign). Thus it is crucial to identify which is the observable related to 
observations of dark energy. 

Abramo and Woodard ||213| , |214|| proposed a local operator (before spatial average) and 
showed that the back-reaction does not vanish. Further calculation | ]2 1 5| , |2 1 6| , |2 1 7| , |2 1 8| , |2 1 9| | 
showed from one-loop and two-loop calculation that the probe scalar field obtains an 
equation of state p < —p, which may drive a period of super-acceleration. Eventually the 
super-acceleration will be turned off by the non-zero renormalized mass. Based on these 
consideration, a scenario of non-local cosmology is proposed to model the above behavior 
220| , [221| , |222|| (see also ||223|| and references therein). 



The back-reaction from scalar field to de Sitter space as a screening mechanism is 
also considered in ||224| , |225| , |226|| . The authors argue that the de Sitter space analytically 



continued from a sphere should not be used to describe realistic cosmology because it 
corresponds to de Sitter space artificially kept at fixed Gibbons-Hawking temperature, 
which corresponds to de Sitter space sourounded by reflecting walls. Instead, They consider 
a period of de Sitter expansion 

ds 2 = _ dt 2 + a( t)2 rfx 2 5 a{f) = e Ttanh(t/T) 5 (10 .1 3 ) 

where T 3> 1/H is a time cutoff such that interactions from de Sitter background are 
removed in the asymptotic past and future. Correspondingly, the vacuum choice is that 
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the field approaches local Minkowski vacuum at t — > — oo, which is different from the 
Bunch-Davies vacuum. 



The calculations in ||226|| show result as follows. For free scalar fields, as a result of 
de Sitter symmetry, the correction to the one point function is a trivial redefinition of 
the cosmological constant. For interacting fields in the FRW patch of the de Sitter space, 
two point function of the scalar correlator produces non-trivial infrared divergence but 
the one point function is not affected. In global de Sitter space, where the accelerating 
expansion of the universe is initialized by a period of contraction, there are non-trivial 
effects for the one point function for scalar fields. However, for the case of dark energy, our 
universe is clearly not contracting right before dark energy domination. Thus it remains 
an open question whether this approach could provide a dynamical explanation for the 
dark energy problem. Nevertheless, the interaction which breaks the de Sitter symmetry 
is an interesting candidate for dynamical dark energy. 

Some other screening mechanisms are also reviewed in the previous section on tuning 
mechanisms, where the issue of back-reaction is not as relevant as discussed here. 

As mentioned above, there is another class of super-Hubble back-reaction models, 
which produce instead of screen a cosmological constant. For example, Kolb, Matar- 
rese, Notari and Riotto considered the effect of super-Hubble perturbations from inflation. 
These perturbations, viewed from a local observer, may look like a source of acceleration. 
There are also counter arguments on these class of mechanisms, see for example [|227 . 



One should also note that there is a whole literature for loop calculation during in- 
flation, which is another accelerating epoch of the universe. The techniques developed in 
those loop calculations are also helpful in understanding the dynamics of dark energy. But 



we shall not introduce these works here. Interested readers are referred to ||228|| and related 
works. 



11. Phenomenological models 

In this section we introduce phenomenological models for dark energy, focusing on 
modification for known matter components. 

11 Sometimes it is difficult to classify whether a model belongs to modifying matter or modifying 
gravity because they are coupled. Indeed some models reviewed in this section can be thought as 
modification of gravity in the infrared. But we include them here anyway because in these models 
gravity is not modified from the action. 
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11.1. Quintessence, phantom and quintom 

The most well-studied parts of phenomeriological models are models with rolling fields. 
This is because such models are direct generalizations of the cosmological constant: When 
the fields are not rolling, their potential energy behaves as a cosmological constant. Here 
we review these models in a logic order and in a brief way. The readers can find more 
details in the reviews [f7|, |229|| , and of course as well as the original papers. 
• Quintessence 



A quintessence field ||230| , |231|1 is a scalar field with standard kinetic term, minimally 
coupled to gravity. The scalar field part action takes the form 



-l g ^d^d u p-V{p) 



;ii.i) 



where the metric convention is (— , +, +, +) such that the scalar field has standard kinetic 
term. Take variation of g^ v , one obtains the stress tensor. 



T^v = d^pd u p - 



11.2) 



The energy density and pressure can be read off from the energy momentum tensor as 

(11.3) 

As usual, the Friedmann equations are 



p=-<f? + V(<p), p= X -p 2 -V{p). 



3M 2 H 2 = p, 



-2M 2 H = p + p, 



(11.4) 



where, as a reminder, M 2 = 1/ (87rG) . In terms of ip, the above equation reads 



3M 2 H 2 



1 



-p 2 + V, 



Another useful equation comes out of a combination of these equations: 



(11.5) 



a 



a 



6M-M- =p + 3p = 2(p z -V). 



(11.6) 



From local energy conservation, the continuous equations is 



p + 3H(p + p) =0, (p + 3Hp + V'{p) = 0, 



;n.7) 
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where the latter is a rewritten of the former in terms of the field tp. 
The equation of state takes the form 



V 

w = — 
P 



V(<p) 



(1U 



Note that the kinetic term has positive pressure and the potential term has negative 
pressure. When the field rolls slowly, the potential dominates thus w approaches —1 from 
above. 

What kind of potentials can give rise to acceleration? A simple answer, is flat poten- 
tials. To make the answer more precise, one can have a look at the critical case: cosmic 
expansion without acceleration or deceleration. For this purpose, consider the power-law 
expansion 

aoct p , (11.9) 

where we have kept p general, keeping in mind that p = 1 corresponds to zero acceleration 
and p > 1 corresponds to acceleration. The potential driving this kind of acceleration can 



be solved from (|11.4[) as 



Thus potentials flatter than 



V = Vq exp I — 



V = Vq exp 



2 (p 

pW p 



(11.10) 



;ii.ii) 



have acceleration solution. Moreover, as the kinetic term has much larger pressure than 
the potential term, the potential domination epoch is an attractor solution as long as the 
potential is flat. 

Quintessence may be realized using axions ||232|| , dilatons [|233| , in QCD ||234j| , in Higgs 
potential ||235|| , or in unparticle theories ||236| . 



Phantom 

A menace from phantom ||237|| can be expressed in terms of the action @ 



S 



-g 



l 



g^d^dvip - V(<p) 



(11.12) 



The action has a "wrong" sign kinetic term: £km oc ~ l f 2 - Here phantom is also often 
referred to as ghost in the literature. 



12 Alternatively, phantom dark energy can be also obtained in the scalar-tensor models [238]. 
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When expressed in terms of p and p, the equations of motion for phantom are identical 
as written in the quintessence case, while now p and p are expressed in terms of (f as 



The equation of state is 



-\<f + V(<p), P=~<?-V(v). (11.13) 



P ±ip 2 - V{<p) 



Now there are two possibilities for the equation of state: w > 1 for the kinetic dominated 
regime and w < — 1 for the potential dominated regime. The latter behaves as a component 
dark energy with super-acceleration. In other words, the universe will have an acceleration 
faster than exponential. The energy density keeps growing until it reaches infinity in finite 
proper time. When the energy density reaches infinity, the expansion rate of the universe 
diverges and every thing is tore off. This is called the "big rip" singularity, a physical 
singularity where all known physical laws break down. To see this, consider for simplicity 
the constant w case. In this case the scale factor can be written as a simple power of time 
as 

a = ao(t - t )3(^. (11.15) 

When w < —1, the power of t— to becomes negative and one concludes a — > oo when t — > to, 
a finite proper time for a comoving observer. The Hubble parameter H oc l/(t — to) also 
blows up in the future. The reason is that when w < —1, to is in the future instead of in 
the past. Similarly the curvature blows up thus the big rip t = to is a physical singularity. 
The big rip is a disaster not only for civilizations but also for physical laws, which need to 
be avoided or resolved. 

Moreover, there is a quantum instability in the phantom models 111. Once the phantom 
quanta interacts with other fields, even through gravity, there will be an instability of the 
vacuum because energy is no longer bounded from below. The ghost busters Cline, Jeon 
and Moore [ |24U| pointed out that for the phantom to be consistent with CMB observations, 
the Lorentz symmetry must be broken for phantom at an energy scale lower than 3MeV. 
Otherwise the effect of vacuum decay into phantom quanta and photons via gravitation 
could have been observed on the CMB. 



13 However, the statement w < — 1 itself does not necessarily mean an instability. For example, 
holographic dark energy has w < — 1 solutions without an instability. Also in some modified 
gravity models such a well-behaved w < — 1 solution may be obtained [ [239fl . 
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• Quintom 

The quintessence field always have w > — 1 while the phantom field (as dark energy) 
always have w < — 1. Is it possible for a field theory model to cross the phantom divide 
w = — 1? The answer is positive, which is known as the quintom dark energy, a combination 
of the quintessence and the phantom ||241|| . 



Before reviewing what can be done, it is helpful to first have a look at what is under 
no-go theorems. As is noticed from the beginning ||241|| , it is not possible to have a single 
scalar field to cross the phantom divide 0. To see this, consider a model with time varying 
kinetic term ~ f(t)(p 2 . If one wants to have the field cross the phantom divide when / 
cross 0, at / = the field will have vanishing kinetic term thus divergent sound speed. 
More general proofs can be found in |244| , |245| , |246| , |247|| . 

Thus one is forced to consider models with at least two fields, with the matter action 



S 



\d^d^ +ld»<Td li <T-V(< P ,<T) 



(11.16) 



The two-component dark energy has an equation of state 

V 



w 



¥ 2 



;n.i7) 



The model is extensively studied in ||248| , [229|| , and the references therein. Depending 
on the potential, the quintom equation of state may across —1 from above, or across —1 
from below liH. For an extended version of quintom scenario, see [250 . 
• Fast oscillating fields 

Here we consider a scalar field with standard kinetic term and potential V cc ip n , 
focusing on the case that the field is oscillating around its minima and the oscillation rate 
is much more quickly than Hubble parameter. Turner calculated the averaged equation of 
state of this case |[251|| , with 

, > n — 2 . . 

( w ) = (11.18) 



n + 2 



When n C 1, the fast oscillating field can drive cosmic acceleration [[252|J253|| . To see this, 
one can apply the virial theorem. Define G = ipip. We have 



G = 2T - (3Hip + V')ip nV, 



(11.19) 



14 For counter examples on classical instabilities, see [242,243]. 

15 Alternatively, in f(R) gravity models, w may also cross —1 1 249 1 . 
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where T = <fi 2 /2 and we have neglected Hubble expansion in the last expression because 
the field is oscillating fast. Note that G is a total derivative thus the time average vanishes 
when we take the averaging time to be a multiple of the oscillation period. Thus 

(T) = |<y>. (11.20) 

Translating to the averaged equation of state, w takes the form of ( |11.18| ). 

Unfortunately, the fast oscillating field has an instability, as discussed in [|254 |. The 



inhomogeneity due to the instability may get this mechanism observationally disfavored. 
• Dark energy interactions 

It is possible that dark energy is dark but not lonely. Especially dark energy may 
decay into dark matter or vice versa. In ||255| , [256| , |167|| , the following class of interaction is 



considered: 

p m + 3Hp m = 5, (11.21) 

p de + 3H(1 + w de ) Pde = -5, (11.22) 

where 5 could take the form 5 = —bp de , 5 = —bp m , or 5 = —b(p m + p de ), etc. The effective 
equations of state for dark energy and matter are what one actually measures and different 
from their actual equations of state: 

<f = -^, wf = w de + ^. (11.23) 

In some sense, the interaction could provide an explanation for the coincidence problem 
because if dark energy decays to dark matter in the future, the fraction of dark energy may 
remain at the value one observes now. However, why dark energy starts to dominate today 
but not much earlier, as the other half of the coincidence problem, remains unsolved. 

Alternatively, the phenomenological time variation of the cosmological constant can 
also be implemented by a A(t)CDM approach, as described in ||257| , |258| . 



11.2. K-essense, custuton, braiding and ghost condensation 

Here we briefly review dark energy models from fields with modified kinetic terms. 
• K-essense 

Chiba, Okabe and Yamaguchi [ |259|| (see also Armendariz-Picon, Mukhanov and Stein- 
hardt | p60| , |261|| ) proposed a more general framework on field theoretic dark energy, named 



k-essense. The idea is to generalize the kinetic term, as long as the perturbations still have 
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second order derivative and Lorentz symmetry is preserved 1 262 , 263 1. The corresponding 
action in this case is 

S = J d 4 Xy/=g~ p(X, ip), X = -g^d^dvip. (11.24) 

Here X is the conventional kinetic term. Note that p appears in the action is exactly the 
pressure, and the energy density takes the form 

p = 2X Px -p, (11.25) 

where px = dp/dX. In terms of p and p, again one have the the Friedmann equations 

(EH)- 

It is shown [|263|| that the perturbations obey second order differential equation of 
motion as usual (except a special case considered later). To be ghost free, the theory 
should have 

px > 0. (11.26) 
To be perturbatively stable, it is also required that the sound speed is real. The condition 



is 



cl = > 0. (11.27) 

px + ZXpxx 

Further, one might also require the sound speed c s < 1. This is satisfied when 

Pxx > 0. (11.28) 

There is a debate on whether c s < 1 is necessary or not | p64|| . Thus the condition that 
k-essense is well behaved is (|11.26|) plus either (|11.27|) or (|11.28|) . 



To make predictions, more concrete forms of p(X, <p) are needed. The simplest case is 
perhaps power law k-essense, where 

^ & = Q n (1 ~?2 + x ") ■ (11-29) 
9(1 + wyip z 

2 

It is shown that the parameter w is indeed w = p/p, with the late time behavior a oc t 3( - 1+w ) . 
Another example is DBI-essense, considered by Martin and Yamaguchi | p65| . Also, k- 



essence may provide a unified framework of inflation and dark energy ||266|| . 
• Cuscuton 
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The Cuscuton dark energy is introduced by Afshordi, Chung and Geshnizjani [|267| . 
as a singular limit of k-essense. In the cuscuton model, the matter action takes the from 



S 



d xy/—g 



fi\\g^d„<pdM-V(<p) 



;il.30) 



The cuscuton model has an infinite propagating speed for linear perturbations. How- 
ever, the phase space volume of linear perturbations is vanishing, thus no information is 
propagating faster than speed of light. 

Cuscuton is inspired by the plant Cuscuta, because cuscuton is parasitic in the sense 
that the cuscuton itself does not have its own dynamics. The equation of motion takes the 
form 

(3/i 2 #) sgn(<^) + V'{<p) = 0. (11.31) 

The evolution of cuscuton follows from other energy components, which can be derived from 
combining (|11.31 ) and the Friedmann equation. When adjusting V(<p), cosmic acceleration 
can be obtained. For example, when V is an exponential potential, the expansion history 
of cuscuton behaves like DGP (but perturbation theory behaves different). Cuscuton may 
be a minimal dynamical dark energy model because it has no dynamics, while remains 
the dynamical feature for dark energy. The cosmic evolution in the cuscuton model is 



investigated in ||268|| . 

• Kinetic gravity braiding 

K-essense is not the most general form for a scalar field with second order equation of 
motion. As shown in [|242| , [269l (see also [|116| , |270|1 in terms of a simplified yet generalized 
Galileon model), a more general form can be written as 



S = J d 4 x [K{<p, X) - G(cp, X)d 2 (p] 



[11.32) 



The cosmological implications are studied in ||271|J272j| . It is shown that the field could 
behave as dark energy. This class of actions can be further generalized to [|273[ , where 
cosmological implications are so far not studied. 
• Ghost condensation 

As another generalization of k-essense, a model of ghost condensation is aimed to cure 
the quantum instability of phantom dark energy, and give an equation of state p = —p 
to drive the late time acceleration of the universe. The idea is that the instabilities come 
form the perturbation level, where the gradient energy plays an important role. On the 
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other hand, the spatial gradient energy does not show up in the homogeneous and isotropic 
background lever. Thus one can propose a ghost like background action, and let the spatial 
gradient terms cure the instability problem. 

To realize this idea, Arkani-Hamed, Cheng, Luty and Mukohyama |[274|| (see also [ 275|Q 
proposed an effective field theory of rolling ghost, preserving cosmo logical symmetries. The 
action takes the form 



S 



-g [p(X) + (VV) 2 + ■■■]■ 



(11.33) 



The cosmological solution of the model is either (p = or px = 0. In interesting models 
analogous to tachyon condensation with a spontaneous symmetry breaking, the ip = 
solution is unstable and the universe is driven to the px = solution dynamically. This 
results in X = constant, with p = —p (recall ( |11.25| ) with px = 0). Thus the condensation 
of ghost behaves like a component of dark energy. The intuitive understanding of a constant 
equation of state is that, there is a shift symmetry of the system, thus rolling of (p with a 
constant speed results in a constant w. The model of ghost condensation can be thought 
of as an IR modification of gravity (for other IR modifications, see 



Inspired by ghost condensation, an inflation model named effective field inflation ||277 



is proposed, with the most general action preserving cosmological symmetries. In addition, 
A ghost dark energy has also been proposed |p78| and has been widely studied ||279j . 



11.3. Higher spin fields 

In cosmology, going to higher spin mostly means spin higher than zero. There are 
attempts to generalize the scalar field dark energy model to spinors ||280| , |281| . |282j , |283|| , 



vectors ||284| , [285| , and p-form fields ||286|J287jl . Here we shall briefly mention the approaches 
of spinors. 

In curved spacetime, the action of a spinor field can be written as 



S 



d x e 



;n.34) 



The r matrices 



where e = dete^ a , and e^ a is the tetrad field satisfying g^ v e^ a e u j, = r} a b. 
are defined as = e^ a ^ a . The covariant derivative is defined as D^if; = (d^ + ^u) fia f ) T, ab ), 
with u^ab = e v aD^e^b, and T, ab — |[7 a ,7 & ]. If the vacuum expectation value transform as 
a scalar instead of a pseudo-scalar, the ^7 5 ^ term can be dropped. 
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When the field is homogeneous, the energy density and pressure can be written as 



p = V, p = V'ipip - V. (11.35) 

It can be shown from the equation of motion that iftifj oc a -3 . In other words the spinor 
field rolls very fast. Thus to have acceleration (and not coming simply from a cosmo logical 
constant), one need a very flat potential. For example, a potential V oc (tfji(j) n with n < 1/2 
will have acceleration. Moreover, as a spinor has multiple degrees of freedom built in, the 
equation of state w may cross —1 without divergence of sound speed [[281 1. 



For vector fields, one have to take care of constraints and large scale anisotropy. Also 
one has to break the conformal invariance otherwise the energy density decays too quickly. 
We shall not discuss these issues in details here. 

11.4- Chaplygin gas and viscous fluid 

The matter components in cosmology are usually and conveniently written in terms 
of fluids. Most dark energy models have fluid description. The fluid models we refer to 
here correspond to the models that originate from direct modification of fluid properties. 
• Chaplygin gas 



The best studied fluid model for dark energy is the Chaplygin gas model ||213| , [288 
The original Chaplygin gas model has the equation of state 

p = . (11.36) 

P 

This kind of fluid is first used in fluid mechanics which describes the air flow near the wing 
of a aircraft. 

In the presence of matter, the energy density and pressure of the Chaplygin gas can 
be written as 

A 



p-y/A + Ba,- 6 , p = - == =. (11.37) 

VA + Ba~ 6 

where B is an integration constant. 

Although the Chaplygin gas itself is not motivated from field theory models, it is 
interesting that the model can be mimicked by a scalar field with a simple potential, in 
the matter or dark energy dominated universe: 

V(<p) = ^ (cosh 3^+ — • (11-38) 
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The Chaplygin gas model is later generalized to 



in ||289|| . Other kinds of generalization are also possible, for example the form |[2Q0 j 



p = -p[l-swc{p /p)], (11.40) 

where sinc(x) = sin(x)/x, and po is a parameter with dimension [mass] 4 . The possible 
theoretical embedding and phenomenology are also considered there. 
• Viscous fluid 

The local energy conservation equation ( |1 1 . 7|) is by far satisfied by all the phenomeno- 
logical models discussed above. However, it is not always the case because energy can leak 
to other forms such as the thermal motion from bulk viscosity. Brevik and Gorbunova 
noticed that the modification due to bulk viscosity leads to a model of viscous cosmology 



291| , |292|1 . The energy momentum tensor of the fluid takes the form 

Tfu, = pu,j,u u + (p-3iIC)V" h^ = g^v + u^u,,. (H-41) 
The local energy conservation equation now reads 

p + 3H{p + p) =9H 2 (, (11.42) 
where £ is the bulk viscosity. The equation for a becomes 

^ = 12nGH( - l -^^H 2 . (11.43) 
a 2 

Thus not only w < —1/3 could drive acceleration, a viscosity ( > will also be able to 
drive acceleration. Meng, Ren and Hu | p93| , |294|| pointed out that the viscous fluid may 



give a unified description of dark energy and dark matter, under the name "dark fluid". 



A more general form of such an inhomogeneous equation of state is proposed in [[295 
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11.5. Particle physics models 

As is well known, non-relativistic particles behave as w = matter and relativistic 
particles behave as w = 1/3 radiation. However, when one considers particles with exotic 
properties, an equation of state other than to = or to = 1/3 may arise. 



For example, DeDeo ||296|| considered a Lorentz violating fermion field 

S = Jd 4 xe [^(iT^D^ - mty - (V^VOM , (H-44) 

where 6 M = (b, 0, 0, 0) is a time like vector, and other notations have been clarified under 
equation ( |11.34|) . It is argued that the Lorentz violation can be achieved by coupling the 
spinor to a condensate of ghost. 

The equation of state of the fermion particle can be calculated as 

w±(k)= r M^L— , (11.45) 
v ; 3[m 2 + (fc T b) 2 }' v ; 

where ± denotes the two helicities of the spinor, and k is the momentum of the particle. 

For the positive helicity particle, w becomes negative and a gas of such particles becomes 

a candidate for dark energy. 

As another approach, Bohmer and Harko [297|| reported that when a minimum length 
is considered, particle excitations could also arise as a component of dark energy. 

Mass variation could also affect the particles' equation of state. For example, Taka- 



hashi and Tanimoto [[298| showed that neutrinos with time varying mass could behave as 
dark energy. 

Alternatively, one could modify the particles' action from symmetry considerations. 
For example, Stichel and Zakrzewski ||299|1 derived an action as a dynamical realization 
of the zero-mass Galilean algebra with anisotropic dilational symmetry. It is shown that 
when the dynamical exponent is z = 5/3, the particles described by such an action (named 
darkons) behave as dark energy. 

By constructing a generalized dynamics for particles, Das et al. [|300|1 presented a new 
framework to generate an effective negative pressure and to give rise to a source for dark 
energy. 

Another particle dark energy model is holographic gas [175], which modifies the dis- 
persion relation of the particle such that the entropy of the gas is holographic. The model 
is already briefly reviewed in the section of Holographic principle. 

In addition, Lima, Jesus and Oliveira ||301|| suggested that the creation of cold dark 
matter can yield a negative pressure and is capable to accelerate the Universe. This model 
can mimics the ACDM model, thereby can provide a good fit to current cosmological data. 
For more studies about this model, see ||302| . 
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11.6. Dark energy perturbations 

A simple cosmological constant will not have perturbations simply because it is a con- 
stant. However, when the cosmological constant is replaced by other fields, fluids or par- 
ticles, dark energy will typically has perturbations, either originating from the fluctuation 
itself or originating from the clustering of other components via gravitational interaction. 

Here we follow the analysis in [|101|| (see also ||303|| ), where the analysis is performed 
in the Newtonian gauge, or equivalently, in terms of the gauge invariant Barden potentials 
3U4j|. The perturbed metric can be written as 



ds 2 = a(r)) 2 [(1 + 2$)cfy 2 - (1 - 2V)6 ij dx i dx j ] , (11.46) 

where after gauge fixing (B = and E = 0), $ and \& are the two remaining scalar 
perturbation variables. Also, perturbe the matter and dark energy components as 

p ->• p + dp = J2(P* + s Pi)i P^P + Sp = ^2(pi + S Pi ), (11.47) 

i i 

6 ee ik'ST^/ip + p), E ee [ST*, - ^S'jST^j /(p + p). (11.48) 

In Fourier space, the perturbed Einstein equations at linear order can be written as 
the following set of equations: 

-k 2 <!> -3H(<I>' + H$) = 4nGa 2 5p, (11.49) 

k 2 (y' + U$) =4nGa 2 (p + p)0, (11.50) 

k 2 

" + n(2W + + {2W + "H 2 )$ + —(*-$) = 4irGa 2 6p, (11.51) 

3 

k 2 (f -$) = 12nGa 2 (p + p)E, (11.52) 

where prime denotes derivative with respect to the conformal time 77, and % = a' /a. The 
above Einstein equations are consistent with a total energy conservation equation. For 
multiple fluids, these equations should be solved together with the equations of motion for 
these fluids to get a closed set of equations. 

For shear- less fluids (such as perfect fluid and scalar fields), the perturbation E van- 
ishes and $ = <&. Then the Newtonian potential $ can be solved as functions of the matter 
energy densities and velocities from ( pL 1 .49| ) and (|11.50| ) as 

3U(p + p)0~ 



k z <5> = -AitGa' 



5p+ *= 



(11.53) 
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The solution of the above perturbation equations depend on detailed models. When w 
crosses the phantom divide, some equations becomes singular and special care are needed 
for proper treatments ||305| , p06|| . In addition, dark energy perturbations also affect the cold 



dark matter power spectrum at large scale ||307|| . This can be used to distinguish dark 
energy models [|3U8|. 



12. The theoretical challenge revisited 



Before getting to the part of experiments and fitting, it is helpful to pause and revisit 
the theoretical challenge. We are having too many dark energy models. On the other 
hand, are we close to a solution to the dark energy problem? 

Let us review the questions discussed by Polchinski ||309|| . Firstly, Polchinski reminds 
us that vacuum will gravitate, at least in a local theory of gravity. 





in vacuum 



Fig. 2: An electron self energy diagram, in the presence of an atom (left panel) 
and in the vacuum (right panel). In a local field theory, graviton g can not tell the 
difference because it only feels the electron. It is not aware of the fields attached 
to the electron. 



The reason is that the gravitational coupling of vacuum energy is already locally 
measured. This comes from the accurate measurement of equivalence principle. Aluminum 
and platinum have the same gravitational mass to inertial mass ratio, up to an experimental 
error of 10 -12 . On the other hand, the electronic loop (see the left panel of Fig. 2) 
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contributions to their mass are about 10 -3 and 3 x 10 -3 to the rest energy of aluminum 
and platinum respectively. Thus up to an accuracy of 10~ 9 , gravity can feel the energy 
from the electronic vacuum loop diagram. 

On the other hand, the left panel and the right panel of Fig. 2 do not look locally 
different. As long as we have a local field theory, gravity can not couple to one but not 
the other. Thus one can not simply ignore the gravitational coupling to vacuum energy in 
a serious consideration of dark energy models. 

Secondly, the universe is not born empty. Instead it has a thermal history. At least 
we expect the abundance of helium and other light elements are explained by the big 
bang nucleosynthesis (BBN). Then the early universe should at least as hot as MeV in 
temperature. How to tune the current value of dark energy at that early time? Or if 
a tracking mechanism operates, how does the mechanism know the energy density will 
eventually fall to a small value instead of a large one? 

Moreover, if we expect the universe originates from a higher energy state, how does 
a solution of dark energy know we will eventually settle down in the present vacuum, not 
the SU(2) x U(l) invariant one or any other one? 

Thirdly, for modifications of gravity or matter on small scales, how does this mod- 
ification affect large scales while passing the solar system tests; On the other hand, for 
modifications of gravity or matter on larger scales, how can such a mechanism avoid the 
universe explode or collapse at small scales? For example, for a UV cutoff ~ lOOGeV, the 
curvature of the universe will be only one meter! 

After these considerations, Polchinski turned to anthropic principle. However, as we 



have discussed in subsection |5.3| , the anthropic principle itself has as many problems as any 
of the other models, if not a lot more. Polchinski also said, "The cosmological constant 
is nonzero, therefore we can calculate nothing." Anthropic principle not only has trouble 
itself, but brings trouble to the whole field of theoretical physics. 

With a great number of principles, ideas, mechanisms and models at hand, we are still 
not able to find a model which can answer all these questions, which at the same time is 
theoretically clearly derived. Thus we are either still very far away from a solution to the 
dark energy problem, or the solution hides behind one of the known mechanisms above, 
but a lot more details need to be understood. Very probably, the status of theoretical 
study needs to be substantially driven by future experiments, as we shall review below. 
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13. Cosmic probes of dark energy 

The most common approach to probe dark energy is through its effect on the expansion 
history of the universe. This effect can be detected via the luminosity distance (1l{z) and 
the angular diameter distance Da(z). In addition, the growth of large-scale structure can 
also provide useful constraints on dark energy. 

Theoretical models and observational data can be related through the x 2 statistic 0. 
For a physical quantity £ with experimentally measured value £ bs, standard deviation a^, 
and theoretically predicted value £ t h, X 2 is given by 

2 (£obs — £th) 2 / 1 q i \ 

Different cosmological observation gives different x|, an d the total x 2 is the sum of all xfs, 
i.e. 

X 2 = J2xl (13-2) 

By minimizing the total x 2 5 the best-fit model parameters can be determined. Moreover, 
by calculating Ax 2 = X 2 — Xmin' one can a ^ so pl°t the la, 2a, and 3a confidence level 
(CL) ranges of a specific model. Statistically, for models with different n p (denoting the 
number of free model parameters), the la, 2a and 3a CL correspond to different Ax 2 - In 
Table I, we list the relationship between n p and Ax 2 - In addition, the x 2 statistics are 



often performed by using the Markov Chain Monte-Carlo technology [|313| , |314 
Table I : Relationship between n p and Ax 2 



Tip 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Ax 2 (la) 


1.00 


2.30 


3.53 


4.72 


5.89 


7.04 


8.18 


9.30 


10.42 


11.54 


Ax 2 (2a) 


4.00 


6.18 


8.02 


9.72 


11.31 


12.85 


14.34 


15.79 


17.21 


18.61 


Ax 2 (3a) 


9.00 


11.83 


14.16 


16.25 


18.21 


20.06 


21.85 


23.57 


25.26 


26.90 



In this section, we will review some mainstream cosmological observations, introduce 
the basic principles of these observations, and describe how they are introduced into the 
X 2 statistics. 



13.1. Type la supernovae 



16 Another alternative is the median statistic. For more details, see [310,311,312] 
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Fig. 3: Discovery data: Hubble diagram of SNIa measured by the Supernova 
Cosmology Project and the High-z Supernova Team. Bottom panel shows residuals 
in distance modulus relative to an open universe with Q m = 0.3 and Qa = 0. From 
@, based on §§. 



Type la supernovae (SNIa) is a sub-category of cataclysmic variable stars that results 
from the violent explosion of a white dwarf star. A white dwarf star can accrete mass 
from its companion star; as it approaches the Chandrasekhar mass, the thermonuclear 
explosion will occur [[315| . Therefore, SNIa can be used as standard candles to measure 
the luminosity distance <1l(z) ||316| , |317| , |318|| , and thus provides a useful tool to measure the 
expansion history of the universe. In 1998, using 16 distant and 34 nearby SNIa from the 
Hubble space telescope (HST) observations, Riess et al. [|I| first discovered the acceleration 
of expanding universe. Soon after, based on the analysis of 18 nearby supernova (SN) 
from the Calan-Tololo sample and 42 high-redshift SN, Perlmutter et al. confirmed the 
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discovery of cosmic acceleration . The discovery of the universe's accelerating expansion 
(see Fig. 3) was another big surprise since Edwin Hubble discovered the cosmic expansion 
in 1929. Because of this great discovery, Adam Riess, Brian Schmidt, and Saul Perlmutter 
win the Nobel prize in physics 2011. For a more detailed history about the discovery of 
cosmic acceleration, see ||319|| . 



In recent years, the surveys of SNIa has drawn more and more attention [10,320,321 



Many research groups focused on this field, such as the Higher-Z Team [|322| , [323|| , the Su- 
pernova Legacy Survey (SNLS) HQ, the ESSENCE (denoting "Equation of State: Su- 
pErNovae trace Cosmic Expansion") [ |326|| , the Nearby Supernova Factory (NSF) [[527), the 



Carnegie Supernova Project (CSP) ||328|| , the Lick Observatory Supernova Search (LOSS) 
52J , and the Sloan Digital Sky Survey (SDSS) SN Survey PDj , etc. In 2008, the Su- 



pernova Cosmology Project (SCP) provided a framework to analyze these SNIa datasets 
in a homogeneous manner. From 414 SNIa samples, they selected 307 high-quality SNIa 
composing the "Union" dataset p31| |. In 2009, the Center for Astrophysics (CfA) SN 
project combined 90 low-redshift SNIa samples with the Union dataset and obtained the 
"Constitution" sample [ |332| . In 2010, the SCP updated their SNIa sample, and released 
the "Union2" dataset ||333|| , which consisting of 557 SNIa. Moreover, in a latest work [ [33 p , 
the Union2.1 SNIa dataset, which consisting of 580 SNIa, was released. 

To constrain dark energy by using the SNIa data, the absolute magnitude of SNIa must 
be determined first. Since the detailed mechanism of SNIa explosions remains uncertain 
||336|| , SNIa are not intrinsically standard candles, with a la spread of order 0.3 mag in 
peak S-band luminosity ||337j| . Fortunately, in 1992 Phillips found that SNIa has a clear 
correlation between their intrinsic brightness at maximum light and the duration of their 



light curve [338]. This so-called "Phillips correlation" was then used to turn SNIa into 



standard candles ||339|| . In 2004, utilizing "stretch" duration-magnitude correction, Kim 
et al. [ |335| reduced the dispersion on SNIa peak magnitude into only 0.10-0.15 mag. 
Fig. 4 shows a comparison of SNIa light curves before and after applying the "stretch" 
duration-magnitude correction. 

After determining the absolute magnitude of the SNIa, one can obtain the observa- 
tional distance modulus 

M, (13.3) 



Z^obs 



m 



where m is the apparent magnitude, and M is the absolute magnitude. On the other hand, 
the theoretical distance modulus can be calculated as 



IJ>th{zi) = 5 log 10 d L (zi) + 25, 



(13.4) 
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Fig. 4: Top panel: B-band light curves of the Calan-Tololo SNIa sample before any 
duration-magnitude correction. Bottom panel: Same light curves after applying the 



"stretch" duration-magnitude correction of Ref. [335 



and the luminosity distance cIl{z) is 

l + Z „ /„ rrrz — r f* dz' 



d L ( Z ) = i+ " f k (H 0V ^h\ r 

-noV «fc v Jo 



where 



(13.5) 



sin(x) if fi fc < (fc = 1); 

/ fc (x) = ^x if O fc =0 (fc = 0); (13.6) 
sinh(x) if Q k > (fe = -1). 

The x 2 f° r the SNIa data is 

.2 //,N [^obs(^») - //thfe; Q)} 2 . . 

x5jvW = 2^ ^ > ( 13 - 7 ) 

i=l * 

where 9 denotes the model parameters, fiobs(zi) and Oi are the observed value and the 

corresponding la error of distance modulus for each SNIa, respectively. It should be 
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mentioned that due to the uncertainty in the absolute magnitude of a SNIa, the degeneracy 
between the Hubble constant and the absolute magnitude implies that one cannot quote 
constraints on either one. Thus, when dealing with SNIa data, people often analytically 
marginalize the nuisance parameter Hq ||3-10|| . 

It should be stressed that the Eq. (|13.7| ) only includes the statistical errors of SNIa, 
while the systematic errors of SNIa are ignored. The systematic errors come from various 
factors including the errors in the photometry, the calibration, the identification of SNIa, 
the selection bias, the intrinsic variation of physical properties of SNIa, the host-galaxy 
extinction, the gravitational lensing, and so on. Currently, the systematic errors in the 
SNIa data have been comparable with the statistical errors, thus they should be considered 
seriously. To include the effect of systematic errors in the analysis, a prescription for using 



the Union2 compilation has been provided in ||341|| . The key of this prescription is a 
557 x 557 systematics covariance matrix, Csn, which captures the systematic errors from 
SNIa (This covariance matrix with systematics can be downloaded from ||341|| ). Utilizing 
Csn, one can calculate the following quantities 



A(9) = (^ bs - = 0, 9))(C^h(nf s - /4 h (Mo = 0, 9)), 

557 

B{9) = (CsNhM hS - A*?(A*o = 0,0)), 



i=l 



557 

where ,u° bs = fi bs(zi), n\ h = (M;h(zi)- The x 2 for the SNIa data is given by 

B{9) 2 



(13.8) 
(13.9) 

(13.10) 



X SN {9) = A(9) - 



C 



;i3.ii) 



In a latest work ||342|| , by using the SNLS3 sample of 472 SN, Conley et al. gave a 
most careful and detailed study about the systematic uncertainties of SNIa. After taking 
into account the systematic uncertainties of SNIa, the x 2 function can be written as 



2 - ArrF • CT 1 • Arrt 



XSN 



(13.12) 



where C is a 472 x 472 covariance matrix capturing the statistic and systematic uncer- 
tainties of the SNIa sample, Arrt = Mb — "t m od is a vector of model residuals, me is the 
rest-frame peak B band magnitude of an SN, and m mo d is the predicted magnitude of the 
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SN given by the theoretical model and two other quantities (stretch and color) describing 
the light-curve of the particular SN. The model magnitude m mo d is given by 

m mod = 5 log 10 D L (z heh z cmh ) - a(s - 1) + f3C + M, (13.13) 

where Dl is the Hubble-constant free luminosity distance, z C mh and Zhei are the CMB 
frame and heliocentric redshifts of the SN, s is the stretch measure for the SN, and C is 
the color measure for the SN. a and /? are nuisance parameters which characterize the 
stretch- luminosity and color-luminosity relationships, respectively. M. is another nuisance 
parameter representing some combination of the absolute magnitude of a fiducial SNIa 
and the Hubble constant. The total covariance matrix C in Eq. (|13.12| ) captures both the 



statistical and systematic uncertainties of the SNIa data. One can factor it as [|342 



C = D stat + C s tat + C sys , (13.14) 

where D stot is the purely diagonal part of the statistical uncertainties, C sta t is the off- 
diagonal part of the statistical uncertainties, and C sys is the part capturing the systematic 
uncertainties. It should be mentioned that, for different a and (3, these covariance matrices 
are also different. Therefore, one has to reconstruct the covariance matrix C for the 
corresponding values of a and (3. Here we do not describe these covariance matrices one 
by one. One can refer to Refs [[342| , |343|| and the public code ||344|| for more details about 



the explicit forms of these covariance matrices and the calculation of Xsn- 

In addition, some interesting methods are also proposed to reduce the systematic 



errors in the SNIa data. For instance, Wang and colleagues ||345|| developed a consistent 
frameworks for the flux-averaging of SNIa to reduce the effect of weak lensing on SNIa 
data. Currently, the systematic errors in the SNIa observations is the major factor that 
confining their ability to precisely measure the properties of dark energy. To enhance 
the precision of SNIa data, improvements on the photometric technique, as well as better 
understandings of the dust absorption and the SN explosions, are needed. For more details 
on SNIa observation and its cosmological applications, see [|3|, |M6| and references therein. 

13.2. Cosmic microwave background 
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Fig. 5: The WMAP 7-year temperature power spectrum, along with the temper- 



ature power spectra from the ACBAR [347| and QUaD [348] experiments. The 
solid line shows the best-fitting ACDM model to the WMAP data, corresponding 
to fU=0.738. From [§49|. 



CMB is the legacy of the cosmic recombination epoch. It contains abundant infor- 
mation of the early universe. In 1964 CMB was firstly detected by Penzias and Wilson 
350|| , who received the Nobel Prize in Physics 1978. Their work provided strong evidence 



that supports the Big Bang theory of the universe |351|1 - In 1989, the first generation of 
CMB satellite, the Cosmic Background Explorer (COBE), was launched. It discovered 
the CMB anisotropy for the first time [|352| and opened the era of the precise cosmol- 
ogy. Two of COBE's principal investigators, Smoot and Mather, received the Nobel Prize 
in Physics 2006. In 1999, the TOCO, BOOMERang, and Maxima experiments [p53 | 
firstly measured the acoustic oscillations in the CMB anisotropy angular power spectrum 
|33^ , |335l , p6|j3^ , |358| , |359| , p60| . In 2001, the second generation of CMB satellite, the Wilkin- 
son Microwave Anisotropy Probe (WMAP) ||361j , |362|| , was launched. It precisely measured 
the CMB spectrum (for latest results of WMAP, see Fig. 5) and probed various cosmo- 
logical parameters with a higher accuracy. Recently, the Planck satellite, as the successor 
to WMAP, was launched in 2009. The early results of Planck have been released recently 

The positions of the CMB acoustic peaks lie on the expand history from the decoupling 
epoch to the present epoch, and contains the information of dark energy thereof. Two 
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distance ratios are often used to constrain dark energy. The first distance ratio is the so- 
called "acoustic scale" I a, which represents the CMB multipole corresponds to the location 
of the acoustic peak. It can be calculated as 

iA = (l + *)^^. (13-15) 

Here z* denotes the redshift of the photon decoupling epoch, whose fitting formula is given 
by |5| 

z* = 1048[1 + 0.00124(fi 6 /i 2 )- - 738 ] [1 + gi (n m h 2 ) 92 ] , (13.16) 

where 

_ 0.0783(O b /i 2 )- - 238 0-560 
91 l + 39.5(0 6 /i 2 ) - 763 ' 92 l + 21.1(ft 6 /i 2 ) 1 - 81 ' 1 ' 

Da(z) is the proper angular diameter distance 

1 fk H oV\^k\ Jo 7T(F) 
d a{z) = jr 1 - - — 1 JJ , (13.18) 

and r s is the comoving sound horizon size 

fUec 1 pl/(l + z) 1 

r s (z)= Cs (i + Z )dt = — tt / ; „ ^ , • (13-19) 



o 



Jo a 2 H(a) y/1 + (30 6 /40 7 )a ' 



One can calculate Xcmb = (^ bs ~^a) 2 I ° 2 A t° include the CMB data into the x 2 statistics. 
The second distance ratio is the so-called "shift parameter" R which takes the form |]36U 



R(z*) = y/f^H Q (l + z*)D A {z+). (13.20) 

One can also calculate Xcmb = (-^obs — Rth) 2 /&R to reflect the contribution of the CMB 
data. 

Table II : Distance Priors from WMAP 7-year Fit 





7 — year ML a 


7 — year Mean b 


Error, a 




302.09 


302.69 


0.76 


R(z*) 


1.725 


1.726 


0.018 




1091.3 


1091.36 


0.91 
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"Maximum likelihood values; b Mean of the likelihood. 
Here we list the values of Ia(z*), R(z*) and z* obtained from the WMAP 7-year 
observations (Table 9 of |49|) in Table II. One can use these data to calculate the Xcmb 



defined above. Moreover, the Xcmb can t> e calculated as 

Xcmb = ( x i hh ~~ x i )(^ oy cMB)ij( x j a 



(13.21) 



where x\ = (I a, R, z*) is a vector, and (CovQ^g) is the inverse covariance matrix. For the 
WMAP 7- year observations , the inverse covariance matrix takes the following forms 

2.305 29.698 -1.333 
6825.270 -113.180 
3.414 

In addition, the presence of dark energy also affects the large scale anisotropy of the 
CMB through the Integrated Sachs Wolfe (ISW) effect |364|j , which is not captured by Eq. 
( |13.21| ). The ISW effect is caused by the variation of the gravitational potential during 
the epoch of the cosmic acceleration. It provides independent evidence for the existence 



of dark energy [|365|J366| . The ISW effect can be detected through the cross correlation 



between the CMB and LSS ||367|| . In ||368|| , Ho et al. reported a 3.7cr detection of ISW 
by cross-correlating the SDSS LSS observations with the WMAP CMB anistropies results. 
For more details on CMB observation and its cosmological applications, see ||369| , [370| , |371 1 
and references therein. 

13.3. Baryon acoustic oscillations 

Baryon acoustic oscillation (BAO) refers to an overdensity or clustering of baryonic 
matter at certain length scales due to acoustic waves which propagated in the early uni- 
354jj372| , |373|| . Similar to SNIa, which provides a "standard candle" for astronomical 



verse 



observations, BAO provides a "standard ruler" for length scale in cosmology to explore 
the expansion history of the universe ||374j| . The length of this standard ruler (~150 Mpc 
in today's universe ||375|| ) corresponds to the distance that a sound wave propagating from 
a point source at the end of inflation would have traveled before decoupling. BAO has a 
characteristic imprint on the matter power spectrum ||355|,|356|1 . So it can be measured at 
low redshifts z < 1 through the astronomical surveys of galaxy clusters [|376| . In addition, 
BAO scales can also be measured through 21 cm emission from reionization, which pro- 



vides abundant information about the early universe at high redshifts 1.5 < z < 20 [377 



The apparent size of the BAO measured from astronomical observations then leads to the 
measurements of the Hubble parameter H(z) and the angular diameter distance D a(z) 
|37|. 
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Fig. 6: The BAO power spectrum for the full luminous red galaxy (LRG) and main 
galaxy samples, measured by the SDSS-II survey. The solid curves correspond to 
the linear theory ACDM fits to WMAP3. The dashed curves include the nonlinear 
corrections. Results from LRG and main galaxy samples are consistent with each 
other, and all provide a confirmation of the predicted large-scale ACDM power 
spectrum. When combined with the WMAP3 result, the LRG data yield to Qa = 

0.761±Hi8 and w = -°- 941 io'.ioi- From B- 



The BAO measurement does not require precision measurements of galaxy mag- 
nitudes, nor did it require that galaxy images be resolved; instead, only their three- 
dimensional positions need to be determined. So the BAO observations are less affected 
by astronomical uncertainties than other probes of dark energy, although it still suffers 
from some systematic uncertainties, such the effects of non-linear gravitational evolution 
and redshift distortions of clustering J2"^| , |380| . The acoustic signature of BAO has already 
been obtained in the galaxy power spectrum at low redshift [ 381 , 374f . For examples, 
the Two-degree-Field Galaxy Redshift Survey (2dFGRS) ||382|| , which was operated by the 
Anglo-Australian Observatory (AAO) between 1997 and 2002, had made public their BAO 
data in 2003. In addition, using a dedicated 2.5-m wide-angle optical telescope, the Sloan 
Digital Sky Survey (SDSS) |[383|| , was launched in 2000. In 2006, the survey entered a new 
phase, the SDSS-II survey. In 3 years, it completes the observations of a huge contiguous 
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region of the northern skies. The BAO power spectrum measured by the SDSS-II survey 
|379H is shown in Fig. 6. In 2009. the SDSS-III survey begun. In January 2011, SDSS-III 



publicly release its eighth Data Release (DR8) ||384[ , which is the latest BAO dataset to 
data. 

From measurements of the galaxy clusters, the BAO scales in both transverse and 
line-of-sight directions are obtained; they correspond to the quantities r(z)/r s (z) and 
r s (z)/H(z), respectively. In addition, three characteristic quantities of BAO, including 
the A parameter, Dy (0.35) / Dy (0.2), and r s (zd)/Dy(z), are often used to constrain dark 
energy parameters. In the following we will provide a rough introduction to these three 
quantities. 



The A parameter is defined as ||374 ] 



A 



th 



(H(z b )/H Q ) 



1 



ZbV\^k\ 



fk ( h vfn* 



dz' 



[13.22) 



where z b = 0.35 is the redshift of the SDSS luminous red galaxies (LRG). The SDSS BAO 

0.469(ti s /0.98)-°- 35 ± 0.017, where the WMAP 7-year 
0.963. One can calculate Xbao — 

(A ohs - A th ) 2 /a 2 A to 



measurements ||374|| gave A ^ s 
results gave a best-fit value n s 
reflect the impact of the BAO data. The A parameter is considered independent of dark 
energy models, and have been widely used in the literature. 

Next, let us turn to the quantity Dy (0.35)/ Dy (0.2). Since the current observations 
are not sufficient enough to measure the BAO scales in both transverse and line-of-sight 
directions independently, alternatively people construct an effective distance ratio Dy (z) , 
which is defined as ||374j| 

1/3 



D v (z)= (l + zYD\(z) 



H(z) 



(13.23) 



In 2005, Eisenstein et al. [374] provided a constraint Dy (z) = 1370±64 Mpc at the redshift 
2=0.35. In 2009, the SDSS DR7 sample fl385[ gives Dy (0.35) j Dy (0.2) = 1.736 ± 0.065. 
One can use this quantity to reflect the impact of the BAO data. 

Lastly, we introduce the quantity r s (zd)/Dy(z). From the SDSS and the 2dFGRS, 
one can extract a quantity r s (zd)/Dy(z) at given z, where Zd denotes the redshift of the 



drag epoch, whose fitting formula is proposed by Eisensten and Hu ||358 

1291(O m /i 2 ) - 251 



Zd 



1 + 0.659(Q m /i 2 ) - 828 
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[i + b 1 (n b h 2 ) b2 ] 



(13.24) 



where 



h = 0.313(O m /i 2 )-°- 419 [1 + 0.607(O m /i 2 ) - 674 ] , b 2 = 0.238(Q m /i 2 ) - 223 . (13.25) 

It is widely believed that r s (zd)/Dy(z) contains more information of BAO than the pre- 
vious two quantities. 

In addition, one can also use the covariance matrix method to construct the x 2 function 
of the BAO data 

X 2 BA o = A Pi [Cov^ AO ( Pi , Pj )]A Pj , Api = pf ata - Pl . (13.26) 

For the latest BAO data from SDSS DR7 [ |385| , P \ = do. 2 and P i = do. 35. The covariance 
matrix is 

„ _i / 30124 -17227 \ 
^ovbao - ^ 86977 J , 

and the vales of P i are 

ta = rf data = g.1905, p^ ta = dg$» = 0.1097. (13.27) 
There are also some other methods capturing the information of BAO from the SDSS data, 



sec 



374| , |386|| and references therein for details. 



13.4- Weak lensing 

Weak lensing (WL) is the slight distortions of distant galaxies' images, due to the 
gravitational bending of light by structures in the Universe (see Fig. 7). Utilizing the 
WL effect, the distribution of dark matter and its evolution with time can be measured, 
thus providing a useful tool to probe dark energy through its influence on the growth of 
structure. 

In the 1990s, WL around individual massive halos was measured [|387| , |388|1 . Soon after, 



WL by LSS was detected by four research groups in 2000 [|389| , |390| , [391| , |392|1 ; from then on 



WL has grown into an increasingly accurate and powerful probe of dark matter and dark 
energy [[3931 • The current survey project of WL is the Canada- France-Hawaii legacy survey 
(CFHTLS) [ |394| , covering ~ 170 square degrees. 



The effect of WL on the distant sources can represent on the distortions in the shapes, 
sizes and brightnesses ||395|| . Current studies mainly focussed on the change in shapes 



(termed as "cosmic shear"), which can be more easily and precisely measured compared 
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Fig. 7: Cosmic shear field (white ticks) superimposed on the projected mass dis- 
tribution from a cosmological iV-body simulation: overdense regions are bright, 
underdense regions are dark. Note how the shear field is correlated with the fore- 
ground mass distribution. From [1C]. 



with the changes in the size and brightness. The lensing effect on the shape of the galaxies 
is by general ~0.01. This effect is much smaller than the typical deviation in the galaxy 
shape, which is about 0.3 — 0.4. So a large number of galaxies are required to detect the 
cosmic shear signal with enough precision. 

The cosmic shear analysis is a widely used method to relate the WL data with the 
dark energy. The cosmic shear field is the weighted mass distribution integrated along 
the line of sight. The Fourier transformed counterpart is the shear power spectrum P K . 
Here k is the "convergence" , which means a magnification of the source of 1 + 2 k at the 
linear order. P K can be inferred from the 3D distribution of matter, and it is related to 
the matter power spectrum P(k, r) through 

Here r(z) = is the comoving distance, n(z) is the mean redshift distribution of the 

source galaxies normalized to unity, rjj is the comoving horizon distance (i.e. the depth 
of the survey). 



dr n(r 



fk(r') 



;i3.28) 
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For statistical analysis of cosmic shear, it is common to use 2-point correlation func- 
tions [|396| . For example, a convenient way is to describe the shear field in terms of 



E/B-modes correlation functions ||397|| 



where 



and 



um = m±m , ue) = M<o-m i (13 . 29) 



poo 11 POO 11 

e+W = y ^lPn(l)MW), = J ^IP K (1)M19), (13-30) 



£>(9) = + 4 dO'^jp- - 126 2 J 



/£-(*') mn2 r de >^)_ (13 31) 



Here Jq and J4 are the zeroth and forth order Bessel functions of the first kind, respectively. 
The shear correlation functions can then be compared with the measurements directly. The 



B-mode correlation functions are expected to be very small [398], so it is common to use 



the E-mode correlation functions to construct the weak lensing % 2 function 

XL = (Cs(0O - m i) (Cov" b %- ((Eidj) - mj ). (13.32) 

Systematic errors in weak lensing measurements arise from a number of sources, in- 
cluding incorrect shear estimates, uncertainties in galaxy photometric redshift estimates, 
intrinsic correlations of galaxy shapes, and theoretical uncertainties in the mass power 
spectrum on small scales ||399| , |400|| . Fortunately, future WL surveys have ability to inter- 
nally constrain the impact of such effects |2"4| , f4Ul| , |4U2[ j. For more details about the WL 
observation and its applications on the probe of dark energy, see |7| J395| , [I03| and references 
therein. 

13.5. Galaxy clusters 

Galaxy clusters (CL) are the largest gravitationally bound objects in the universe. 
They typically contain 50 to 1,000 galaxies and have a diameter from 2 to 10 Mpc. CL can 
be detected through the following approaches ||4U5||: (1) optical or infrared imaging and 



spectroscopy (see Fig. 8). In comparison with optical surveys, infrared searches are more 
useful for finding higher redshift clusters. (2) X-ray imaging and spectroscopy (see Fig. 8). 
CL with active galactic nucleus (AGN) are the brightest X-ray emitting extragalactic ob- 



jects, so they are quite prominent in X-ray surveys. (3) Sunyaev-Zel'dovich effect [ j373| . 



The hot electrons in the intracluster medium scatter radiation from the cosmic microwave 



99 



Fig. 8: Massive galaxy cluster 2XMM J083026+524133 detected by the X-ray 
telescope XMM-Newton (the green contours) and its the optical image observed 
by the Large Binocular Telescope (LBT) in the Arizona desert. From [404]. 



background through inverse Compton scattering. This produces a "shadow" in the ob- 
served CMB at some radio frequencies. (4) gravitational lensing. CL bend the light from 
distance galaxies and distort the observed images. The observed distortions can be used 
to detect the masses of clusters. 

In principle, the number density of cluster-sized dark halos n(z, M) as a function of 



redshift z and halo mass M can be accurately predicted from N-body simulations [[406 
Comparing these predictions to large area cluster surveys can provide precise constraints 
on the cosmic expansion history P07| , |408| , |409 



In a survey that selects clusters according to some observable O with redshift- 
dependent selection function f(0,z), the redshift distribution of CL is given by |L(J 



d 2 N(z) _ r\z) f°° w ^ [°° ^dn(z 

dzdtt ~ Wz) 



J™ f(0, z)dO J™ p(0\M, z) < ^dM. (13.33) 



Here N is the number of CL, O is the solid angle, r(z) is the comoving distance, H(z) 
is the Hubble parameter, dn(z)/dM is the space density of dark halos in comoving co- 
ordinates, and p(0\M,z) is the mass-observable relation, the probability that a halo of 
mass M at redshift z is observed as a cluster with observable property O. The utility of 
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this probe depends on the ability to robustly associate cluster observables such as cluster 
galaxy richness, X-ray luminosity, Sunyaev-Zeldovich effect flux decrement, or weak lens- 
ing shear with cluster mass [|T0| , f405| . As seen in this equation, the sensitivity of cluster 



counts to dark energy arises from two factors: cosmic expansion history, r 2 (z)/H(z) is 
the comoving volume element that contains information of the cosmic expansion history; 
growth of structure, dn(z)/dM depends on the evolution of density fluctuations, and the 
cluster mass function is also determined by the primordial spectrum of density perturba- 
tions. One can see ||410|| and references therein for more details about the surveys of CL 



and their applications on the probe of dark energy. 

13.6. Gamma-ray burst 

Gamma-ray bursts (GRB) are flashes of gamma rays associated with extremely ener- 
getic explosions in distant galaxies. They are the most luminous electromagnetic events in 
the universe. 

GRBs were first detected in 1967 by the U.S. Vela satellites. They can be classified 
into "long GRBs" (their durations are longer than 2 seconds) and "short GRBs" (their 
durations are shorter than 2 seconds). At the beginning, most astronomers believed that 
GRBs originate from inside the Milky Way galaxy, only Paczynski insisted that GRBs 
originate from external galaxy ||411|| . In 1997, an X-ray astronomy satellite BeppoSAX 



detected the "afterglow" of GRB 970228 jH^I , and verified that GRBs indeed originate 



from external galaxy. This discovery opened up a new era in the history of GRBs studies. 

Many satellites had been launched to probe GRBs in the past decade. The Swift 
mission |41ij| j was launched in 2004. This satellite is equipped with a very sensitive gamma 



ray detector as well as on-board X-ray and optical telescopes, which can be rapidly and 
automatically slewed to observe afterglow emission following a burst. In 2008, the Fermi 
|414j| was launched. Its main instruments include the Large Area Telescope (LAT) 



mission 



and the Gamma-Ray Burst Monitor (GBM); the former is used to perform an all-sky 
survey, and the latter is used to detect sudden flares of gamma-rays. Meanwhile, on the 
ground, numerous optical telescopes have been built or modified to incorporate robotic 
control software that responds immediately to signals sent through the Gamma-ray Burst 
Coordinates Network (GCN) | }415|| . This allows telescopes to rapidly repoint towards a 



GRB, often within seconds of receiving the signal and while the gamma-ray emission itself 
is still ongoing. 
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Fig. 9: The Hubble diagram for 69 GRBs, out to the redshift of z > 6. The curve 
is the luminosity distance in a flat ACDM model with Q m =0.27 and w = —1. It 
can be seen that the observational data smoothly follows the curve. From [4161. 



GRBs have been proposed to be a complementary probe to the SNIa in the studies of 
dark energy. For example, in |[416|| , Schaefer presented 69 GRBs data points (see Fig. 9) 
over a redshift range of z = 0.17 to z > 6 with half the bursts having a redshift larger 
than 1.7, and showed that the GRB Hubble diagram is consistent with the existence of 
cosmo logical constant. Compared to the SNIa, an advantage of GRB is that the high 
energy photons in the gamma-ray band are almost immune to dust extinction. Moreover, 
the redshifts of observed GRBs are much higher than SNIa: there have been many GRBs 
observed at 1 < z < 8, whereas the maximum redshift of GRBs is expected to be 10 or even 



larger |[417|| . Therefore, GRBs are considered to be a promising probe to fill the redshift 
desert between the redshifts of SNIa and CMB. For more details on the GRB cosmology, 
see 



418 , 419 ] and references therein. 



However, there still exist some troubles in the application of GRB data to the probe 
of dark energy. A big problem is that since our knowledge on the mechanisms underlying 
the GRB emission is still limited, treating them as standard candles is still suspicious. 
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Another well-known problem is the so called "circularity" problem in the GRB calibration 
|419j , [420fl , mainly due to the lack of low-redshift GRBs at z < 0.1 which are cosmology- 
independent. To alleviate this problem, various statistical methods have been proposed. 



For examples, in [[42 If , Wang summarized the GRB data by a set of model-independent 
distance measurements and provided a convenient method to use the GRB data in cosmo- 
logical analysis. In |[422||, Liang et al. proposed a new method to estimate the distance 



modulus of GRBs by interpolating from the Hubble diagram of SNIa. Using this method, 
in 



423|| , Wei obtained 59 calibrated high-redshift GRBs from a total number of 109 long 



GRBs ||416| , [424|| and the Union2 SNIa sample ||333|| . For more statistical methods to cal 



ibrate the GRB data, see e.g. [ |419| , |420|j422| ]. Due to the lack of a large amount of well 
observed GRBs, the current GRB data are still not able to provide forceful constraint on 
dark energy. Therefore, although the GRBs are considered to be a promoting probe to fill 
the redshift desert between SNIa and CMB, there is still a long way to use them extensively 
and reliably to probe dark energy. 

13. 7. X-ray observations 

X-ray is an important observational branch of astronomy which deals with the study 
of X-ray emission from celestial objects. In 1962, utilizing a US Army V-2 rocket, Giacconi 
and colleagues firstly discovered cosmic X-ray source. In addition, they also discovered the 
existence of astronomical X-ray background [[425]. Due to his great contributions to this 
field, Giacconi received the Nobel Prize in Physics 2002. 

Since Giacconi' s great discovery, many X-ray satellites have been launched. The first 



orbiting X-ray astronomy satellite, the Uhuru satellite [4261, was launched in 1970. Soon 



after, in 1978, the first fully imaging X-ray telescope, the Einstein Observatory |427fl , was 



launched. In 1999, two important X-ray satellites, the Chandra observatory ||428|| and 



the XMM- Newton observatory ||429|| , were launched. The Chandra observatory has high 



space resolution (less than 1 arc-second) and a wide wave band (0.1-1 keV), while the 
XMM-Newton observatory has very high spectrum resolution. 
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Fig. 10: The 68.3% and 95.4% confidence constraints in the fl m -Q,A, using Chan- 
dra measurements of the X-ray gas mass fraction / gas data (red contours) . A non- 
flat ACDM model is assumed. Also shown are the independent results obtained 
from the CMB data (blue contours) and SNIa data (green contours). A combined 
analysis using all the three data sets yields a result of Qa = 0.735 ± 0.023 and 



-0.010 ± 0.011 at the 68.3 % CL (the inner, orange contours). From pC] 



In galaxy clusters, the X-ray emitting intracluster gas is the dominate component of 
baryonic mass content (it exceeds the mass of optically luminous material by a factor ~6 
||431|| ). In addition, galaxy clusters also have the dark matter content |432j , ^33|| . So the 
measurements of universal baryonic mass function in the clusters, /j, = f^/S^m, have been 
widely used to determine the cosmic matter fraction fi m [f434| . In late 1990s, Sasaki [[135 1 
and Pen p36|| proposed that the dependence of fb on the angular diameter distances to the 
clusters can also be used to constrain the geometry of the universe. In 2002, Allen et al. 
137] . r 138||439|1 carried out such a test and obtained fi m = 0.30 ± 0.04 using a small sample 



of X-ray luminous galaxy clusters. In 2008, based on an improved analysis using the data 
of 42 clusters from Chandra X-ray observations, Allen et al. [|430|1 found a detection of 
dark energy at 4a CL, with fl A = 0.86 ± 0.21 for a non-flat ACDM model (see Fig. 10). 

The results of Refs. ||437| , [438|1439| . |430|| are obtained from the measurements of the 
apparent evolution of the cluster X-ray gas mass fraction, / gas = M gas /M to t- As showed 
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in ||430||, the x 2 function for / gas data can be calculated as 



42 



EU gas Jt 



gas,ij 



(13.34) 



gas, i 



where the fitting formula of fg^ I)M (z) is given by 



ACDM, 
gas 1 



KAjb(z) ( O 
l + s(z) 



jACDM 
a A 

[Da(z) 



1.5 



(13.35) 



Here D A and L>^ CDM are the angular diameter distances to the clusters in the current 
test model and reference cosmologies, respectively. The factor b(z) = bo(l + a^z) is a 
"biased factor" with 0.65 < bo < 1.0, —0.1 < ab < 0.1; parameter s(z) = so(l + a s z) 
models the baryonic mass fraction in stars, and sq = (0.16 ± 0.05)h®Q, —0.2 < a s < 0.2. 7 
models non-thermal pressure support in the clusters (1.0 < 7 < 1.1). K is a "calibration" 
constant arises from the residual uncertainty in the accuracy of the instrument calibration 
and X-ray modelling (K = 1 ± 0.1). The factor A is 



A 



//iACDM\ V 

(irH 

V y 2500 / 



H(z)D A {z) 



[H{z)D a {z)} aci >m 



;i3.36) 



where r\ = 0.214 ± 0.022. One can refer to [430] for details about the origins of the 
parameters b(z), s(z), K, and A. 

Like the SNIa and BAO data, the / gas measurements can also probe the redshift- 



43011 . It has been shown that 



distence relation, with the dependence / gas oc cIl(z)Da(z) ' 5 
|430| , |440| , |441|| the / gas data is useful in breaking the degeneracy between some cosmological 
parameters such as O m , w and Hq when combined with other cosmological observations. 
In addition to / gas , measurements of the amplitude and evolution of matter fluctuations 
using X-ray observations have also been applied to probe dark energy. One can refer to 
g3|,|42[ for more details. 



13.8. Hubble parameter measurements 



In 1929 Hubble [|443|1 discovered a linear correlation between the apparent distances 
to galaxies D and their recessional velocities v 



(13.37) 
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where Hq is the so-called Hubble constant. This discovery provided strong evidence that 
our universe is in a state of expansion and opened the era of the modern cosmology. Soon 
after, people find that this Hubble's law is just an approximate formula, and Hq should 
be replaced by H(z), a function of the redshift z. As mentioned above, H(z) describes 
the expansion history of the universe, and plays a central role in connecting dark energy 
theories and observations. 

At the beginning, Hubble and Humason ||444|| measured a value for Hq of 500 
km/s/Mpc, which is much higher than the currently accepted value due to errors in their 
distance calibrations. Since the launch of the Hubble Space Telescope (HST) | )445|| , the 
value of Hq was estimated to be between 50 and 100 km/s/Mpc. For example, using the 
HST key project, Freedman et al. [f!46| obtained Hq = 72 ± 8 km/s/Mpc, and Sandage 
et al. |47| advocate a lower Hq = 62.3 ± 6.3 km/s/Mpc. In 2009, Riess et al. [@48| gave 
Hq = 74.2 ± 3.6 km/s/Mpc with a 4.8 % uncertainty. In a latest work [ 449|| , Riess et al. 
obtained Hq = 73.8 ± 2.4 km/s/Mpc, corresponding to a 3.3 % uncertainty. In the future 
with more precise observations from HST, Spitzer 1 45U j , Global Astrometric Interferometer 
for Astrophysics satellite (GAIA) ||451|| and James Webb Space Telescope (JWST) [|52 



an uncertainty of 1% in the Hubble constant will be a realistic goal for the next decade 

EH- 

The precise measurements of Hq will be helpful to break the degeneracy between some 



cosmological parameters ||453|| . For example, in ||454|| , Hu pointed out that a measurement 
of Hq to the percent level, when combined with CMB measurements with the statistical 
precision of the Planck satellite, offers one of the most precise measurements of dark energy 
EOS at z ~ 0.5. 



Table III : Values of Hq Measured in p55jg56| 



z 


0.1 


0.17 


0.27 


0.4 


0.48 


0.88 


0.9 


1.3 


1.43 


1.53 


1.75 


H(z) 


69 


83 


77 


95 


97 


90 


117 


168 


177 


140 


202 




12 


8 


14 


17 


60 


40 


23 


17 


18 


14 


40 



Table IV: Values of Hq Measured in [|57 
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z 


0.24 


0.34 


0.43 


H(z) 


79.69 


83.80 


86.45 


0~H,st 


2.32 


2.96 


3.27 


0~H,sys 


1.29 


1.59 


1.69 



The inferred H(z) with its statistical and systematical errors for each redshift slice. 



In addition to the direct measurements for the Hq, the precise measurements of H(z) 
are also useful in studying the cosmic acceleration. In 2005, Simon et al. §455|| measured 
the Hubble parameter H(z) at nine different redshifts from the differential ages of passively 



evolving galaxies. In 2009, Stern et al. |[456|| extended this dataset to eleven data points (see 
Table III). Soon after, by studying the clustering of LRG galaxies in the latest spectroscopic 



SDSS data releases, Gaztanaga et al. [j457| obtained three more H(z) data points (see 
Table IV). Utilizing these H(z) data, it is straightforwad to put constraint on dark energy 
parameters by calculating the corresponding % 2 as 



Xh 



14 

E 

1=1 



[H ohs (zi) - H(zi 



(71 



(13.38) 



As shown in ||458|| , the current H(z) data from direct measurements can provide valuable 
constraints on dark energy. In the future, with the developments in the observational tech- 
nique of LRGs, the H{z) measurements can provide useful complements to other cosmic 
observations [|459|J460|| ■ 



13.9. Cosmic age tests 

The cosmic age problem is a longstanding issue in cosmology and provides an impor- 
tant tool for constraining the expanding history of the Universe f461| • Before the great dis- 
covery that our universe is undergoing an accelerated expansion, the most popular SCDM 
model (i.e., a flat universe with Q m = 1) is always plagued by a longstanding puzzle: in 
this model the present age of the universe is to = ~ 9 Gyr, while astronomers have 
already discovered that many objects are older than 10 Gyr |f462| , |463|| . For example, based 
on the study of white dwarf cooling, a lower limit of cosmic age to = 12.7 ± 0.7 Gyr have 



been obtained [464]. The cosmic age problem becomes more acute if one considers the age 
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of the universe at a high redshift. For instance, a 3.5 Gyr-old galaxy 53W091 at redshift 
z = 1.55 and a 4 Gyr-old galaxy 53W069 at z = 1.43 are more difficult to accommodate in 
the SCDM model ||465|1 . Along with the discovery of accelerated expansion of the universe 
and the return of the cosmo logical constant A, the cosmic age problem has been greatly 
alleviated. The 7-year WMAP observations show that in the ACDM model the present 
cosmic age is £g bs = 13.75 ±0.11 Gyr. Besides, it is shown that the ACDM model can also 
easily accommodate galaxies 53W091 and 53W069 |[466|| . So the cosmic age problem is a 
"smoking-gun" of evidence for the existence of dark energy. 

However, the cosmic age problem has not been completely removed by the introduction 
of dark energy. By comparing photometric data acquired from the Beijing- Arizona- Taiwan- 



Connecticut system with up-to-date theoretical synthesis models, Ma et al. ||467|| obtained 
the ages of 139 globular clusters (GCs) in the M31 galaxy, in which 9 extremely old 
GCs are older than the present cosmic age predicted by the 7-year WMAP observations 
468H . In addition, the existence of high-z quasar APM 08279+5255 at z = 3.91 []469[| is 



also a mystery |470| , |471||472|| . Using the maximum likelihood values of the 7-year WMAP 
observations fi m = 0.272 and h = 0.704, the ACDM model can only give a cosmic age 
t = 1.63 Gyr at redshift z = 3.91, while the lower limit of this quasar's age is 1.8 Gyr. 
To accommodate these anomalous objects, some authors suggested that a lower Hq should 



be advocated [ |447| , f473|| , while some authors suggested that more complicated cosmological 
model should be taken into account 1 474 , 475 1. Therefore, the cosmic age puzzle still remains 
in the standard cosmology. 

In addition, one can also use the ages of old galaxies to perform the best-fit analysis 
on dark energy. In [ f455| , Simon et al. established these so-called "lookback time-redshift" 
(LT) data by estimating the age of 32 old passive galaxies distributed over the redshift 
interval 0.11 < z < 1.84 and the total age of the universe t° hs . The galaxy samples 
of passively evolving galaxies are selected with high-quality spectroscopy, and the method 
used to determine ages of galaxy samples indicates that systematics are not a serious source 
of error for these highredshift galaxies |[455|| . Utilizing these LT data, one can calculate the 
corresponding x 2 as 

= » + es^oM! (1 , 39) 

1=1 ± l o 

where t(z; 9) is the age of the universe at redshift z, given by 

r°° dz 

/ rr-TT^r. (13.40) 

Jz (l + z)H(z;9) 
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Here t (9) = t(0;6) is the present age of the universe. o\ = af + cr^ oba , where Oi is the 
uncertainty in the individual lookback time to the i galaxy of the sample, o" t ob a = 0.7Gyr 
stands for the uncertainty in the total expansion age of the universe (t° bs ), and r means 
the time from Big Bang to the formation of the object. 



13.10. Growth factor 



In addition to the expansion history of the universe, the growth of large-scale structure 
can also provide important constraints on dark energy and modified gravity. The growth 
rate of large scale structure is derived from matter density perturbation 5 = 5p m /p m in 
the linear regime that satisfies the simple differential equation ||476| 



6 + 2H5 -AirG effPm 5 = 0, 



(13.41) 



where the effect of dark energy is introduced through the "Hubble damping" term 2H5, 
and the effect of modified gravity is introduced via the effective gravitational "constant" 
G e ff. Eq. ( P-3.41 ) can be written in terms of the logarithmic growth factor / = dlnS/dlna 



df 
din a 



+ f + 
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W 



+ 2 / 



3G 
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2 G 



(13.42) 



There is no analytical solution to the Eq. ( |13.42|) , and much efforts have been paid 
to solve this equation ||477| , [178| , f479|| . The growth function / was shown to be well approx- 
imated by the ansatz IQgLgggggl 



/ = «; 



(13.43) 



where 7 is the so called "growth index" . Based on this ansatz, the approximate expressions 
of / for various gravitational theory have been obtained | 482 , 183 |. Besides, there has been 



much interest in exploring parameterization of the growth index as a function of the redshift 



I Hi- 



Table V : The Growth Factor Data 
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z 



fohs References 



0.15 


0.49 ±0.1 


485 


0.35 


0.7 ±0.18 


379 


0.55 


0.75 ±0.18 


486 


0.77 


0.91 ±0.36 


485 


1.4 


0.9 ±0.24 


487 


3.0 


1.46 ±0.29 


488 


2.125 - 2.72 


0.74 ±0.24 


489 


2.2 - 3 


0.99 ±1.16 


490 


2.4 - 3.2 


1.13 ±1.16 


490 


2.6 - 3.4 


0.99 ±1.16 


490 


2.8-3.6 


0.99 ±1.16 


490 


3.0-3.8 


0.99 ±1.16 


490 



In Table V, we list the growth data that was converted in the work of |[491|J492| , |485| , 
from either measurement of redshift distorsion parameter or from various power spectrum 
amplitudes from Lyman-a Forest data. The list of the respective original references is 
also given in the Table V. A caveat in using this data to constrain other cosmological 
models is that in various steps in the process of analysing or converting the data, the 
ACDM model was assumed. So if one wants to use the data to constrain other models, 
and in particular modified gravity models, one than should redo all the steps assuming 
that model, starting from original observations. For more studies concerning the growth 



factor and its application on the probe of dark energy, see e.g. ||491| , f492| , |493| , |194 
13.11. Other cosmological probes 



Sandage-Loeb test |[495| , |496|1 is a method which directly measures the evolution of the 
universe. The idea is to measure the drifts of the comoving cosmological sources caused 
by the cosmic acceleration/deceleration. For an object at redshift z s , after a period Ato 
(to stands for the cosmic age at our position), there would be small variation of its redshift 
due to the drift caused by the evolution of the universe 



Az 



a(t ) - a(t s ) 
a(t s ) 



At . 



13.44) 
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The variations of redshifts can be obtained by direct measurements of the quasar Lyman-a 
absorption lines at sufficiently separated epochs (e.g., 10-30 yrs). The Sandage-Loeb test 
is unique in its coverage of the "redshift desert" at 2 < z < 5, where other dark energy 
probes are unable to provide useful information about the cosmic expansion history. Thus, 
this method is expected to be a good complementary to other dark energy probes [[497 . 



Gravitational waves (GW) observations also have potential to make interesting con- 
tributions to the studies of dark energy. In 1986, Schutz [f!98| found that the luminosity 
distance of the binary neutron stars (or black holes) can be independently determined by 
observing the gravitational waves generated by these systems. If their redshifts can be 
determined, then they could be used to probe dark energy through the Hubble diagram 
499j| . This so-called GW "standard sires" has drawn a lot of attentions ||500|| . Besides, 



dark energy can also leave characteristic features on the spectrum of primordial gravita- 



tional waves H503-II , which may be detected by future ground-based and space-based GW 
detectors [|02| . 

In addition, through astronomical observations, the validity of general relativity can 
be tested from observations of solar systems, BBN, CMB, LSS, gravitational waves, and so 
on. In the following we will present a brief introduction to the tests coming from the solar 
system and BBN. For more details about the tests of gravity theories, see ||504j,|505|,|506 



and the references therein. 

Solar system tests of alternate theories of gravity are commonly described by the 
Parametrized Post-Newtonian (PPN) formalism [ 5Q7|| . In this formalism, the metric is 
written as a perturbation about the Minkowski metric 



ds 2 = -(1 + 2* - 2f3^ 2 )dt 2 + (1 - 2^)dx 1 



(13.45) 



here the potential \I/ = -GM/r for the Schwarzschild metric. The parameter {3 stands for 
the nonlinearity in the superposition law of gravity, and 7 describes the spacetime curvature 
induced by a unit mass. The parameter 7 is the most relevant PPN parameter for the 
modified gravity theories of interest. So far, 7 has been tightly constrained from various 
observational methods. The tightest constraint comes from time-delay measurements in 
the solar system, specifically the Doppler tracking of the Cassini spacecraft, which gives 
7 — 1 = (2.1 ± 2.3) x 10 -5 ||508|| . Other tests include the observed perihelion shift of 



Mercury's obit [|509| and light deflection measurement [|510 |, which can constrain 7 at the 
10~ 3 and 10~ 4 level, respectively (see Fig. 11). 
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Fig. 11: Measurements of the coefficient (f + 7)/2 from light deflection and time 
delay measurements. Its value in the Einstein gravity is unity. The arrows at the 
top denote anomalously large values from early eclipse expeditions. The Shapiro 
time-delay measurements using the Cassini spacecraft yielded on agreement with 
Einstein gravity to 1CP 3 percent, and VLBI (very-long-baseline radio interferome- 
try) light deflection measurements have reached 0.02 percent. Hipparcos denotes 



the optical astrometry satellite, which reached 0.1 percent. From 504] 



Moreover, solar system observations can also test the modified gravity models through 
their violations of the Strong Equivalence Principle (SEP). That will result in a difference 
in the free-fall acceleration of the earth and the moon towards the sun, named as the 
Nordtvedt effect |509|j . This effect is detectable in the Lunar Laser Ranging (LLR). Current 
LLR data have constrain PPN deviations from the Einstein gravity at the 10 -4 level [ p04 



In the near future, the bound is expected to be improved by an order of magnitude by the 



Apache Point Observatory for Lunar Laser-ranging Operation (APOLLO) project | 511 



The BBN happened when the universe was ~ 10-100 seconds old. The abundance of 
light elements produced during BBN is very sensitive to the Hubble parameter H and the 
temperature T. On the other hand, from the Friedmann Equation H and T are related to 
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the gravitational constant G by (assume radiation domination and zero curvature) 

H 2 ~ Gg*T 4 , (13.46) 



where g* is the number of relativistic species at BBN. Based on this relation, a constraint 



on G at the level of better than 10% have been inferred from BBN ||512 |. 

In addition, the influence of modified gravity on H can also affect the epoch of re- 
combination, resulting a detectable effect on the damping of the CMB power spectrum 
at high I. From current CMB observations, a constraint on G at ~ 10% level have been 



obtained [|513| . The constraint from CMB can help to break the high degeneracy of Hubble 
parameter H and temperature T in the BBN observations and improve the constrain to 
~ 3% ||513|| . In the future, it is expected that the Plank satellite can constrain G at the ~ 



1.5% level IBI3 



14. Dark energy projects 

Here we provide a brief overview of present and future projects to probe dark energy. 
Although some high energy physics experiments (such as LHC ||514|| ) might shed light 
on dark energy through discoveries about super symmetry or dark matter, here we only 
introduce the projects involving cosmological observations. According to the DETF report 



23], the dark energy projects can be classified into four stages: completed projects are 



Stage I; on-going projects, either taking data or soon to be taking data, are Stage II; 
intermediate-scale, near-future projects belong to Stage III; larger-scale, longer-term future 
projects belong to Stage IV. 



Moreover, to compare various dark energy projects, DETF J24| also proposed a quan- 
tity called figure of merit (FoM). Utilizing the famous Chevallier-Polarski-Linder (CPL) 
parameterization 

w(a) = wo + (1 - a)w a (14.1) 

a FOM is defined as the reciprocal of the area of an error ellipse in the wq — w a plane. A 
conventional normalization takes for the FoM the square root of the determinant of the 
2x2 Fisher matrix for wq and w a . Soon after, an extended version of FOM was proposed by 



Wang ||515|| . More advanced stages are expected to deliver tighter dark energy constraints, 
and give larger FoMs. For examples, Stage III experiments are expected to deliver a factor 
~ 3 improvement in the FoM compared to the combined Stage II results, while Stage IV 
experiments will improve the FoM by a factor of 10 compared to Stage II. 
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In this section, we will introduce some most representative projects of Stage II, Stage 
III and Stage IV. The corresponding dark energy projects are shown in Fig 7. The Stage 
I experiments that have already reported results, such as SNLS, ESSENCE, and WMAP, 
will not be discussed here. Notice that some dark energy projects have already launched 
since 2006, our classifications of dark energy projects are slightly different from that in the 
DETF report. For a comprehensive list of the dark energy experiments, see Ref. ITO . 
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Fig. 12: Dark energy projects introduced in this work. The black ones are stage 
II projects, the blue ones are stage III projects, and the red ones are stage IV 
projects. 



14-1- On-going projects 

In this subsection we introduce the on-going projects. A list of some most represen- 
tative projects of Stage II is given in Table VI. 

Table VI : On-Going Projects 
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Survey 


Location 


Description 


Probes 


SDSS BOSS 


Sacramento (USA) 


Optical, 2.5-m 


BAO 


LAMOST 


Xinglong (China) 


Optical, 4-m 


BAO 


Pan-STARRSl 


Hawaii (USA) 


Optical, 1.8-m 


SN, WL, CL 


SPT 


South Pole 


Submillimeter, 10-m 


CL 


HST 


Low Earth orbit 


Optical / near-infrared 


SN 


Planck 


Sun-Earth L2 orbit 


SZE 


CL 



The SDSS Baryon Oscillation Spectroscopic Survey (BOSS) ||516|| is one of the four 
surveys of SDSS-III, using a dedicated wide-field, 2.5-meter optical telescope at Apache 
Point Observatory in Sacramento Mountains. With a 5-year survey (2009 - 2014) of 1.5 
million luminous red galaxies (LRGs), BOSS will achieve the first measurement of the BAO 
absolute cosmic distance scale with 1% precision at redshifts z = 0.3 and z = 0.6. It will 
also measure the distribution of quasar absorption lines at z = 2.5, yielding a measurement 
of the angular diameter distance at that redshift to an accuracy of 1.5%. 



The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) |5T7 



is 



a National Major Scientific Project (NMSP) built by the Chinese Academy of Sciences 
(CAS). It is a wide-field, 4-meter ground-based optical telescope located at the Xinglong 
observing station of National Astronomical Observatories (NAO). The main construction 
of this instrument was finished in 2008. After its commissioning stage, LAMOST will 
study dark energy through the BAO technique. 

The Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) ||518 



is an international collaboration program led by University of Hawaii. Its first stage, Pan- 
STARRSl, is a 1.8-m wide-field telescope located at Haleakala in Hawaii. The mirror has 
a 3 degree field of view and be equipped with a CCD digital camera with 1.4 billion pixels. 
The regular observing of Pan-STARRSl has already started in 2009. As one of science 
goals, Pan-STARRSl study dark energy through the SN, the WL and the CL techniques. 

The South Pole Telescope (SPT) ||518 l is an international collaboration program oper- 
ating at the USA NSF South Pole research station. As the largest telescope ever deployed 
at the South Pole, SPT is a 10-m submillimeter-wave telescope with a 1000-element bolo- 
metric focal plane array. Constructed between 2006 and 2007, it is conducting a survey of 
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galaxy clusters over 4000 deg 2 using the Sunyaev - Zel'dovich effect (SZE). About 20,000 
clusters are expected to be discovered by the SPT in recent years. 

The Hubble Space Telescope (HST) | j520| is one of the most famous telescopes in the 
world. It was built by the USA National Aeronautics and Space Administration (NASA), 
with contributions from the European Space Agency (ESA). HST has two mirrors, the 
primary mirror diameter is 2.4 m, and the secondary mirror diameter is 0.3 m. Since 
launched in 1990, many HST observations have led to great breakthroughs in astrophysics, 
such as the discovery of the currently observed cosmic accelaration |l|J^]. So far, utilizing 
the HST, several SNIa datasets have been obtained, such as Gold ||322|J323|| , Union [ |331| , 



Constitution ||332|| and Union2 ||333|| . The present SN survey project of HST, HST Cluster 



Supernova Survey ||521| , is targeting supernovae in high redshift galaxy clusters at 0.9 < 



z < 1.5. It is expected that more than 100 high redshift SNIa will be found in recent years. 
The Planck |j522[ is a CMB satellite, which is created as the third Medium- Sized 



Mission of ESA's Horizon 2000 Scientific Programme. The telescope is an off-axis, aplanatic 
design with two elliptical reflectors and a 1.5 m diameter. It is designed to image the CMB 
anisotropies and polarization maps over the whole sky, with unprecedented sensitivity and 
angular resolution. While the WMAP data reach 200 billion samples after its nine-year 
mission, just a single year of observing by Planck will yield 300 billion samples. Besides 
pinning down the cosmological parameters of dark energy, Planck will also detect thousands 
of galaxy clusters using the SZE. It has launched in 2009, and is expected to yield definitive 
data by 2012. 

14-2. Intermediate- scale, near-future projects 

In this subsection we introduce the intermediate-scale, near-future projects. A list of 
some most representative projects of Stage III is given in Table VII. 

Table VII : Intermediate-Scale, Near-Future Projects 



Survey 


Location 


Description 


Probes 


HETDEX 


Davis Mountains (USA) 


Optical, 9.2-m 


BAO 


DES 


Cerro Pachon (Chile) 


Optical, 4-m 


BAO, SN, WL, CL 


Pan-STARRS4 


Hawaii (USA) 


Optical, 1.8-mx4 


SN, WL, CL 


ALPACA 


Cerro Pachon (Chile) 


Optical, 8-m 


SN 


BigBOSS 


Tucson and Cerro Pachon 


Optical, 4-m 


BAO 
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The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) ||523|| is an inter- 
national collaboration program led by University of Texas, Austin. It will use the 9.2-m 
Hobby-Eberly Telescope at McDonald Observatory to measure BAO over two areas, each 
100 deg 2 , using one million galaxies over the redshift range 1.8 < z < 3.7. It is expected 
that the survey will begin in 2012 and last for 3 years. 

The Dark Energy Survey (DES) [|524| is an international collaboration program led 
by Fermilab. It is a new 570 megapixel digital camera (DECam) mounted on the 4-m 
Blanco Telescope of the Cerro Tololo Inter-American Observatory (CTIO) in Chile. The 
DES plans to obtain photometric redshifts in four bands, and the planned survey area is 
5000 deg 2 . It will probe dark energy by using the BAO, the CL, the SN, and the WL 
technologies. The proponents will begin at 2011 and take 5 years to complete. 

The Pan-STARRS4 is the updated version of Pan- STARRS 1 |5l8^ . It is a large optical 



survey telescope to be sited on Mauna Kea in Hawaii. It consists of an array of four 1.8-m 
telescopes, each mirror will have a 3 degree field of view and be equipped with a CCD 
digital camera with 1.4 billion pixels. The dark energy science goals of Pan-STARRS4 
include SN, WL and CL surveys. The survey will continue for ten years. 

The Advanced Liquid-mirror Probe for Astrophysics, Cosmology and Asteroids (AL- 
PACA) ||525|| is a proposed wide-field telescope employing an 8-meter rotating liquid mirror. 



It brings together the technologies of liquid-mirrors, lightweight conventional mirrors, and 
advanced detectors, to make a novel telescope with uniquely-powerful capabilities. Scan- 
ning a long strip of sky every night at CTIO, ALPACA will discover ~ 50, 000 SNIa and 
~ 12, 000 type II SN each year to a redshift z ~ 0.8. 

The Big Baryon Oscillation Spectroscopic Survey (BigBOSS) ||526|| is an USA 



NSF/DOE collaboration program. As the updated version of SDSS BOSS, BigBOSS will 
probe dark energy through the BAO and the redshift distortions techniques. It will build 
a new 4000-fiber spectrograph covering a 3-degree diameter field. This instrument will be 
mounted on the 4-meter Mayall telescope at Kitt Peak National Observatory (KPNO) for 
a 6-year run, then will be moved to the 4-m CTIO Blanco Telescope for another 4-year 
run. After 10-year operation, BigBOSS will complete the survey of 50 million galaxies and 
1 million quasars from 0.2 < z < 3.5 over 24,000 deg 2 . The construction of the instrument 
will begin in 2011, and the operation of the survey will start in 2015. 
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lJf.,3. Larger- scale, longer-term future projects 

Table VIII : Larger-Scale, Longer- Term Future Projects 



Survey 


Location 


Description 


Probes 


LSST 


Cerro Pachon (Chile) 


Optical, 8.4-m 


BAO, SN, WL 


SKA 


Australia or South Africa 


Radio, km 2 


BAO, WL 


WFIRST 


Sun-Earth L2 orbit 


Infrared, 1.5-m 


BAO, SN, WL 


Euclid 


Sun-Earth L2 orbit 


Optical/NIR, 1.2-m 


BAO, WL 


1X0 


Sun-Earth L2 orbit 


X-ray 


CL 



In this subsection we introduce the larger-scale, longer-term future projects. A list of 
some most representative projects of Stage IV is given in Table VIII. 

A most ambitious ground-based dark energy survey project is the Large Synoptic 
Survey Telescope (LSST) ||527|| , which is an USA NSF/DOE collaboration program. LSST 
is an 8.4-meter ground-based optical telescope (with a 9.6 deg 2 field of view and a 3.2 
Gigapixel camera) to be sited in Cerro Pachon of Chile. Over a 10-year lifetime, it will 
obtain a database including 10 billion galaxies and a similar number of stars. As one of 
the most important scientific goals, LSST will study dark energy through a combination 
of the BAO, the SN and the WL techniques. Since its compelling scientific capacity and 
relatively low technical risk, LSST was selected as the top priority large-scale ground-based 



project for the next decade of astronomy in the Astro2010 report ||528|| . The appraised 
construction cost of LSST is 465 million U.S. dollars. The project is scheduled to have 
first light in 2016 and the beginning of survey operations in 2018. 

Another ambitious large-scale ground-based telescope is the Square Kilometer Array 
(SKA) ||529 |, which is an international collaboration program. SKA has a collecting area 
of order one square kilometer and a capable of operating at a wide frequency range (60 
MHz - 35 GHz) . It will be built in the southern hemisphere (either in Australia or in South 
Africa), and the specific site will be determined in 2012. SKA will probe dark energy by 
BAO and WL techniques via the measurements of the Hydrogen line (HI) 21-cm emission 
in normal galaxies at high redshift. The total Budget of SKA is 1.5 billion U.S. dollars. 
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Construction of the SKA is scheduled to begin in 2016 for initial observations by 2020 and 
full operation by 2024. 

A most exciting space-based dark energy project is the Wide Field Infrared Sur- 



vey Telescope (WFIRST) J53(J, which is an USA NASA/DOE collaboration program. 
WFIRST is a 1.5-meter wide- field near-infrared space telescope, orders of magnitude more 
sensitive than any previous project. The design of WFIRST is based on three separate 
inputs (JDEM-Omega, the Microlensing Planet Finder, and the The Near-Infrared Sky 



Surveyor) to Astro2010 [531]. It will enable a major step forward in dark energy under- 
standing, provide a statistical census of exoplanets and obtain an ancillary data set of great 
value to the astronomical community. The objective of the WFIRST dark energy survey is 
to determine the nature of dark energy by measuring the expansion history and the growth 
rate of large scale structure. A combination of the BAO, the SN and the WL techniques 
will be used to probe dark energy. Since its compelling scientific capacity and relatively low 
technical risk, WFIRST was selected as the top priority large-scale ground-based project 
for the next decade of astronomy in the Astro2010 report [ |528|| . The appraised construc- 
tion cost of WFIRST is 1.6 billion U.S. dollars. The project is scheduled to launch in 2020 
and has a 10-years lifetime. 

Another exciting space-based dark energy project is the Euclid |532(| , which is an 
ESA project. Euclid is 1.2 m Korsch telescope, with optical and near-infrared (NIR) 
observational branch. The primary goal of Euclid is to map the geometry of the dark 
universe, and it will search galaxies and clusters of galaxies out to z ~ 2, in a wide 
extragalactic survey covering 20,000 deg 2 , plus a deep survey covering an area of 40 deg 2 . 
The mission is optimised for two primary cosmological probes: BAO and WL. It will 
also make use of several secondary cosmological probes such as the ISW, CL and redshift 
space distortions to provide additional measurements of the cosmic geometry and structure 
growth. After 5 years' survey, Euclid will measure the DE EoS parameters wq and w a to 
a precision of 2% and 10%, respectively, and will measure the growth factor exponent 7 
with a precision of 2%. 

The last Stage IV project is the International X-ray Observatory (IXO) ||533|| , which 
is a partnership among the NASA, ESA, and the Japan Aerospace Exploration Agency 
( J AXA) . IXO is a powerful X-ray space telescope that features a single large X-ray mirror 
assembly and an extendible optical bench with a focal length of ~ 20 m. With more than 
an order of magnitude improvement in capabilities, IXO will study dark energy through 
the WL technique. Because of IXO's high scientific importance, it was selected as the 



fourth-priority large-scale ground-based project in the Astro2010 report ||528|| . The total 
Budget of IXO is 5.0 billion U.S. dollars. The project is scheduled to launch in 2021. 
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15. Observational constraints on specific theoretical models 

In this section we will briefly review some research works concerning the cosmological 
constraints on the specific theoretical models. These models can be classified into three 
classes: (i) models that modify the energy-momentum tensor on the r.h.s. of the Einstein 
equation, i.e. dark energy models. We will briefly review some numerical results of scalar 
field models, Chaplygin gas models, and holographic models, (ii) models that modify 
the l.h.s. of the Einstein equation, i.e. modified gravity models. We will briefly review 
some numerical results of the DGP braneworld model, the f(R) gravity, the Gauss-Bonnet 
gravity, the Brans-Dicke theory and the f(T) gravity, (iii) models that attempt to explain 
the apparent cosmic acceleration by assuring the inhomogeneities in the distribution of 
matter. We will briefly review some numerical works on the inhomogeneous Lemaitre- 
Tolman-Bondi model and backreaction model, Lastly, we will summarize and compare 
these theoretical models. 



15.1. Scalar field models 



As is well known, the most popular phenomenological models are models with rolling 
scalar fields (p ||231| , |366|| . These models are the direct generalization of the cosmological 



constant and have been well-studied both theoretically and numerically. In Sec. 11, we 
have introduced the theoretical studies on the scalar field models. In this subsection, we 
will discuss the scalar field models from the aspect of observations. 
• Reconstructing the scalar filed models from observations 

After the proposal of quintessence field as a candidate of dynamical dark energy 
|231| , |366| |, a lot of numerical studies have been carried out to distinguish this model from the 



cosmological constant (see [[108| , |536| , |537| , ^35| , |538| , |539| , p40|| for some early studies) . A widely 
used approach is directly reconstructing dark energy from the observational quantities like 
di,(z), H(z), and so on. In the quintessence model, the Einstein equations can be rewritten 
as p6| , |535| , ^Tl 



3H, 



o 



IP 

~h 2 



x dH 2 1 



6H 2 dx 



8nG fd(j)\ 
3H 2 [dx) 



2 dhiH QmX 



3H 2 x dx 



H 2 



x 



1 + z. 



:i5.D 



(15.2) 
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Fig. 13: The reconstructed potential V(z) and the kinetic energy term <fi 2 in 
units of the critical density p cr = 3MpH 2 , based on 54 SNIa data of the SCP 
including the low-z Calan Tololo sample [ 534 1 . Also plotted are the two forms of 
V{4>) for this V(z). The value of </>, known up to an additive constant, is plotted 
in units of the Planck mass. The solid line corresponds to the best-fit values of 
the parameters, the shaded area/dotted lines covers the range of 68%/90% errors, 
respectively. The hatched area represents the unphysical region <j) 2 < 0. From 
|535|. 



One can determine <p(z), and thus its inversion z(<p) by integrating Eq. (|15.2| ). The Hubble 
parameter H(z) together with its first derivative can be determined from observations, such 
as the measurements of the luminosity distance cLl(z) ||536| , p37| , |535| , |538| , |541 



H(z) 



d_ ( d L {z) 
dz \l + z 



15.3) 



Substituting H(z) and z(<p) into Eq. ( |15.1|) , one can reconstruct the potential V((f>). In 
2000 Saini et al. [|535| reconstructed V((f>) and wq(z) based on 54 SNIa given by Perlmutter 
et al. 0. Their results are shown in Fig. 13. 

In the same way one can also reconstruct other physical quantities from observations. 
For instance, in [ |535| , Saini et al. reconstructed the EOS w(z) using the relation 



w(z) 



{2/3){l + z)d1nH/dz- 1 
l-(H /H) 2 n m (l + z) 3 



(15.4) 
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They obtained the la constraint 



-Kw< -0.86, 



-l<w< -0.66, 



;i5.5) 



corresponding to the value of w at z = and z = 0.83, respectively. So the result allows 
the possible evolution of dark energy, while the cosmo logical constant with w = — 1 is still 
consistent with the data. 

In addition to the above method, a more widely used approach is to reconstruct the 



quantities by a fitting ansatz which relies on a small number of free parameters |[455|| [|542 
We will discuss this issue in the section 16. 

Broadly speaking, by making use of the scalar fields one can reconstruct a dark energy 
component with any property. So the issue of the observational tests of scalar field models 
is somewhat similar as the observational probe of the dynamical behavior of dark energy. 
One can see 



, [544| , |545| , |546| , [547| and references therein for more studies on the scalar 
field models from the aspect of observations. In addition, there are also some theories 
using vector fields to describe dark energy ||284| , |548| , |549| , |550 . 



• Quintessence, phantom or quintom 

Another interesting issue concerning the scalar filed models is the future evolution 
of dark energy and the fate of our universe. Besides the quintessence field satisfying 
w > — 1, another well-studied scalar field model is the phantom model satisfying w < 
— 1. In this model, the dark energy density will reach infinity and lead to the "big-rip" 
singularity. The WMAP7 measurements give w = —1.10 ± 0.14 ||349 ], the analysis of the 



Union2 SNIa dataset give w = — 1. 035^.0^7 ||333|| , and the analysis of the SDSS DR7 give 
w = —0.97 ± 0.10 ||385|1 . So the current observational data are still consistent with the 



cosmological constant, although the possibilities of quintessence and phantom all exist. 
Besides, the possibility of the quintom, where the EOS can cross w = — 1, can also provide 



a consistent fit to the observations. We refer to ||551|| for related numerical studies. 
• Interaction between dark sectors 

It is also worthwhile to consider the possible interaction between the scalar field and 



the matter component. In ||255|| , Amendola proposed the coupled quintessence (CQ) model, 
in which the scalar field (p and the dark matter fluid with each other through a source term 
in their respective covariant conservation equation 



;i5.6) 
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where T^,^ and T^,, represent the energy momentum tensors of </> and matter. In 
|255| , |552|| , it was proposed that the source term Q u assumes the form 



Q v = -K/3(0)T( m )V„< 



(15.7) 



where T( m ) is the trace of T^s, and (3 ((f)) (hereafter 0) is the coupling function that sets 
the strength of the interaction. 



SNLS + BAO + CMB 




SNLS + BAO + CMB 




Fig. 14: Observational constraints on the coupling between dark energy and dark 



matter using 71 SNIa from the fist year SNLS [324], the CMB shift parameter from 
the three- year WMAP observations [ |361| ], and the BAO peak found in the SDSS 
[ 374 1 . The constant coupling 5 = T/H is considered, which leads to the Friedmann 
equation p m + 3Hp m = 8Hp m . The left and right panels shows observational 
contours in the (wx,S) and (Qxo,S) planes, respectively (here "X" stands for 
the dark energy component). The best-fit model parameters correspond to 5 = 
-0.03, w x = -1-02 and Qxo=0.73. A stringent constraint -0.08 < 5 < 0.03 at 
95% CL is obtained. From pi. 



In ||255|| , Amendola also investigated the evolution of the perturbations as well as the 
observational signature of this model. Hereafter, a lot of studies were performed to con- 
strain the CQ model utilizing the various observational techniques ||554| , |555| , |556| , |553| , |557 l . 



For instance, in | 555 ], from first-year WMAP observations, Amendola et al. obtained a 
upper constraint ~ 0.1 for the coupling parameter (3. In [556] , Maccio et al. performed the 
A^-body simulations in CQ model and found that this model is consistent with the growth 
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of structure in the non-linear regime. In [ 553 ], Guo et al. considered the following kinds 
of interaction 



Pm + 3Hp m = 5Hp T 



:i5. 



From a combination of SNIa, WMAP and BAO data, they put a stringent constraints on 
the constant coupling mode 



-0.08 < 5 < 0.03, 



15.9) 



at the 2a CL (see Fig. 14). In [ 557 ] , Bean et al. investigated a variety of CQ models 
and found that a combination of SNIa, LSS and CMB data can constrain the strength of 
coupling between dark sectors to be less than 7% of the coupling to gravity. In all, the 
current observational data have already given tight constraints on the interaction between 
dark sectors. One can see ||558| , [559| and references therein for more studies on this topic. 
In addition, for the studies of dynamical system descriptions about the scalar dark energy, 
see pOl . 



15.2. Chaplygin gas models 

In ||288|| , Kamenschchik et al. explained dark energy as a kind of fluid called Chaplygin 
gas (CG), characterized by the EOS 

P=-~, (15.10) 
P 



where A is a constant. Later, Bilic et al. ||561|| and Bento et al. ||289|| proposed an extension 



of the original Chaplygin gas model, called generalized Chaplygin gas (GCG) model, with 
the EOS 

p=-4- ( i5 - n ) 

p a 

• Observational inspection of the Chaplygin gas models 
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Fig. 15: Left Panel: A constraint on the standard Chaplygin gas using SNIa, 



BAO and CMB. Clearly this model is a very poor fit to the data. From [562]. 
Right Panel: GCG as the unification of dark matter and dark energy will produce 
inconsistent oscillations in the mass power spectrum. The data points are the power 
spectrum of the 2df galaxy redshift survey, and the curves from top to bottom are 
GCG models with a = -1(T 4 , -1(T 5 , (ACDM), 1(T 5 and 1(T 4 , respectively. 
From |563|. 



The observational inspection of the CG and GCG models have been extensively inves- 
tigated using various observational methods ||562| , |563| , |564j , |565|| , including the observations 



of SNIa, BAO, CMB, WL, X-ray, GRB, X-ray, Fanaroff- Riley type lib radio galaxies [gggj. 
and so on. A common result of these studies is that the CG model has been ruled out 



by the observations. For example, in ||562|| , Davis et al. showed that the CG model is 
strongly disfavored by a combination of SNIa+BAO-fCMB data, since this model gives 
a Xmin ~ larger than that given by the ACDM model (see Fig. 15 for more details). 
Similar result was obtained in ||565|| , where Zhu investigated the GCG model from the 
measurements of X-ray, the SNIa and Fanaroff- Riley typellb radio galaxies, and got a 
constraint 

a = -0.09±°;1* (15.12) 

at the 95% CL. Therefore, the CG model, which corresponds to a=l, is ruled out at a 
99% CL. 

• The GCG model as a unification of dark energy and dark matter 

An attractive feature of the GCG model is that it can explain both dark energy and 
dark matter in terms of a single component and has been refered to as "unified dark 
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matter" (UDM) ||293| , 294 , 567 ]. From Eq. ( |15.11 ) one can obtain the energy density of 
GCG JjB7| 

i 

B 



Pit) 



A + 



fl3(l+a) 



l + a 



where B is an integration constant. Defining 

B 

P* 



then the Eq. ( |15.11|) becomes 

P(z) = 



A + B 



(A + B)^ 



Po 



l + a 



(15.13) 



[15.14) 



15.15) 



Obviously, the flat ACDM scenario with matter ratio is recovered with a = 0. 

However, a problem with the UDM scenario is that it will produce oscillations or 
exponential blowup of the matter power spectrum, which is inconsistent with observation. 
In |[563|| , Sandvik et al. found that the density fluctuation 5k with wave vector k evolves 
as 



SZ + [2 + (-3(2w-c 2 s )}S' k 



~(1 - 6ct + 8w- 3w 2 ) - ( ^ 
1 1 aH 



5k 



[15.16) 



where the quantity £, the EOS w and the squared sound speed c 2 take the form [|563 
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—aw, 
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1 — 

and prime denotes partial differentiation with respect to In a. For Chaplygin gas, a > 
will result in c 2 s < 0, which leads to a non-zero Jeans length Aj = ^n\c 2 s \/Gp. In this 
regime the fluctuations will oscillate instead of grow polynomially. On the other hand, 
when a < 0, perturbations will grow exponentially, which is also ruled out. 

In ||563| 1 , Sandvik et al. further confirmed the above arguments by numerical solving 
the equations (see Fig. 15). By performing a x 2 fit of the theoretically predicted power 
spectrum against that observed by the 2dFGRS ||382|| , they found that this inconsistency 
excludes most of the previously allowed parameter space of a, leaving essentially only 
the standard ACDM limit. This topic was also studied in |[568| , |569|| . In [[569]], based on 
an analysis including the SNIa, the baryonic matter power spectrum, the CMB and the 
perturbation growth factor, Park et al. found that the allowed region for a is 

-5 x 10" 5 < a < 10" 4 . (15.19) 

The allowed parameter space is extremely close to the ACDM model, so the possibility of 
the unification of dark matter and dark energy in the GCG scenario has been ruled out 
by the current cosmological observations. 
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15.3. Holographic dark energy models 

In the following, we will introduce some numerical works about the holographic dark 
energy (HDE) model, which arises from the holographic principle. 
• The HDE model with the future event horizon as the cutoff 

The HDE has the form of energy density 



Pde 



(15.20) 



where c is a number introduced in ||162|| . In 2004, Li ||162|| proposed to take L = 1/Rh, 
where Rh is the future event horizon defined as 



Rt. 



a 



da 
~Ha~ 2 



That yields an EOS 



w 



1 2 yOde 

3 ~ 3 c 



(15.21) 



(15.22) 



which satisfies w < —1/3 and can accelerate the cosmic expansion. 



In ||570|| , Huang and Gong first performed a numerical study on the HDE model. 
Making use of the 157 gold SNIa data, they obtained fi m = 0.25iS;SJ, w = -0.91 ± 0.01 
at la CL for the c = 1 case. In [ p71| , Zhang and Wu constrained the HDE model by 
performing a joint analysis of SNLa+CMB+LSS. They found that the best fit results are 
c = 0.81, = 0.28, and h = 0.65, which implies that the HDE behaves as a quintom- 
type dark energy. In |572fl , Yi and Zhang tested the HDE model by using the H(z) data 
and the LT data. They also verified that the HDE behaves like a quintom-type at la 
CL. In addition, a lot of numerical studies were performed to test and constrain the HDE 
model [|573|1 . These works showed that the HDE model can provide a good fit to the data. 
For example, in [|574| , by using the combined Constitution+BAO+CMB data, Li et al. 
obtained the following Xmin s f° r the ACDM and HDE models 



xlcDM = 467.775, 



Xhde 



465.912. 



(15.23) 



So the HDE model is consistent with the current observations. Similar results have been 



obtained in e.g. ||575| , |576| , |577| , |578| . Therefore, from the perspective of current observations, 
HDE is a competitive model. 
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Fig. 16: Left Panel: Constraints on the HDE model in Q m -c plane from ob- 



servational data including Constitution SNIa sample | 332 |, the BAO data from 



SDSS p7i] , and the CMB data from WMAP5 measurements pnfl . The data favor 
c < 1. For the constraint on c, the analysis gives 0.8181q'^54 at the 68.3% CL, and 
0.818lg:o97 at the 95.4% CL. From f74}. Right Panel: The contour maps of a vs. 
c for interacting HDE (IHDE) with 68%, 95.5% and 99.7% CL, obtained from a 
joint analysis including the Golden06+BAO+X-ray+GRB+CMB data. The black 
dot-dashed curve denotes w e ff = —1 when z — > —1 with Qde = 0.73, and the 
region below (over) it means w e ff will (not) cross -1 during infinite time. From 
[5761 



The parameter c plays an essential role in determining the evolution of the HDE. If 
c = 1, the dark energy EOS will be asymptotic to that of a cosmological constant and the 
universe will enter the de Sitter phase in the future; if c > 1, the EOS is always greater 
than -1, and HDE behaves as quintessence dark energy; if c < 1, initially the EOS of HDE 
is greater than -1, then it will decrease and eventually cross the w = —1 line, leading 
to a phantom universe with big rip as its ultimate fate. The numerical studies on HDE 



showed that the cosmological observations favors c < 1. For example, Ref. [|574| gave 
c = 0.818^^097 a t the 68% CL (see the left panel of Fig. 16), while some later works 



423| , |578|| with the improved data tighten the constraint to c < 0.9 at the 95% CL. 



Interactions between HDE and dark matter 



In ||167| , |579|| , Wang et al. first studied the interaction between dark matter and the 



HDE. The introduction of interactions not only alleviates the cosmic coincidence prob- 
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lem but also avoids the future big-rip singularity ||47ip75J57(: |, |581 , |582| , ^83| , |584|| . Phe- 
nomenologically, the interaction can be introduced by 



Pm + 3Hp m = Q, p de + 3H(p de + p de 



-Q. 



(15.24) 



where Q is the interaction term. Current observations have put tight constraints on the 
interactions. In [[5761 , with the assumption of a flat universe, Ma considered the following 
Q and obtained the corresponding constraint 



Q = 3aH Pde , a = -0.006±g;gli, 



;i5.25) 



from a joint analysis (see the right panel of Fig. 16). In ||583|1 , Feng et al. considered 
another class of Q and obtained the corresponding constraint 



q = 3bH( Pm + Pde ), b = -o.oo3lg;8}l, 



(15.26) 



from a combination of SNIa+BAO+CMB+Lookbacktime data. Later, in ||584jl , Li et al. 
revisted these two models and obtained 



a 



;-e.i x io- e )±g : gg, 



b = (-1.6 x io- 4 )l° ; °o 8 9- 



;i5.27) 



from a joint analysis of Constitution+BAO+CMB. They also found that there exists sig- 
nificant degeneracy between the phenomenological interaction and the spatial curvature in 
the HDE model. 
• The ADE and RDE models 

In addition to the HDE model with future event horizon as the cutoff, the Agegraphic 
dark energy (ADE) model [|178| , |179| , |585|1 and the Holographic Ricci dark energy (RDE) 
model ||182|| are also motivated by the holographic principle (the ADE model can also be 
obtained from the Karolyhazy relation; see [|178| for details). In these two models, the IR 
cutoff length scale is given by the conformal time r\ and the average radius of the Ricci 
scalar curvature |7£| -1//2 , respectively. There have been some numerical studies on these 
two models [|574| , |586| , |587|,|588| , |589| , |590|1 . In general, these studies showed that the ADE and 
RDE models are not favored by current observations. For example, in |574j| , Li et al. 
obtained 

Xade = 481.694, X rde = 483.130. (15.28) 

The vi- s of the ADE and RDE model are much larger than that of the ACDM and HDE 
model listed in Eq. ( |15.23| ), showing that these two models are not favored by observations. 
The results have been further confirmed in some later works [[423| , |578|1 . 



17 coincidence problem may also be solved by the anthropic constraints [ 580 1 
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15.4- Dvali-Gabadadze-Porrati model 



The Dvali-Gabadadze-Porrati (DGP) braneworld model is a theory where gravity 
is altered at immense distances by the slow leakage of gravity off our three-dimensional 
universe 



1 1-4|| - In this model, the Firedmann equation is modified as 

87rG 



r r 



-Pr 



(15.29) 



where r c = {Hq{1 — O m )) _1 is the length scale beyond which gravity leaks out into the 
bulk. At early times, Hr c 3> 1, the Firemann equation of general relativity is recovered. 
In the future, H — > l/r c , the expansion is asympototically de Sitter. 

There is also a generalized phenomenological DGP model characterized by the Fried- 



mann equation 591 



1 c 



8tvG 



(15.30) 



where r c = H — fl r 



\a-2 



This model interpolates between the pure ACDM model 



and the DGP model with an additional parameter a. a = 1 corresponds to the DGP 
model and a = corresponds to the ACDM model. 
• The DGP model is disfavored by the observations 

Although DGP is an attractive model allowing a self acceleration, many research 
works show that it is disfavored by observations [ 592 , 593 , |594 , |595| , [)96| , 597 1 . For examples, 
592]] , from a joint analysis of SNIa+BAO+CMB, Rubin et al. found that the DGP 



m 



was disfavored by the data, with a Ax 2 =15 compared with the ACDM model (see the left 
panel of Fig. 17). In [ |594| , from a combination of SNIa and BAO measurments, Guo et al. 
provide the constraints to the model parameters 



n 



+0.018 r, 
u - z ' -0.017' ^ L r c 



0.216™, Q k = - .350^- 080 



083' 



(15.31) 



at 99.73% CL. This result is in contradiction to the WMAP results indicating a flat uni- 
verse. Moreover, the constraints to the generalized DGP model gives a small a, indicating 



that the DGP is incompatible with the observations. In a recent work ||598|| , Xia performed 
a joint analysis including SNIa, BAO, CMB, GRB and the linear growth factor of matter 
perturbations, and found a constraint 



a = 0.254 ±0.153, 



(15.32) 
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Fig. 17: Le/t Panel: Combined SNIa+BAO+CMB observational constraint shows 
that DGP does not achieve an acceptable fit. The areas of intersection of any 
pair are distinct from other pairs. This is a strong signal that the DGP model is 



incompatible with the observations. From [592]. Right Panel: Predictions for the 
power spectra of the CMB temperature anisotropics Cj T of the best-fit DGP model 
(solid), a quintessence model with the same expansion history as DGP (short- 
dashed), and the ACDM model (dashed, conincident with the quintessence model 
at low I). Obtained by fitting to SNLS+ WMAP5+HST assuming a flat universe. 
Bands represent the 68% and 95% cosmic variance regions for the DGP model. 
Points represent WMAP5 measurements. It is clear that the best-fit DGP model 
over predicts the law-l modes anisotropy. From [593]. 



at the 68% CL, manifesting that this model tends to collapse to the cosmological constant 

when confronted with current observations. 

• Testing the DGP model from the growth of structure 

As a modified gravity scenario, the growth of structure in the DGP gravity differs 



from that in the ACDM scenario. This can be used to test the DGP model ||477| , |599| | ■ The 
perturbation theory in the DGP model has been studied [|600| , |601| , |602| , |603|| . These studies 
showed that the DGP gravity is disfavored by the observational data. For examples, in 



601j| , Song et al. showed that the constraints from SNIa+CMB+i^o exclude the simplest 



flat DGP model at about 3cr. Even including spatial curvature, best-fit open DGP model 
is a marginally poorer fit to the data than the flat ACDM model. In ||593 |, Fang et al. 
showed that the DGP model is excluded at 4.9a and 5.8a levels with and without curvature 
respectively (see the right panel of Fig. 17). The corresponding Xmin s f° r ^ ne DGP and 
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ACDM model are 

Xacdm = 2777.8, X dgp = 2805.6. (15.33) 

The result is mainly due to the earlier beginning of the acceleration and the additional 
suppression of growth in the DGP scenario. In [ |604|| , by performing cosmological iV-body 
simulations of the DGP model, Schmidt found that, independently of CMB constraints, 
the self-accelerating DGP model is strongly constrained by WL and cluster abundance 
measurements. Compared with the ACDM model, the abundance of halos above 10 14 Mq 
is suppressed by more than a half in the DGP model. In all, when confronted with 
experiments, a lot of problems will emerge in the DGP model, indicating that this theory 
is strongly disfavored by the observations 

15.5. f(R) models 

f(R) gravity is a simplest modification to the general relativity with the replacement 
R — > f(R). In this section, we will focus on the observational tests of this theory. One can 
see [ P^2] , p8] , |5Tl5] , |6U6| , |6Tr7|J6U8[j and references therein for more details. 
• f(R) gravity and its viable conditions 

At the beginning, it was proposed to explain the cosmic acceleration using the model 
with f(R) = R — a/R n (a > 0,n > 0) ||609| , |610| , [611|| . However, later on, a lot of problems 



emerged in this model. In [|612|| , it was shown that this model will lead to the matter 
instability. In ||613|| , it was also found that this model is unable to satisfy local gravity 
constraints and pass the solar system tests of gravity. Much attentions have been paid 
to the analysis of the viable conditions of f(R) models, and a lot of valuable results have 



been obtained ||614| , |615| , |616| , |617| , |618| , |619|| . For examples, to have a stable perturbation, the 



condition /^rr = d 2 f/dR 2 < is required; to have a stable late-time de Sitter point, 
the condition < Rj \rr/ f \r < 1 is also necessary. In summary, the conditions for the 



viability of f(R) dark energy models include [|12] , |6 1 8| , |619|] : 



(1) f >R > and f,R R > for R < R ; 

(2) f(R) ^R-2A for i?<i? ; (15.34) 

(3) < Rf,RR/f,R < 1 at the de Sitter point satisfying Rf^R = 2f. 

To be acceptable, an f(R) model must satisfy the following conditions. Some viable 
models satisfying all these requirements have been proposed, and one can refer to e.g. 
12l, |606| , |6T8|j619| , |620| , |621| , |622l for more details about the viable conditions of f(R) gravity. 



Cosmological tests of the f(R) gravity 



132 



§ 10 4 



- 


: 








i \voo\. 




0.1 \w 




o.oi ysv\ 




0.001 \\ 




(ACDM) V 



all data 

with CA 

with-glSW 
(B -> 10 L! So) 



0.001 



(Mil 



0.1 




k (/j/Mpc) 

Fig. 18: Left panel: Linear matter power spectrum for several values of -Bo in the 
ACDM expansion history in the model [ 623 , 624 . \fB is of order the Compton 
wavelength of the new gravitational degree of freedom, measured in the unit of 
Hubble length. For length scale larger than \fB/H, the modified gravity effect is 
suppressed because of the mass of the scalar gravitaton. For length scales smaller 
than t/B/H, gravity is modified. The modification of gravity leads to an enhance- 
ment in the growth of perturbations and a corresponding suppression for the decay 



of the gravitational potential. From [524]. Right panel: Contours of 2D marginal- 
ized 68%, 95% and 99% confidence boundaries using observational data from the 
measurements of CMB, BAO, Hubble constant, and so on. For the constraints 
from galaxy-ISW (gISW) cross correlation, the Compton wavelength is rescaled as 
5 ->■ 1(T 2 -B . From 



We have listed some general conditions required for an f(R) model to be valid. An 
f(R) model satisfying these basic requirements, furthermore, should be confronted with 
the cosmological observations. A natural method is to test the f(R) theories from the 
observations about the growth of structure, which depends on the theory of gravity. This 
issue attracted a lot of interests [|616| , |623| , |624| , |626|| . Some observational signatures of the 
f(R) models have been presented. For example, in [|623| , |624|1 , Song et al. studied an f(R) 
model parameterized by "Compton wavelength parameter" B, which is proportional to 
the second order derivative f >R R 

f )RR HdR/d\na , 1 , Q ,, 
B = ~f JWTT\ • 15.35 

The B < branch violates the third condition of Eq. (|15.34|) . For the stable B > branch, 
it was found that this model will predict a lower large-angle CMB anisotropy by reducing 
the ISW effect, qualitatively change the correlations between the CMB and galaxy surveys, 
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and alter the shape of the linear matter power spectrum (see the left panel of Fig. 18 for 
details) . 

To test of the modified gravity theory, it is worthwhile to include various observational 
techniques probing gravity in different scales ||505|| . For example, in [|627 |, Schmidt et al. 
conduct the first, simulation calibrated, cluster abundance constraints on a two-parameter 
modified action model 



f(R) =R-2a- 



R 



(15.36) 



R + V 2 ' 

They found that the local cluster abundance, when combined with the data from CMB, 
SNIa, Hq and BAO, can lead to a very tight constraint to the model, improving the 
previous constraints by 3-4 orders of magnitude. The reason is that, the inclusion of 
cluster abundance data improves the bounds on the range of force modification from the 



several Gpc scale to tens of Mpc scale. In ||625|| , Lombriser et al. revisited the model 
studied in [|623| , |624jl and reported a strong constraint to the current value of the Compton 
wavelength 

B Q < 1.1 x 1(T 3 , (15.37) 

at the 2a CL (see the right panel of Fig. 18), mainly due to the inclusion of data from 
cluster abundance. 

Some other observational techniques have also been used to study the f(R) models, 
auch as the cosmic shear experiments [|628| , |629|| , the "21 cm intensity mapping" (detection 
of LSS in three dimensions without the detection of individual galaxies) [ |377| , |63U|| , the 
variation of the fine structure "constant" ||631| , |632|| , and so on. Nowadays, the observational 
tests of f(R) gravity has drawn increasingly attention. For more details on the test of 
modify gravity theories from observations, see [p|, |505| , |633l and the references therein. 
• The curvature singularity problem in f(R) gravity 



In ||634| , |635|| , It was proposed that there is curvature singularity problem in f(R) 
theories, due to the dynamics of the effective scalar degree of freedom in the strong gravity 
regime. This problem leads to the contradiction with the existence of the relativistic stars 
(like neutron stars) ||636|| . Furthermore, in ||637|| , it was shown that this problem can be 
cured via the addition of R 2 term. One can see Ref. [ b38 for some studies on this topic. 
• f(R) theories with Palatini approach 

The observational constraints on f(R) theories within Palatini approach was firstly 



performed by Amarzguioui [639] et a/., where they investigated the parameterization with 
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the form of f(R) = R + aR p using the cosmological measurements SNIa+BAO+CMB. 
The best fit values is found to be j3 = 0.09 with the allowed values rage 



|0| < 0.2, 



(15.38) 



at la CL, and the previously commonly considered 1/R model corresponding to (3 = — 1 
is ruled out by the constraint. In |640[ j, using the same form of parameterization, Koivisto 
calculate the power spectrum in the Palatini formulation of f(R) gravity. By comparing 



the results to the SDSS data 641 , it was found that the observational constraints reduce 



the allowed parameter space to 



<~ 10" 



;i5.39) 



which is a tiny region around the ACDM. Besides, the Palatini f(R) gravity is also faced 
with some problems associated with non-dynamical nature of the scalar-field degree of 
freedom. One can refer to ||642| , |643| , |644|| for more studies about the Palatini f(R) gravity. 



15.6. Other modified gravity models 

We have introduced the observational aspects of the DGP scenario and the f(R) 
models. In this section, we will discuss some other modified gravity theories, including the 
Gauss-Bonnet gravity, the Brans-Dicke gravity and the f(T) gravity. 
• Gauss-Bonnet gravity 

When compared with observations, many problems arise in the Gauss-Bonnet model 
|645|| . In ||646| , |647|| , from the data analysis including the constraints from BBN, LSS, BAO 



and solar system data, Koivisto and Mota found that this model is strongly disfavored 
by the observations. In [[348 and [|649 , Li et al. and De Felice et al. investigated the 
Gauss-Bonnet gravity, and found that the growth of perturbations gets stronger on smaller 
scales. This is incompatible with the observed galaxy spectrum, unless the deviation from 
the Einstein gravity is very small. Thus, the Gauss-Bonnet models have been effectively 
ruled out. 

• Brans-Dicke theory 

Another well-known modified gravity model is the Brans-Dicke theory [ [117 . For the 
Brans-Dicke theory, the PPN parameter 7 (see Eq. ( |13.45|) ) takes the form 



7 = (l + w)/(2 + £j), 
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where uj is the constant in Eq. ( 7.79 ). From solar system and binary pulsar observations, 
the parameter uj has been tightly constrained to [ |5U4 



uj > 4 x 10 4 



(15.41) 



and the cosmological effects of the scalar field are rendered insignificant. Later on, 
the Brans-Dicke theory is generalized to the scalar tensor theory [ p.18 . One can see 
]T^ , pU5| , |6Tl6] , |65Cl| , |651| and the references therein for more details on the test of this model. 
• f(T) theory 

There have been some numerical studies on the recently proposed f(T) gravity 
||126| , |127|| as an explanation of the cosmic acceleration. In ||652|| , Wu and Yu examined 
the following models from the Union2 SNIa dataset together with the BAO and CMB 
data 



f(T) = a(-T) n , f(T) = -aT(l - e pTo/T ) 



15.42) 



and obtained the constraint n = O.O^g;^, p = — 0.02^;^, at the 95% CL. They also 
compared the two models with the ACDM by using the Xmm/^ / (dof: degree of freedom) 
criterion. The results showed that ACDM is mildly favored by the data. Later, the power- 
law model was revisited by Bengochea [|653|| , with the inclusion of GRB and H$ data into 
consideration. In [ 654j| , Wei, Ma, and Qi tried to constrain /(T) theory by using the 
varying fine structure "constant"; it is found that the observational Aa/a data make f(T) 



theory almost indistinguishable from the ACDM model. In addition, in [|655|| , they also 
constrained f(T) theory by using the varying gravitational "constant"; it is found that the 
allowed model parameter n has been significantly shrunk to a very narrow range around 
0. 

In ||656|| , Bamba el al. studied the cosmological evolutions of the EoS for dark energy 
Wde in the exponential and logarithmic as well as their combination f(T) theories. They 
found that the crossing of the phantom divide line of Wd e = — 1 can be realized in the 
combined f(T) theory even though it cannot be in the exponential or logarithmic f(T) 
theory. Moreover, the crossing is from Wd e > — 1 to Wd e < — lj which is favored by the 
recent observational data. 

The perturbations in f(T) gravity has been studied in ||657| , |658| , |659|| . In [ p59| , Zheng 



and Huang derived the evolution equation of growth factor for matter over-dense pertur- 
bation in f(T) gravity. In addition, a problem in f(T) gravity was pointed out by Sotiriou 
et al. in [|660|| , where they showed that the Lorentz symmetry can not be restored in / (T) 
theories due to sensible dynamics. 



136 



15.7. Inhomogeneous LTB and backreaction models 

The inhomogeneous models have gathered significant interest in recent years as a 
scenario to explain the cosmological observations without invoking dark energy . In this 
section, we will briefly discuss the observational signature of the LTB model and the 
backreaction model. 
• LTB models 

The LTB mo dels [[HI are the most commonly considered inhomogeneous scenario 
to explain the cosmic acceleration without introducing the dark energy component. For 
LTB models the structure of our universe is described by the inhomogeneous isotropic 
Lemaitre-Tolman-Bondi (LTB) metric [[198yi99| , |200f 



ds 2 = -dt 2 + Rl{ri !)\ dr 2 + R 2 (r,t)(d6 2 + sin 2 ## 2 ), (15.43) 
l + p(r) 

where the prime denotes partial differentiation with respect to r. (3(r) is a function of 
r. Notice that the FRW metric is recovered by requiring R = a(t)r and p = -kr 2 . The 
Hubble parameters at the transverse and radial direction are expressed as 

H± = §, #,, = 5, (15-44) 

and the apparent cosmic acceleration can be explained by choosing suitable form of R(r, t), 
without the introduction of dark energy ||662|| . In the simplest class of such models we live 
close to the center of a huge, spherically symmetric Gpc scale void. Due to the spatial gra- 
dients in the metric, our local region has a larger Hubble parameter than the outer region 
663| ] . In some more complicated scheme, it have been proposed to reconstruct the cosmo- 



logical constant in an inhomogeneous universe [ |664| | . The idea is that, since cosmological 
observations are limited on the light cone, it is possible to reconstruct an inhomogeneous 
cosmological model (indistinguishable from the homogeneous ACDM model) to explain the 
cosmic acceleration without a cosmological constant. So far there have been some studies 
on the observational tests of the LTB model. In the following, we will briefly review some 
related works. 

It has been shown that the void model is able to provide a good fit to the SNIa data. 



For example, in ||665fl , Sollerman et al. tested two kinds of LTB models and compared them 
with the ACDM model. The models they considered has the following matter distribution 

On(r-) = O out + (O in - O out )e- (r / ro)2 , (15.45) 
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O m (r) = O out + (O in - O 



out j 



1 + e -r /Ar 



;i5.46) 



1 -|- e (r-r )/Ar 

where 0; n is the matter density at the center of the void, O ou t is the asymptotic value of 
the matter density outside the void, and ro is the size the underdensity. The second model 
|666|| has a much sharper transition of matter density than the first one, with the extra 



parameter A r characterizes the transition width. Using the first-year SDSS-II SNIa data 
||330|| together with the BAO and CMB measurements, they found that x 2 values for the 
LTB fits are comparable to that of the ACDM model 



XACDM — 233.2, XlTB modell 



235.5, Xltb modei2 — 237.6. 



[15.47) 



while the extra parameters in the LTB models make the models fare poorly in the infor- 
mation criteria tests. 
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Fig. 19: Left panel: An off-center galaxy cluster in a void will observe a dipole 
in the CMB. Since the expansion rate inside the void is higher, photons arriving 
through the void (from the right in the figure) will have a larger redshift (Azi n ) 
than photons did not pass through the void (left, with redshift Az out ). From J667| . 
Right panel: The changes in the \ 2 values as a function of the observers position. 
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The SDSS-II data [ 330 combined with the CMB dipole requirement is used. The 
X 2 value quickly increases as the observer is displaced away from the center. See 
3681 for details about the meaning of the diamonds, stars, lines and arrows. From 



mi. 



However, a lot of problems arise in the LTB model when confronted with some other 



cosmo logical observations. In [669], Alnes et al. pointed out that a problem of the void 
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model is that it requires us to live precisely near the center of the void. If there is a 
deviation between our position and the center of the void, the observed CMB dipole would 
become much larger than that allowed by observations (see the left panel of Fig. 19 for a 
brief description). Currently, the maximum distance to the center have been constrained 
to be very small p0| , |671| , |672l 

r obs <~ 20 Mpc, (15.48) 



which leads to a fine-tuning problem. For example, in [ |668| Blomqvist and Mortsell tested 
the LTB described by Eq. (|15.45| ) by using the SNIa and CMB dipole data. They found 
that the position of the observer has been confined to within about one percent of the void 
scale radius (see the right panel of Fig. 19). 

Besides, even if we happen to live very close to the center of the void, there should 
be some off-centered galaxy clusters where a large CMB dipole can be observed in their 
reference frame. For us, the relative motion between the CMB frame and the matter frame 



manifests itself observationally as a kinematic Sunyaev-Zeldovich effect. In [|667|1 , Bellido 
et al. demonstrated that the limited observations of only 9 clusters with large error bars 
already rule out LTB models with void sizes greater than ~ 1.5 Gpc and a significant 
underdensity. 

Another problem with the void models is that when fitted to the data, they slow local 
expansion values ||671| , |672| , |673| , |S74| , |675 1 



ho ~ 0.45-0.6, 



(15.49) 



which seems to be in contradiction with the measurements of the local Hubble constant. 
Recently, in |[449|| , Riess et al. obtained 



H = 73.8 ± 2.4 km/s/Mpc, 



;i5.50) 



corresponding to a 3.3 % uncertainty. They found that the void models with ho ~ 0.6 is 
ruled out by the measurements in more than 5 a. 

There have been some other useful methods proposed to distinguish the void models 
between other inhomogeneous dark energy models. In ||676|| , Uzan et al. presented the 
redshift drift z in a general spherically symmetric spacetime, and demonstrated that its 
observation would allow the test of Copernican principle. In addition, Yoo [|677|| and 
Quartin [|678|| showed that the dz/dt in void modes is always negative, which is greatly 
different from other dark energy models. Another interesting idea is, the ionized universe 
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severs as a mirror to reflect CMB photons in other regions of the universe to us, and thus 
can tell us deviation from the Coper nican principle ||679|,680 l. Utilizing this method, some 
models with largest voids have be excluded [|S80| . There are also some other methods of 
testing the LTB model, such as the scalar perturbations [|674|1 , the slope of low z SNIa 



distance moduli 6811, the constant curvature condition 



the cosmic neutrino background [K84|, the cosmic age test 



81, the small scale CMB [§83 
and so on. 



In all, in recent years the research of the LTB model has drawn a lot of interests. 
Although this model can provide a good explanation of the SNIa data without introducing 
the mysterious dark energy component, when confronted with the CMB, Hq and some 
other cosmological observations, a lot of problems emerged. For more research works 
about the LTB model see [|686| , |687| , |688|| . For a recent review, see ||689| . 
• Backreaction model 

Like the LTB model, the backreaction model JTgg[I^gggggI].gggJ is another model 
in which the cosmic expansion is due to the effect of the inhomogeneities of the universe. 
There have been some works concerning the possible observational signature of this model 
693| , |694| , |695|| . For example, in ||693|| , Li et al. showed a non-trivial scale dependence of 



the Hubble rate in this model. In ||695|| , Larena et al. proposed to use a template metric 
to deal with observations in backreaction context and found that averaged inhomogeneous 
models can reproduce the observations of SNIa data and the position of the CMB peaks. 
To provide evidence for the backreaction mechanism, further studies are needed. One can 
see Refs. ||686|,|696|| for mode studies concerning the backreaction model. 



15.8. Comparison of dark energy models 

Facing so many dark energy candidates, it is very important to decide which one is 
more favored by the observational data. So far there have been many works on the com- 
parison of various dark energy models P23|J562| , |578| , |592|J697| , |698| , |699| , |70C1|| . We will briefly 
review the topic of model comparison in this subsection. 
• Model selection and the information criteria 

The x 2 statistics alone cannot provide an effective way to make a comparison between 
competing dark energy models. To do this, we should take into account the relative 
complexity of the models. To give a blatant example, a lOth-order polynomial will always 
give an equal or better fit than a straight line to any data set, but this does not mean 
that any of the extra eight coefficients have any significance. It juts means that a model 
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with more parameters will generally give an improved fit (always, if the simpler model is 
a subclass of the more complex one) ||562 |. 



To enforce a model comparison, a general way is to employ the information criteria 
(IC) to assess different models [ [701| , |702| , |703|| . These statistics favor models that give a good 



fit with fewer parameters. The most frequently used IC including the Bayesian information 
criterion (BIC) [[701 1 and the Akaike information criterion (AIC) ||702|| . They are defined as 



BIC = -21n£ max + klnN (15.51) 

AIC = -2 In £ max + 2k (15.52) 

where £ max is the maximum likelihood (under Gaussian assumption x min = —2 ln£ max ), k 
is the number of parameters, and N is the number of data points used in the fit. According 
to these criteria, models that give a good fit with fewer parameters will be more favored. 
So these criteria embody the principle of Occam's razor, "entities must not be multiplied 
beyond necessity". Generally, a ABIC of more than 2 (or 6) is considered positive (or 



strong) evidence against a model ||704|| . It should be noted that the IC alone can at most 
say that a more complex model is not necessary to explain to current data, since a poor 
information criterion result might arise from the fact that the current data are too limited 
to constrain the extra parameters in this complex model, and it might become preferred 
with improved data. 

In addition to the AIC and BIC, a more sophisticated method for model selection is 
the so-called Bayesian evidence (BE), which considers the increase of the allowed volume 
in the data space due to the addition of extra parameters rather than simply counting 
parameters. So it requires an integral of the likelihood over the whole model parameter 
space 

BE = J £(d\0,M)p(9\M)dd (15.53) 

In practice, the integral in Eq. ( |15.53| ) can be calculated using the nested sampling method 
||705|| . BE has already been applied in cosmology ||706|J707j| . Since in most cases the 



simpler AIC and BIC methods are sufficient to employ a model comparison, they are more 

commonly used compared with BE. 

• Comparison of different dark energy models 
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Fig. 20: Graphical representation of the results of the IC values of some popular 
dark energy models. AAIC and ABIC are represented by the light and dark grey 
bars, respectively. From left to right, the models are given in order of increasing 
AAIC. The crosses mark the number of free parameters in each model (right-hand 
ordinate). The "unsupported" or "strongly unsupported" lines stand for a ABIC 
of 2 and 6. Clearly the flat ACDM is the most preferred model. From [562]. 



In [[562], based on the observational data of SNIa, BAO and CMB, Davis et al. scruti- 
nized and compared a number of dark energy models by using AIC and BIC. They found 
that the ACDM model almost achieves the best fit of all the models despite its economy 
of parameters. A series of models, including the XCDM model with constant EOS of dark 
energy w, the Cardassian expansion model and the CPL parametrization, can also provide 
comparably good fits but have more free parameters. The DGP model and the standard 
CG model with a = 1 are clearly disfavored. Based on their AIC and BIC values, one can 
determine the "rank" of these dark energy models (see Fig. 20) . From the figure, it is clear 
that the ACDM model is best favored, while the DGP model and the standard CG model 
is strongly disfavored. 

By comparing the values of IC, similar results are obtained by Szydlowski et al. [ B98 |, 
A. Kurek et al. ||699|| , Li et al. ||578|| and Wei et al. [123 ] . In the above works, the authors 
all found that the one-parameter flat ACDM performs best in the set of models considered 
in the context. These results further solidified the status of the ACDM scenario as the 
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standard paradigm in modern cosmology, although there are still some puzzling conflicts 
between ACDM predictions and current observations |[708|| . 



To the contrary of the good performance of the ACDM model, another common 
results of these works is that the DGP model and the standard CG model usually performs 
badly and ranks worst in the dark energy models considered thereof. So these two models 
faced with high crisis when confronted with observations. Besides the DGP and the CG 
model, some other models like the ADE and RDE models also performed badly in a model 
comparison. For example, in ||578|| , Li et al. obtained the following results of Xmin s fr° m 



the combination of Constitution+BAO+CMB+i^o data 

xicDM = 468.461, X rde = 493.772, xIde = 503.039, X dgp = 530.443. (15.54) 

Compared with the ACDM model, the last three models have much larger BIC values 

ABICrde = 31.308, ABIC ade = 34.578, ABICdgp = 61.982. (15.55) 

What should be mentioned is that a somewhat different result was obtained by Soller- 
man et al. in [|665[| , where they found that the flat DGP model performed even better 



than the flat ACDM model from first-year SDSS-II SNIa dataset [|330| analyzing using the 



MLSCS2k2 light-curve fitter. Notice that in this work, the authors also took two kinds 
of LTB models into consideration, which were not included in the model-comparison by 
the numerical studies listed above. In addition, they showed that the extra parameters 
required by these two models are not supported by the IC tests (their IC values are ~ 5 — 25 
larger compared with ACDM model). 

Some other models, such as the general DGP model, the HDE model and the param- 
eterizations like XCDM and CPL, can also provide rather good fits to the observational 
data. From current observational data it is hard to discriminate these models. Since they 
have more free parameters, these models are all less favored by the ACDM under the IC 
tests, indicating that given the current quality of the data there is no reason to prefer more 
complex models. The resulted ABIC values (with the ACDM model as a reference) of 
these models given in |578|j are (in a flat universe) 



ABICxcdm = 5.862, ABICgdgp = 5.897, (15.56) 
ABICcpl = 11-195, ABICrde = 8.048. (15.57) 
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Here the GDGP model refers to the generalized DGP model described by Eq. ( |15.30| ). 
Notice that the CPL model does not achieve a good performance due to its large number 
of free parameters. Another interesting phenomenon is that most of the above models 
(like the HDE model) can reduce to the ACDM model; from their best-fit parameters, it 



is found that they do tend to collapse to ACDM model [|562| , |578|1 . Therefore, the current 
observational data are still too limited to distinguish which theoretical model is better. 



16. Model-independent dark energy reconstructions 

In the last section, we have introduced some representative numerical works on the 
specific dark energy models. Due to the lack of a compelling fundamental theory to explain 
the dark energy, another route, the model-independent dark energy reconstructions, have 
drawn more and more attentions [|56W0^ , fn^fni| , |TT^ , • 



The dark energy reconstruction is a classic statistical inverse problem for the Hubble 
parameter 

H(z) = H y/n m0 (l + z^ + (l-n m0 )f(z), (16.1) 

where f(z) = Pde( z ) / Pde(®) is the dark energy density function. Different reconstruction 
method will give different f(z). An ideal dark energy reconstruction should be sufficiently 
versatile to accommodate a large class of dark energy models. The main target of a dark 
energy reconstruction is to detect the dynamical property of dark energy, i.e. determine 
whether the accelerating expansion is consistent with a cosmo logical constant. 

To begin with one should choose an appropriate quantity characterizing dark energy. 
It is widely believed that the EOS of dark energy w = Pde/ Pde holds essential clues for 
the nature of DE |718| j. It should be mentioned that w(z) is related with f(z) through an 



integration 719,720,721 



f(z) = exp ( 3 I <h 



,l + w(z') 
1 + z' 



(16.2) 



Since the ACDM model always satisfy w = — 1, the deviation from this constant EOS will 
reveal the variation of dark energy density. Therefore, most researchers have chosen to 
study dark energy by constraining w from observations. 



144 




_ p x (^) parametrization 




1 m 

i 


^t^LLLLWXLL^-^T — - J- 

(Wang *c MukJiarjeo 200-3) 



Fig. 21: A comparison between the u; parametrization and the pde parametriza- 
tion using the same observational data. The regions inside the solid and dashed 
lines correspond to la and 2a confidence regions, respectively. Clearly the pde 
reconstruction can give much tighter constraint on dark energy compared with the 
w parametrization. From [721]. 



However, in a series of works [ p5]|71^ , [7^ , [r^|7^ , [723| , [724l , Wang and collaborators 
argued that, due to the smearing effect |712fl arising from the multiple integrals relating 
w(z) to the luminosity distance of SN dz,(z), it is difficult to constrain w using the SN data 
711|| . On the contrary, since using the dark energy density pd e can minimize the smearing 



effect by removing one integral, pa e can be constrained more tightly than w given the same 
observational data (see Fig. 21 for details). It should be mentioned that there is still a 
debate on which quantity, w or pd e , is better in describing dark energy ||725j . 

In addition to w and pd e , some other quantities, such as the deceleration parame- 
ter q |[710| , |726| , [727| , |728|| , the redshift of acceleration-deceleration transition z c ||729|| , the 
statefmder diagnostic (r,s) [[73QJ , the jerk parameter j |[731| , |732|| and the diagnostic Om 
||733| , [734|| , can also provide very useful information for the study of dark energy. Although 
these quantities can accurately reconstruct some dark energy models, they have difficulty 



to discriminate between different models of dark energy ||735| . So in this section, we will 
introduce the model-independent methods based on the reconstructions of w and pa e - 
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The model-independent dark energy reconstructions can be divided into four classes: 
(i) Specific Ansatz: assuming a specific parameterized form for w(z) and estimating the 
associated parameters, (ii) Binned Parametrization: dividing the redshift range into differ- 
ent bins and using a simple local basis representation for w(z) or pd e (z)- (hi) Polynomial 
Fitting: treating the dark energy density function f(z) = Pde( z )/ P deify as a f ree function 
of redshift and representing it by using the polynomial, (iv) Gaussian Process modeling: 
using a distribution over functions that can represent w(z) and estimating the statistical 
properties thereof. These four classes of reconstruction methods and the related research 
works will be introduced in the following. 

16.1. Specific ansatz 




The "specific ansatz" is the most popular approach currently. The key idea is assuming 
a specific parameterized form for w(z) and estimating the associated parameters. A simple 
and widely used ansatz is the XCDM ansatz, in which the EOS of dark energy is a constant, 
i.e. w = const. This yields a simple form of f(z) 

f(z) = (l + z) 3{1+w) . (16.3) 
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In [j362| , by combining the WMAP5 observations with BAO and SN data, Komatsu et 
al. obtained w = — 0.992±g:ggJ at the la CL, while in |349|] , a combination analysis of 
WMAP7+BAO+SN gave w = -0.980^.^11. A more recent constraint on w by the SCP 



team [[333] also presented the consistent result (see Fig. 22). So the current observations 
still favor w = -1 (i.e. ACDM model). 
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Fig. 23: Left panel: 68.3%, 95.4%, and 99.7% confidence regions of the (wo,w a ) 
plane from Union2 SNIa sample combined with the constraints from BAO and 
CMB both with (solid contours) and without (shaded contours) systematic er- 
rors, for a flat universe. Points above the dotted line (wo + w a = 0) vio- 
late early matter domination and are implicitly disfavored in this analysis by 
the CMB and BAO data. From [333]. Right panel: Joint constraints on the 
CPL model from the WMAP7 observations. The contours show 68.3% and 
95.4% CL from WMAP+ff +SNIa (red), WMAP+BAO+# +SNIa (blue) and 
WMAP+BAO+Ho+^At+SNIa (black), for a flat universe. "Da*" denotes the 
time-delay distance to the lens system B1608+656 at z = 0.63 measured by [736]. 



From [349]. 



Besides, one can also assume that the EOS of dark energy is not a constant. The most 
popular parametrization with dynamical w, which assume w(z) = wq + w a z/(l + z), was 
firstly proposed by Chevallier and Polarski [ 737 |, then was used to explore the expansion 
history of the universe by Linder ]|738|] . So this ansatz is often called CPL parametrization. 
The corresponding f(z) is given by ||737| , [739 



f(z) = (l + z) 



3(l+w +w a ) 



exp 



3WgZ 

1 + z 



(16.4) 



Here wq denotes the value of the present EOS, while w a denotes the variation of the EOS. 
Because of its bounded behavior at high redshift and high accuracy in reconstructing many 
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scalar field EOS ||738||, the CPL parametrization has become one of the most popular 



methods to study dark energy |[714| , |740|| . In Fig. 23 we show the constraint on (wo,w a ) in 
the CPL ansatz given by the Union2 SNIa dataset p33| | and by the WMAP7 observations 
349|| . Obviously, the current data favor the result of wq = — 1 and w a = 0, which is 



consistent with the cosmological constant. 

In addition, some other ansatzs have also been proposed. Using the principal compo- 
nent analysis, Linder and Huterer [ [74 If argued that 2 parameters, involving a measure of 
the EOS value at some epoch (e.g. wq) and a measure of the change in EOS (e.g. it/), are 
most realistic in projecting dark energy parameter constraints. Therefore, most ansatzs 
contain 2 parameters. In ||742|| , Huterer and Turner proposed a linear parametrization 
w(z) = wq + w\z to study the evolution of dark energy. This ansatz can fit the low red- 
shifts data well, but its dark energy component grows increasingly unsuitable at redshifts 
z > 1. In ||714|| Jassal et al. proposed a more general form of the CPL parametrization 
w(z) = wq + w a z/(l + z) p and investigated the case of p = 1,2. In ||743|| , Efstathiou 
introduced another parametrization w(z) = wq + w\ ln(l + z). It should be mentioned 
that the parametrizations listed above have difficult to fit the rapidly varying dark energy 
models, and some other parametrizations ||713| , [744| , |745| , |746|1 have been proposed to fit a 
fast transition of w{z). For more research works of various parametrization forms, see 
|5^ , TOTOI^ , P^rara and references therein. 



16.2. Binned parametrization 

In addition to the specific ansatz, another popular approach is the binned parametriza- 
tion. The binned parametrization was firstly proposed by Huterer and Starkman [[750 



based on the principal component analysis (PCA) [[750| , |751|1 . It often used to measure the 
EOS w and the density pd e of dark energy. The key idea is dividing the redshift range 
into different bins and picking a simple local basis representation for w(z) or pde(z). The 
simplest way is setting w(z) or pde(z) as piecewise constant in redshift. For the case where 
w is piecewise constant in redshift, f(z) can be written as |752[ 

n-1 



ffa-i <z<z n ) = (l + Z )^ 1+w ^ + 



(16.5) 



i=0 



where Wi is the EOS parameter in the ith redshift bin defined by an upper boundary at 
Zi. This parametrization has been extensively studied [|733|J753| , [754|[ . For the case where 
dark energy density pd e is piecewise constant in redshift, f(z) can be written as 

1 0<z< Zl ; 
fi Zi-x < z < Zi (2 < i < n). 



;i6.e) 
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Here fa is a piecewise constant, and from the relation E(0) = 1 one can easily obtain fi = 1. 
For same number of redshift bins, the number of free parameters of piecewise constant pde 
parametrization is one fewer than that of piecewise constant w parametrization. 



N 



"(0 




Fig. 24: Uncorrelated estimates of the expansion history. Using nAz = 40, 20, 
and 15, respectively, 3, 4, or 5 bins' independent measurements of H(z) from the 
gold sample are plotted in the top panel. The solid black line in this plane denotes 
the prediction of the ACDM model with Q m o = 0.29. The bottom panel shows 
the independent measurements of the kinematic quantity a versus redshift. In this 
plane a positive or negative sign of the slope of the data indicates deceleration or 
acceleration of the expansion, respectively. From J323[| . 



It should be mentioned that the optimal choice of redshift bins is still in debate. In 
, Riess et al. proposed an uniform, unbiased binning method, in which the number 
of SNIa in each bin times the width of each bin is a constant (i.e. nAz = const). Using 
nAz = 40, 20, and 15, respectively, they derived 3, 4, or 5 bins' independent measurements 
of H(z) and d from the gold sample ||323|| . (See Fig. 24). This binning method has drawn 
a lot of attention. For examples, by setting nAz ~ 30 and using the piecewise constant w 
parametrization, Gong et al. explored the Constitution dataset in [ 755 1 , and analyzed the 
Union2 dataset in [ [746|| . 



In ||75611 ) Wang argued that one should choose a constant Az for redshift slices. This 
is because for a galaxy redshift survey, the observables are H and 1/Da (length scales 
extracted from data analysis). Since these scales are assumed to be constant in each 
redshift slice, the redshift slices should be chosen such that the variation of H and 1 / Da 
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Fig. 25: Constrains on the dark energy density from a joint data set of SN, BAO, 
CMB, and Hq. The left panel is the case for three bins, the right panel is the case 



for four bins. From [332]. 



in each redshift slice remain roughly constant with z. This binning method have also 



been adopted by some experimental groups. For example, in [[333| , the SCP SNIa group 
explored the Union2 dataset by using this constant binning method. They find that the 
current small sample of SNIa cannot constrain the existence of dark energy above redshift 
1 (see Fig. 25). 

In |[757| , |758|| , we presented a new binned parametrization method. Instead of choosing 
the discontinuity points Zj by hand, one can treat Zi as models parameters and let them 
run freely in the redshift region of SNIa samples. Using the piecewise constant w and the 
piecewise constant pd e parametrization, respectively, the Constitution SNIa dataset has 



been explored [|758|1 . In addition, utilizing the piecewise constant pa e parametrization, the 
Union2 SNIa dataset has also been analyzed [[757|| (the corresponding results are given in 
Fig. 26). These works show that the Constitution dataset favors a dynamical dark energy, 
while the Union2 dataset is still consistent with a cosmological constant. Comparing with 
those two binning methods listed above, the advantage of this binning method is that it 
can achieve much smaller Xmin^- 

Besides the piecewise constant parametrization, some other local basis representa- 



tions for w(z) or p<ie(z) are also proposed, such as wavelet [[760|| and numerical derivatives 
566| , [f6T |. 



18 



A simple comparison of these three binning methods can be seen in [759|. 
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Fig. 26: The relationship between the Xmin an d the discontinuity points of redshift 
(zi and Z2) for the 3 bins piecewise constant pd e parametrization. The left panel is 
plotted by using the Union2 SNIa sample alone, and the right panel is plotted by 
using the combined SNIa+CMB+BAO data. The x-axis represents the redshift of 
the first discontinuity point zi, while the y-axis denotes the redshift of the second 
discontinuity point Z2- Notice that the light-colored region corresponds to a big 
X 2 , and the dark-colored region corresponds to a small % 2 . Since z\ < Z2 must be 
satisfied, the bottom-right region of the figure is always blank. From [757]. 



16.3. Polynomial fitting 



The third approach is the polynomial fitting method. The key idea is treating the dark 
energy density function f(z) as a free function of redshift and representing it by using the 
polynomial. Compared with the binned parametrization, the advantage of the polynomial 
fitting parametrization is that the dark energy density function f(z) can be reconstructed 
as a continuous function in the redshift range covered by the observational data. 



A simple polynomial fit to f(z) was proposed by Alam et al. ||763|| , which is a truncated 
Taylor expansion 

f(z) =A + A 1 (l + z) + A 2 (l + z) 2 . (16.7) 

This ansatz has only three free parameters (O m o, A±, A2) since Aq + A\ + A2 = 1 — O m o 
for a flat universe. By using this ansatz, Alam et al. argued that the Tonry/Barris SNIa 
sample ||320| , |321|| appear to favour dark energy which evolves in time [[762] (see Fig. 27). 
This conclusion will be modified if the effect of the CMB/LSS observations are taken into 
account [764] . It should be mentioned that there was a debate about the reliability of this 



ansatz 765,766 
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Fig. 27: The evolution of w(z) with redshift for different values of f2 m . The left 
panel corresponds to the case of Q m = 0.2, the middle panel corresponds to the 
case of Qm = 0.3, and the right panel corresponds to the case of Q m = 0.4, In each 
panel, the thick solid line shows the best-fit, the light grey contour represents the 
ltx confidence level, and the dark grey contour represents the 2a confidence level 
around the best-fit. From (176^ 



Another interesting polynomial fit is the polynomial interpolation, which was proposed 
by Wang [P15j|345| , |7ll| , ff2rjiff2l , [72^ , ff2^^ . It choose different redshift points z { = i * 



z max /n(i = 1,2, ••-,ra), and interpolate f(z) by using its own values at these redshift 
points. This yields 

f /N f ( z ~ zi)...(z - Zi-i)(z - z i+1 )...(z - gn) , . 

^ (Zi - Z!)...{Zi ~ - Z i+1 )...( Zi - Zn) 

Here fi = f(zi) and z n = z max - Based on the relation /(0) = 1, one parameter can 
be fixed directly, and only n — 1 model parameters need to be determined by the data. 
In [720] , Wang and Tegmark made an accurate measurement of the dark energy density 



function f(z) by using the spectacular of the high redshift supernova observations from 
the HST/GOODS program and previous supernova. In ||721|| , Wang and Freese demon- 



strated that pde(z) can be constrained more tightly than w(z) given the same observa- 
tional data by using the Tonry/Barris SNIa sample ||32q , |321|| . In |[723|| , by utilizing the 



nearby+SDSS+ESSENCE+SNLS+HST set of 288 SNIa and the Constitution set of 397 
SNIa, Wang showed that flux-averaging of SNIa can be used to test the presence of un- 
known systematic uncertainties, and yield more robust distance measurements from SNIa 
(see Fig. 28 for details). The latest Union2 set of 557 SNIa have also been explored by 



using this polynomial interpolation method ||757 
16.4- Gaussian process modeling 
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Fig. 28: Dark energy density function f{z) = px (z) / px (0) measured from combin- 
ing SN la data with CMB, BAO, GRB data, and imposing the SHOES prior on Ho. 
The 68% (shaded) and 95% confidence level regions are shown. A flat universe is as- 
sumed. The left panel is plotted by using the Constitution set of 397 SNIa, while the 
right panel is plotted by using the nearby+SDSS+ESSENCE+SNLS+HST data 
set of 288 SNIa. As seen in this figure, using the nearby+SDSS+ESSENCE+SNLS 
+HST dataset gives much more stringent constraints on dark energy than using 
the Constitution dataset, and gives measurements that are closer to a cosmologi- 
cal constant. Besides, flux-averaging has larger impact on the results from using 
the Constitution set, and brings the measurements closer to that predicted by a 



cosmological constant. From [ 723 1 . 



The fourth approach is the Gaussian Process (GP) modeling, which is proposed by 
Holsclaw et al. ||767| , |768|| . GP is a stochastic process, which is indexed by z. The defining 



property of a GP is that the vector that corresponds to the process at any finite collection 
of points follows a multivariate Gaussian distribution ||769|| . GPs are elements of an infinite 



dimensional space, and can be used as the basis for a nonparametric reconstruction method. 
They are characterized by a mean and a covariance function, defined by a small number 
of hyperparameters ||769| , [770|| . 



Based on the definition of a GP, one can assume that, for any collection z%, z n , 
w(zi), w(z n ) follow a multivariate Gaussian distribution with a constant negative mean 
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Fig. 29: Reconstruction of w(z) based on GP modeling combined with MCMC. 
The left panel uses a Gaussian covariance function (a ~ 2), while the right panel 
uses an exponential covariance function (a = 1). These two results are very similar, 
both are very close and in agreement with a cosmological constant (black dashed 
line). The dark blue shaded region indicates the 68% CL, while the light blue 



region extends it to 95 % CL. From [767]. 



and exponential covariance function written as 

K(z,z') = K 2 p\ z - z '\ a . (16.9) 

The hyperparameters p G (0, 1) and k, and the parameters defining the likelihood, are 
determined by the data. The value of a G (0, 2] influences the smoothness of the GP 
realizations: for a = 2, the realizations are smooth with infinitely many derivatives, while 
a = 1 leads to rougher realizations suited to modeling continuous non-(mean-squared)- 
differentiable functions. Moreover, one can set up the following GP for w 

w(u) ~ GP(-1, K(u, it'))- (16.10) 

Making use of the Eq. ( |16.10|) , one can take advantage of the particular integral structure 
of luminosity distance di,(z) expressed by w(z) (see ||767| , [768|| for details). 

This new, nonparametric reconstruction method has the following advantages: it 
avoids artificial biases due to restricted parametric assumptions for w(z), it does not lose 
information about the data by smoothing it, and it does not introduce arbitrariness (and 
lack of error control) in reconstruction by representing the data using a certain number of 
bins, or cutting off information by using a restricted set of basis functions to represent the 
data [|767| , |768| j . In |767[j , using this reconstruction method, Holsclaw et al. reconstructed 
w(z) utilizing the Constitution dataset ||332|| . The obtained results are consistent with the 
cosmological constant, with no evidence for a systematic mean evolution in w with redshift 
(see Fig. 29). 
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17. Concluding Remarks 

We have reviewed theoretical models as well as observational technologies and exper- 
imental projects and numerical studies of dark energy. Numerous works and papers have 
been done and written since the discovery of the accelerating expansion of the universe, it 
is impossible to cover even a small part of the heroic endeavors of our community in any 
review article. 

However, the problem of understanding the nature of dark energy is as daunting as 
ever, or perhaps some already hold the key to this understanding without being commonly 
accepted yet. Clearly, there is a long long way to go for both theorists and experimentalists. 

It is without any doubt that the process of detecting the nature of dark energy and 
understanding its origin will prove to be one of the most exciting stories in modern science. 
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